The landscape of RNA polymerase II transcription initiation in C. elegans reveals promoter and enhancer architectures

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.

Assembling the RNA therapeutics toolbox

From the approval of COVID-19 mRNA vaccines to the 2023 Nobel Prize awarded for nucleoside base modifications, RNA therapeutics have entered the spotlight and are transforming drug development. While the term "RNA therapeutics" has been used in various contexts, this review focuses on treatments that utilize RNA as a component or target RNA for therapeutic effects. We summarize the latest advances in RNA-targeting tools and RNA-based technologies, including but not limited to mRNA, antisense oligos, siRNAs, small molecules and RNA editors. We focus on the mechanisms of current FDA-approved therapeutics but also provide a discussion on the upcoming workforces. The clinical utility of RNA-based therapeutics is enabled not only by the advances in RNA technologies but in conjunction with the significant improvements in chemical modifications and delivery platforms, which are also briefly discussed in the review. We summarize the latest RNA therapeutics based on their mechanisms and therapeutic effects, which include expressing proteins for vaccination and protein replacement therapies, degrading deleterious RNA, modulating transcription and translation efficiency, targeting noncoding RNAs, binding and modulating protein activity and editing RNA sequences and modifications. This review emphasizes the concept of an RNA therapeutic toolbox, pinpointing the readers to all the tools available for their desired research and clinical goals. As the field advances, the catalog of RNA therapeutic tools continues to grow, further allowing researchers to combine appropriate RNA technologies with suitable chemical modifications and delivery platforms to develop therapeutics tailored to their specific clinical challenges.

Genome-wide analysis of transcription-coupled repair reveals novel transcription events in Caenorhabditis elegans

Bulky DNA adducts such as those induced by ultraviolet light are removed from the genomes of multicellular organisms by nucleotide excision repair, which occurs through two distinct mechanisms, global repair, requiring the DNA damage recognition-factor XPC (xeroderma pigmentosum complementation group C), and transcription-coupled repair (TCR), which does not. TCR is initiated when elongating RNA polymerase II encounters DNA damage, and thus analysis of genome-wide excision repair in XPC-mutants only repairing by TCR provides a unique opportunity to map transcription events missed by methods dependent on capturing RNA transcription products and thus limited by their stability and/or modifications (5'-capping or 3'-polyadenylation). Here, we have performed the eXcision Repair-sequencing (XR-seq) in the model organism Caenorhabditis elegans to generate genome-wide repair maps from a wild-type strain with normal excision repair, a strain lacking TCR (csb-1), or one that only repairs by TCR (xpc-1). Analysis of the intersections between the xpc-1 XR-seq repair maps with RNA-mapping datasets (RNA-seq, long- and short-capped RNA-seq) reveal previously unrecognized sites of transcription and further enhance our understanding of the genome of this important model organism.

Quantitative analysis of transcription start site selection reveals control by DNA sequence, RNA polymerase II activity and NTP levels

Transcription start site (TSS) selection is a key step in gene expression and occurs at many promoter positions over a wide range of efficiencies. Here we develop a massively parallel reporter assay to quantitatively dissect contributions of promoter sequence, nucleoside triphosphate substrate levels and RNA polymerase II (Pol II) activity to TSS selection by 'promoter scanning' in Saccharomyces cerevisiae (Pol II MAssively Systematic Transcript End Readout, 'Pol II MASTER'). Using Pol II MASTER, we measure the efficiency of Pol II initiation at 1,000,000 individual TSS sequences in a defined promoter context. Pol II MASTER confirms proposed critical qualities of S. cerevisiae TSS -8, -1 and +1 positions, quantitatively, in a controlled promoter context. Pol II MASTER extends quantitative analysis to surrounding sequences and determines that they tune initiation over a wide range of efficiencies. These results enabled the development of a predictive model for initiation efficiency based on sequence. We show that genetic perturbation of Pol II catalytic activity alters initiation efficiency mostly independently of TSS sequence, but selectively modulates preference for the initiating nucleotide. Intriguingly, we find that Pol II initiation efficiency is directly sensitive to guanosine-5'-triphosphate levels at the first five transcript positions and to cytosine-5'-triphosphate and uridine-5'-triphosphate levels at the second position genome wide. These results suggest individual nucleoside triphosphate levels can have transcript-specific effects on initiation, representing a cryptic layer of potential regulation at the level of Pol II biochemical properties. The results establish Pol II MASTER as a method for quantitative dissection of transcription initiation in eukaryotes.

Emerging insights into enhancer biology and function

Cell type-specific gene expression is coordinated by DNA-encoded enhancers and the transcription factors (TFs) that bind to them in a sequence-specific manner. As such, these enhancers and TFs are critical mediators of normal development and altered enhancer or TF function is associated with the development of diseases such as cancer. While initially defined by their ability to activate gene transcription in reporter assays, putative enhancer elements are now frequently defined by their unique chromatin features including DNase hypersensitivity and transposase accessibility, bidirectional enhancer RNA (eRNA) transcription, CpG hypomethylation, high H3K27ac and H3K4me1, sequence-specific transcription factor binding, and co-factor recruitment. Identification of these chromatin features through sequencing-based assays has revolutionized our ability to identify enhancer elements on a genome-wide scale, and genome-wide functional assays are now capitalizing on this information to greatly expand our understanding of how enhancers function to provide spatiotemporal coordination of gene expression programs. Here, we highlight recent technological advances that are providing new insights into the molecular mechanisms by which these critical cis-regulatory elements function in gene control. We pay particular attention to advances in our understanding of enhancer transcription, enhancer-promoter syntax, 3D organization and biomolecular condensates, transcription factor and co-factor dependencies, and the development of genome-wide functional enhancer screens.

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.

Dynamic interplay between non-coding enhancer transcription and gene activity in development

Non-coding transcription at the intergenic regulatory regions is a prevalent feature of metazoan genomes, but its biological function remains uncertain. Here, we devise a live-imaging system that permits simultaneous visualization of gene activity along with intergenic non-coding transcription at single-cell resolution in Drosophila. Quantitative image analysis reveals that elongation of RNA polymerase II across the internal core region of enhancers leads to suppression of transcriptional bursting from linked genes. Super-resolution imaging and genome-editing analysis further demonstrate that enhancer transcription antagonizes molecular crowding of transcription factors, thereby interrupting the formation of a transcription hub at the gene locus. We also show that a certain class of developmental enhancers are structurally optimized to co-activate gene transcription together with non-coding transcription effectively. We suggest that enhancer function is flexibly tunable through the modulation of hub formation via surrounding non-coding transcription during development.

Enhancer-promoter entanglement explains their transcriptional interdependence

Enhancers not only activate target promoters to stimulate messenger RNA (mRNA) synthesis, but they themselves also undergo transcription to produce enhancer RNAs (eRNAs), the significance of which is not well understood. Transcription at the participating enhancer-promoter pair appears coordinated, but it is unclear why and how. Here, we employ cell-free transcription assays using constructs derived from the human GREB1 locus to demonstrate that transcription at an enhancer and its target promoter is interdependent. This interdependence is observable under conditions where direct enhancer-promoter contact (EPC) takes place. We demonstrate that transcription activation at a participating enhancer-promoter pair is dependent on i) the mutual availability of the enhancer and promoter, ii) the state of transcription at both the enhancer and promoter, iii) local abundance of both eRNA and mRNA, and iv) direct EPC. Our results suggest transcriptional interdependence between the enhancer and the promoter as the basis of their transcriptional concurrence and coordination throughout the genome. We propose a model where transcriptional concurrence, coordination and interdependence are possible if the participating enhancer and promoter are entangled in the form of EPC, reside in a proteinaceous bubble, and utilize shared transcriptional resources and regulatory inputs.

SAMS-1 coordinates HLH-30/TFEB and PHA-4/FOXA activities through histone methylation to mediate dietary restriction-induced autophagy and longevity

Dietary restriction (DR) is known to promote autophagy to exert its longevity effect. While SAMS-1 (S-adenosyl methionine synthetase-1) has been shown to be a key mediator of the DR response, little is known about the roles of S-adenosyl methionine (SAM) and SAM-dependent methyltransferase in autophagy and DR-induced longevity. In this study, we show that DR and SAMS-1 repress the activity of SET-2, a histone H3K4 methyltransferase, by limiting the availability of SAM. Consequently, the reduced H3K4me3 levels promote the expression and activity of two transcription factors, HLH-30/TFEB and PHA-4/FOXA, which both regulate the transcription of autophagy-related genes. We then find that HLH-30/TFEB and PHA-4/FOXA act collaboratively on their common target genes to mediate the transcriptional response of autophagy-related genes and consequently the lifespan of the animals. Our study thus shows that the SAMS-1-SET-2 axis serves as a nutrient-sensing module to epigenetically coordinate the activation of HLH-30/TFEB and PHA-4/FOXA transcription factors to control macroautophagy/autophagy and longevity in response to DR.Abbreviations: ChIP: chromatin immunoprecipitation; ChIP-seq: chromatin immuno precipitation-sequencing; COMPASS: complex of proteins associated with Set1; DR: dietary restriction; GO: gene ontology; SAM: S-adenosyl methionine; SAMS-1: S-adenosyl methionine synthetase-1; TSS: transcription start site; WT: wild-type.

Transcriptional enhancers at 40: evolution of a viral DNA element to nuclear architectural structures

Gene regulation by transcriptional enhancers is the dominant mechanism driving cell type- and signal-specific transcriptional diversity in metazoans. However, over four decades since the original discovery, how enhancers operate in the nuclear space remains largely enigmatic. Recent multidisciplinary efforts combining real-time imaging, genome sequencing, and biophysical strategies provide insightful but conflicting models of enhancer-mediated gene control. Here, we review the discovery and progress in enhancer biology, emphasizing the recent findings that acutely activated enhancers assemble regulatory machinery as mesoscale architectural structures with distinct physical properties. These findings help formulate novel models that explain several mysterious features of the assembly of transcriptional enhancers and the mechanisms of spatial control of gene expression.

Evolutionary Invariant of the Structure of DNA Double Helix in RNAP II Core Promoters

Eukaryotic and archaeal RNA polymerase II (POL II) machinery is highly conserved, regardless of the extreme changes in promoter sequences in different organisms. The goal of our work is to find the cause of this conservatism. The representative sets of aligned promoter sequences of fifteen organisms belonging to different evolutional stages were studied. Their textual profiles, as well as profiles of the indexes that characterize the secondary structure and the mechanical and physicochemical properties, were analyzed. The evolutionarily stable, extremely heterogeneous special secondary structure of POL II core promoters was revealed, which includes two singular regions—hexanucleotide “INR” around TSS and octanucleotide “TATA element” of about −28 bp upstream. Such structures may have developed at some stage of evolution. It turned out to be so well matched for the pre-initiation complex formation and the subsequent initiation of transcription for POL II machinery that in the course of evolution there were selected only those nucleotide sequences that were able to reproduce these structural properties. The individual features of specific sequences representing the singular region of the promoter of each gene can affect the kinetics of DNA-protein complex formation and facilitate strand separation in double-stranded DNA at the TSS position.

pop-1/TCF, ref-2/ZIC and T-box factors regulate the development of anterior cells in the C. elegans embryo

Patterning of the anterior-posterior axis is fundamental to animal development. The Wnt pathway plays a major role in this process by activating the expression of posterior genes in animals from worms to humans. This observation raises the question of whether the Wnt pathway or other regulators control the expression of the many anterior-expressed genes. We found that the expression of five anterior-specific genes in Caenorhabditis elegans embryos depends on the Wnt pathway effectors pop-1/TCF and sys-1/β-catenin. We focused further on one of these anterior genes, ref-2/ZIC, a conserved transcription factor expressed in multiple anterior lineages. Live imaging of ref-2 mutant embryos identified defects in cell division timing and position in anterior lineages. Cis-regulatory dissection identified three ref-2 transcriptional enhancers, one of which is necessary and sufficient for anterior-specific expression. This enhancer is activated by the T-box transcription factors TBX-37 and TBX-38, and surprisingly, concatemerized TBX-37/38 binding sites are sufficient to drive anterior-biased expression alone, despite the broad expression of TBX-37 and TBX-38. Taken together, our results highlight the diverse mechanisms used to regulate anterior expression patterns in the embryo.

Regulating specificity in enhancer-promoter communication

Enhancers are cis-regulatory elements that can activate transcription remotely to regulate a specific pattern of a gene's expression. Genes typically have many enhancers that are often intermingled in the loci of other genes. To regulate expression, enhancers must therefore activate their correct promoter while ignoring others that may be in closer linear proximity. In this review, we discuss mechanisms by which enhancers engage with promoters, including recent findings on the role of cohesin and the Mediator complex, and how this specificity in enhancer-promoter communication is encoded. Genetic dissection of model loci, in addition to more recent findings using genome-wide approaches, highlight the core promoter sequence, its accessibility, cofactor-promoter preference, in addition to the surrounding genomic context, as key components.

Genomic and functional conservation of lncRNAs: lessons from flies

Over the last decade, the increasing interest in long non-coding RNAs (lncRNAs) has led to the discovery of these transcripts in multiple organisms. LncRNAs tend to be specifically, and often lowly, expressed in certain tissues, cell types and biological contexts. Although lncRNAs participate in the regulation of a wide variety of biological processes, including development and disease, most of their functions and mechanisms of action remain unknown. Poor conservation of the DNA sequences encoding for these transcripts makes the identification of lncRNAs orthologues among different species very challenging, especially between evolutionarily distant species such as flies and humans or mice. However, the functions of lncRNAs are unexpectedly preserved among different species supporting the idea that conservation occurs beyond DNA sequences and reinforcing the potential of characterising lncRNAs in animal models. In this review, we describe the features and roles of lncRNAs in the fruit fly Drosophila melanogaster, focusing on genomic and functional comparisons with human and mouse lncRNAs. We also discuss the current state of advances and limitations in the study of lncRNA conservation and future perspectives.

Ssl2/TFIIH function in transcription start site scanning by RNA polymerase II in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, RNA polymerase II (Pol II) selects transcription start sites (TSSs) by a unidirectional scanning process. During scanning, a preinitiation complex (PIC) assembled at an upstream core promoter initiates at select positions within a window ~40-120 bp downstream. Several lines of evidence indicate that Ssl2, the yeast homolog of XPB and an essential and conserved subunit of the general transcription factor (GTF) TFIIH, drives scanning through its DNA-dependent ATPase activity, therefore potentially controlling both scanning rate and scanning extent (processivity). To address questions of how Ssl2 functions in promoter scanning and interacts with other initiation activities, we leveraged distinct initiation-sensitive reporters to identify novel ssl2 alleles. These ssl2 alleles, many of which alter residues conserved from yeast to human, confer either upstream or downstream TSS shifts at the model promoter ADH1 and genome-wide. Specifically, tested ssl2 alleles alter TSS selection by increasing or narrowing the distribution of TSSs used at individual promoters. Genetic interactions of ssl2 alleles with other initiation factors are consistent with ssl2 allele classes functioning through increasing or decreasing scanning processivity but not necessarily scanning rate. These alleles underpin a residue interaction network that likely modulates Ssl2 activity and TFIIH function in promoter scanning. We propose that the outcome of promoter scanning is determined by two functional networks, the first being Pol II activity and factors that modulate it to determine initiation efficiency within a scanning window, and the second being Ssl2/TFIIH and factors that modulate scanning processivity to determine the width of the scanning widow.

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.

Integrator is recruited to promoter-proximally paused RNA Pol II to generate Caenorhabditis elegans piRNA precursors

Piwi-interacting RNAs (piRNAs) play key roles in germline development and genome defence in metazoans. In C. elegans, piRNAs are transcribed from > 15,000 discrete genomic loci by RNA polymerase II (Pol II), resulting in 28 nt short-capped piRNA precursors. Here, we investigate transcription termination at piRNA loci. We show that the Integrator complex, which terminates snRNA transcription, is recruited to piRNA loci. Moreover, we demonstrate that the catalytic activity of Integrator cleaves nascent capped piRNA precursors associated with promoter-proximal Pol II, resulting in termination of transcription. Loss of Integrator activity, however, does not result in transcriptional readthrough at the majority of piRNA loci. Taken together, our results draw new parallels between snRNA and piRNA biogenesis in nematodes and provide evidence of a role for the Integrator complex as a terminator of promoter-proximal RNA polymerase II during piRNA biogenesis.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

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.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

Enhancer RNAs are an important regulatory layer of the epigenome

Noncoding RNAs (ncRNAs) direct a remarkable number of diverse functions in development and disease through their regulation of transcription, RNA processing and translation. Leading the charge in the RNA revolution is a class of ncRNAs that are synthesized at active enhancers, called enhancer RNAs (eRNAs). Here, we review recent insights into the biogenesis of eRNAs and the mechanisms underlying their multifaceted functions and consider how these findings could inform future investigations into enhancer transcription and eRNA function.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae

BACKGROUND

The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function.

RESULTS

To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model.

CONCLUSIONS

Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

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.

The full-length transcriptome of C. elegans using direct RNA sequencing

Current transcriptome annotations have largely relied on short read lengths intrinsic to the most widely used high-throughput cDNA sequencing technologies. For example, in the annotation of the Caenorhabditis elegans transcriptome, more than half of the transcript isoforms lack full-length support and instead rely on inference from short reads that do not span the full length of the isoform. We applied nanopore-based direct RNA sequencing to characterize the developmental polyadenylated transcriptome of C. elegans Taking advantage of long reads spanning the full length of mRNA transcripts, we provide support for 23,865 splice isoforms across 14,611 genes, without the need for computational reconstruction of gene models. Of the isoforms identified, 3452 are novel splice isoforms not present in the WormBase WS265 annotation. Furthermore, we identified 16,342 isoforms in the 3' untranslated region (3' UTR), 2640 of which are novel and do not fall within 10 bp of existing 3'-UTR data sets and annotations. Combining 3' UTRs and splice isoforms, we identified 28,858 full-length transcript isoforms. We also determined that poly(A) tail lengths of transcripts vary across development, as do the strengths of previously reported correlations between poly(A) tail length and expression level, and poly(A) tail length and 3'-UTR length. Finally, we have formatted this data as a publicly accessible track hub, enabling researchers to explore this data set easily in a genome browser.

ORSO (Online Resource for Social Omics): A data-driven social network connecting scientists to genomics datasets

High-throughput sequencing has become ubiquitous in biomedical sciences. As new technologies emerge and sequencing costs decline, the diversity and volume of available data increases exponentially, and successfully navigating the data becomes more challenging. Though datasets are often hosted by public repositories, scientists must rely on inconsistent annotation to identify and interpret meaningful data. Moreover, the experimental heterogeneity and wide-ranging quality of high-throughput biological data means that even data with desired cell lines, tissue types, or molecular targets may not be readily interpretable or integrated. We have developed ORSO (Online Resource for Social Omics) as an easy-to-use web application to connect life scientists with genomics data. In ORSO, users interact within a data-driven social network, where they can favorite datasets and follow other users. In addition to more than 30,000 datasets hosted from major biomedical consortia, users may contribute their own data to ORSO, facilitating its discovery by other users. Leveraging user interactions, ORSO provides a novel recommendation system to automatically connect users with hosted data. In addition to social interactions, the recommendation system considers primary read coverage information and annotated metadata. Similarities used by the recommendation system are presented by ORSO in a graph display, allowing exploration of dataset associations. The topology of the network graph reflects established biology, with samples from related systems grouped together. We tested the recommendation system using an RNA-seq time course dataset from differentiation of embryonic stem cells to cardiomyocytes. The ORSO recommendation system correctly predicted early data point sources as embryonic stem cells and late data point sources as heart and muscle samples, resulting in recommendation of related datasets. By connecting scientists with relevant data, ORSO provides a critical new service that facilitates wide-ranging research interests.

New sequencing technologies have rapidly transformed biomedical research. Public data repositories now contain millions of datasets, which have the potential to accelerate and bolster research projects. However, the sheer magnitude of available data makes navigation difficult. We created ORSO (Online Resource for Social Omics) to address these challenges. ORSO is a social network where entries are not status updates or tweets, but biological datasets. Users may add their own data to ORSO, joining 30,000 validated datasets that are already hosted, and other users may find these data through intuitive search functions and informative analytics. Users can then favorite datasets relevant to their interests or follow contributing users. ORSO also uses a recommendation system like those used on commercial websites to automatically recommend data to users based on user interactions and dataset similarities. By making data more accessible and by connecting users to relevant data, we anticipate that ORSO will be an important resource for scientists. ORSO may be the first of many applications that use methods originating in social media and ecommerce to enhance and further research projects in the life sciences.

DAF-16/FOXO requires Protein Phosphatase 4 to initiate transcription of stress resistance and longevity promoting genes

In C. elegans, the conserved transcription factor DAF-16/FOXO is a powerful aging regulator, relaying dire conditions into expression of stress resistance and longevity promoting genes. For some of these functions, including low insulin/IGF signaling (IIS), DAF-16 depends on the protein SMK-1/SMEK, but how SMK-1 exerts this role has remained unknown. We show that SMK-1 functions as part of a specific Protein Phosphatase 4 complex (PP4^SMK-1). Loss of PP4^SMK-1 hinders transcriptional initiation at several DAF-16-activated genes, predominantly by impairing RNA polymerase II recruitment to their promoters. Search for the relevant substrate of PP4^SMK-1 by phosphoproteomics identified the conserved transcriptional regulator SPT-5/SUPT5H, whose knockdown phenocopies the loss of PP4^SMK-1. Phosphoregulation of SPT-5 is known to control transcriptional events such as elongation and termination. Here we also show that transcription initiating events are influenced by the phosphorylation status of SPT-5, particularly at DAF-16 target genes where transcriptional initiation appears rate limiting, rendering PP4^SMK-1 crucial for many of DAF-16’s physiological roles.

The transcription factor DAF-16/FOXO mediates a wide variety of aging-preventive responses by driving the expression of stress resistance and longevity promoting genes. Here the authors show that transcriptional initiation at many DAF-16/FOXO target genes requires the dephosphorylation of SPT-5 by Protein Phosphatase 4.

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 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.

Transcriptional control by enhancers and enhancer RNAs

The regulation of gene expression is a fundamental cellular process and its misregulation is a key component of disease. Enhancers are one of the most salient regulatory elements in the genome and help orchestrate proper spatiotemporal gene expression during development, in homeostasis, and in response to signaling. Notably, molecular aberrations at enhancers, such as translocations and single nucleotide polymorphisms, are emerging as an important source of human variation and susceptibility to disease. Herein we discuss emerging paradigms addressing how genes are regulated by enhancers, common features of active enhancers, and how non-coding enhancer RNAs (eRNAs) can direct gene expression programs that underlie cellular phenotypes. We survey the current evidence, which suggests that eRNAs can bind to transcription factors, mediate enhancer-promoter interactions, influence RNA Pol II elongation, and act as decoys for repressive cofactors. Furthermore, we discuss current methodologies for the identification of eRNAs and novel approaches to elucidate their functions.

Enhancer RNAs: Insights Into Their Biological Role

Enhancers play a central role in the transcriptional regulation of metazoans. Almost a decade ago, the discovery of their pervasive transcription into noncoding RNAs, termed enhancer RNAs (eRNAs), opened a whole new field of study. The presence of eRNAs correlates with enhancer activity; however, whether they act as functional molecules remains controversial. Here we review direct experimental evidence supporting a functional role of eRNAs in transcription and provide a general pipeline that could help in the design of experimental approaches to investigate the function of eRNAs. We propose that induction of transcriptional activity at enhancers promotes an increase in its activity by an RNA-mediated titration of regulatory proteins that can impact different processes like chromatin accessibility or chromatin looping. In a few cases, transcripts originating from enhancers have acquired specific molecular functions to regulate gene expression. We speculate that these transcripts are either nonannotated long noncoding RNAs (lncRNAs) or are evolving toward functional lncRNAs. Further work will be needed to comprehend better the biological activity of these transcripts.

Identification of functional long non-coding RNAs in C. elegans

BACKGROUND

Functional characterisation of the compact genome of the model organism Caenorhabditis elegans remains incomplete despite its sequencing 20 years ago. The last decade of research has seen a tremendous increase in the number of non-coding RNAs identified in various organisms. While we have mechanistic understandings of small non-coding RNA pathways, long non-coding RNAs represent a diverse class of active transcripts whose function remains less well characterised.

RESULTS

By analysing hundreds of published transcriptome datasets, we annotated 3392 potential lncRNAs including 143 multi-exonic loci that showed increased nucleotide conservation and GC content relative to other non-coding regions. Using CRISPR/Cas9 genome editing, we generated deletion mutants for ten long non-coding RNA loci. Using automated microscopy for in-depth phenotyping, we show that six of the long non-coding RNA loci are required for normal development and fertility. Using RNA interference-mediated gene knock-down, we provide evidence that for two of the long non-coding RNA loci, the observed phenotypes are dependent on the corresponding RNA transcripts.

CONCLUSIONS

Our results highlight that a large section of the non-coding regions of the C. elegans genome remains unexplored. Based on our in vivo analysis of a selection of high-confidence lncRNA loci, we expect that a significant proportion of these high-confidence regions is likely to have a biological function at either the genomic or the transcript level.

Comparative Epigenomics Reveals that RNA Polymerase II Pausing and Chromatin Domain Organization Control Nematode piRNA Biogenesis

Piwi-interacting RNAs (piRNAs) are important for genome regulation across metazoans, but their biogenesis evolves rapidly. In Caenorhabditis elegans, piRNA loci are clustered within two 3-Mb regions on chromosome IV. Each piRNA locus possesses an upstream motif that recruits RNA polymerase II to produce an ∼28 nt primary transcript. We used comparative epigenomics across nematodes to gain insight into the origin, evolution, and mechanism of nematode piRNA biogenesis. We show that the piRNA upstream motif is derived from core promoter elements controlling snRNA transcription. We describe two alternative modes of piRNA organization in nematodes: in C. elegans and closely related nematodes, piRNAs are clustered within repressive H3K27me3 chromatin, while in other species, typified by Pristionchus pacificus, piRNAs are found within introns of active genes. Additionally, we discover that piRNA production depends on sequence signals associated with RNA polymerase II pausing. We show that pausing signals synergize with chromatin to control piRNA transcription.

•

Nematode piRNA transcription evolved from small nuclear RNA biogenesis

•

Clustered piRNAs are produced from regulated (H3K27me3) chromatin domains

•

Dispersed piRNAs are produced from active (H3K36me3) chromatin domains

•

RNA polymerase II pausing determines the short (∼28 nt) length of piRNA precursors

Highly efficient 5' capping of mitochondrial RNA with NAD+ and NADH by yeast and human mitochondrial RNA polymerase

Bacterial and eukaryotic nuclear RNA polymerases (RNAPs) cap RNA with the oxidized and reduced forms of the metabolic effector nicotinamide adenine dinucleotide, NAD⁺ and NADH, using NAD⁺ and NADH as non-canonical initiating nucleotides for transcription initiation. Here, we show that mitochondrial RNAPs (mtRNAPs) cap RNA with NAD⁺ and NADH, and do so more efficiently than nuclear RNAPs. Direct quantitation of NAD⁺- and NADH-capped RNA demonstrates remarkably high levels of capping in vivo: up to ~60% NAD⁺ and NADH capping of yeast mitochondrial transcripts, and up to ~15% NAD⁺ capping of human mitochondrial transcripts. The capping efficiency is determined by promoter sequence at, and upstream of, the transcription start site and, in yeast and human cells, by intracellular NAD⁺ and NADH levels. Our findings indicate mtRNAPs serve as both sensors and actuators in coupling cellular metabolism to mitochondrial transcriptional outputs, sensing NAD⁺ and NADH levels and adjusting transcriptional outputs accordingly.

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.

Eukaryotic core promoters and the functional basis of transcription initiation

RNA polymerase II (Pol II) core promoters are specialized DNA sequences at transcription start sites of protein-coding and non-coding genes that support the assembly of the transcription machinery and transcription initiation. They enable the highly regulated transcription of genes by selectively integrating regulatory cues from distal enhancers and their associated regulatory proteins. In this Review, we discuss the defining properties of gene core promoters, including their sequence features, chromatin architecture and transcription initiation patterns. We provide an overview of molecular mechanisms underlying the function and regulation of core promoters and their emerging functional diversity, which defines distinct transcription programmes. On the basis of the established properties of gene core promoters, we discuss transcription start sites within enhancers and integrate recent results obtained from dedicated functional assays to propose a functional model of transcription initiation. This model can explain the nature and function of transcription initiation at gene starts and at enhancers and can explain the different roles of core promoters, of Pol II and its associated factors and of the activating cues provided by enhancers and the transcription factors and cofactors they recruit.

Young genes have distinct gene structure, epigenetic profiles, and transcriptional regulation

Species-specific, new, or "orphan" genes account for 10%-30% of eukaryotic genomes. Although initially considered to have limited function, an increasing number of orphan genes have been shown to provide important phenotypic innovation. How new genes acquire regulatory sequences for proper temporal and spatial expression is unknown. Orphan gene regulation may rely in part on origination in open chromatin adjacent to preexisting promoters, although this has not yet been assessed by genome-wide analysis of chromatin states. Here, we combine taxon-rich nematode phylogenies with Iso-Seq, RNA-seq, ChIP-seq, and ATAC-seq to identify the gene structure and epigenetic signature of orphan genes in the satellite model nematode Pristionchus pacificus Consistent with previous findings, we find young genes are shorter, contain fewer exons, and are on average less strongly expressed than older genes. However, the subset of orphan genes that are expressed exhibit distinct chromatin states from similarly expressed conserved genes. Orphan gene transcription is determined by a lack of repressive histone modifications, confirming long-held hypotheses that open chromatin is important for new gene formation. Yet orphan gene start sites more closely resemble enhancers defined by H3K4me1, H3K27ac, and ATAC-seq peaks, in contrast to conserved genes that exhibit traditional promoters defined by H3K4me3 and H3K27ac. Although the majority of orphan genes are located on chromosome arms that contain high recombination rates and repressive histone marks, strongly expressed orphan genes are more randomly distributed. Our results support a model of new gene origination by rare integration into open chromatin near enhancers.

A large transcribed enhancer region regulates C. elegans bed-3 and the development of egg laying muscles

Gene expression is regulated by the interaction of the RNA polymerase with various transcription factors at promoter and enhancer elements. Transcriptome analyses found that many non-protein-coding regions are transcribed to produce long non-coding RNAs and enhancer-associated RNAs. Production of these transcripts is associated with activation of nearby protein-coding genes, and at least in some cases, the transcripts themselves mediate this activation. Non-coding transcripts are also reported from large enhancers or clusters of enhancers. However, not much is known about the function of large transcribed enhancer regions during organismal development. Here we investigated a transcribed 10.6 kb intergenic region located upstream of the C. elegans bed-3 gene. We found that parts of this region exhibit tissue-specific promoter and enhancer activities. Deletion of the region disrupts egg laying, a phenotype also observed in bed-3 mutants, but with the severity correlating with the size of the deletion. This phenotype is not caused by overall reduction in bed-3 expression. Rather, deletions reduce bed-3 expression specifically in the mesoderm lineage. We found that bed-3 has a previously unknown function in the generation of sex myoblast (SM) cells from the M lineage, and deletions cause loss of SM cells leading to loss of vulval muscles required for egg laying. Furthermore, injection of dsRNA targeting non-coding transcripts from this region disrupted egg laying in the wild type but not in RNAi-defective mutants. Therefore, the region upstream of bed-3 is required for robust expression of bed-3 in a specific tissue, and non-coding transcripts may mediate this interaction.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

The degree of enhancer or promoter activity is reflected by the levels and directionality of eRNA transcription

Here, Mikhaylichenko et al. investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. The authors demonstrate that while the timing of enhancer transcription is correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers and conclude that this is likely an inherent sequence property of the elements themselves.

Gene expression is regulated by promoters, which initiate transcription, and enhancers, which control their temporal and spatial activity. However, the discovery that mammalian enhancers also initiate transcription questions the inherent differences between enhancers and promoters. Here, we investigate the transcriptional properties of enhancers during Drosophila embryogenesis using characterized developmental enhancers. We show that while the timing of enhancer transcription is generally correlated with enhancer activity, the levels and directionality of transcription are highly varied among active enhancers. To assess how this impacts function, we developed a dual transgenic assay to simultaneously measure enhancer and promoter activities from a single element in the same embryo. Extensive transgenic analysis revealed a relationship between the direction of endogenous transcription and the ability to function as an enhancer or promoter in vivo, although enhancer RNA (eRNA) production and activity are not always strictly coupled. Some enhancers (mainly bidirectional) can act as weak promoters, producing overlapping spatio–temporal expression. Conversely, bidirectional promoters often act as strong enhancers, while unidirectional promoters generally cannot. The balance between enhancer and promoter activity is generally reflected in the levels and directionality of eRNA transcription and is likely an inherent sequence property of the elements themselves.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

Chromatin accessibility dynamics reveal novel functional enhancers in C. elegans

Chromatin accessibility, a crucial component of genome regulation, has primarily been studied in homogeneous and simple systems, such as isolated cell populations or early-development models. Whether chromatin accessibility can be assessed in complex, dynamic systems in vivo with high sensitivity remains largely unexplored. In this study, we use ATAC-seq to identify chromatin accessibility changes in a whole animal, the model organism Caenorhabditis elegans, from embryogenesis to adulthood. Chromatin accessibility changes between developmental stages are highly reproducible, recapitulate histone modification changes, and reveal key regulatory aspects of the epigenomic landscape throughout organismal development. We find that over 5000 distal noncoding regions exhibit dynamic changes in chromatin accessibility between developmental stages and could thereby represent putative enhancers. When tested in vivo, several of these putative enhancers indeed drive novel cell-type- and temporal-specific patterns of expression. Finally, by integrating transcription factor binding motifs in a machine learning framework, we identify EOR-1 as a unique transcription factor that may regulate chromatin dynamics during development. Our study provides a unique resource for C. elegans, a system in which the prevalence and importance of enhancers remains poorly characterized, and demonstrates the power of using whole organism chromatin accessibility to identify novel regulatory regions in complex systems.

Genome-wide discovery of active regulatory elements and transcription factor footprints in Caenorhabditis elegans using DNase-seq

Deep sequencing of size-selected DNase I–treated chromatin (DNase-seq) allows high-resolution measurement of chromatin accessibility to DNase I cleavage, permitting identification of de novo active cis-regulatory modules (CRMs) and individual transcription factor (TF) binding sites. We adapted DNase-seq to nuclei isolated from C. elegans embryos and L1 arrest larvae to generate high-resolution maps of TF binding. Over half of embryonic DNase I hypersensitive sites (DHSs) were annotated as noncoding, with 24% in intergenic, 12% in promoters, and 28% in introns, with similar statistics observed in L1 arrest larvae. Noncoding DHSs are highly conserved and enriched in marks of enhancer activity and transcription. We validated noncoding DHSs against known enhancers from myo-2, myo-3, hlh-1, elt-2, and lin-26/lir-1 and recapitulated 15 of 17 known enhancers. We then mined DNase-seq data to identify putative active CRMs and TF footprints. Using DNase-seq data improved predictions of tissue-specific expression compared with motifs alone. In a pilot functional test, 10 of 15 DHSs from pha-4, icl-1, and ceh-13 drove reporter gene expression in transgenic C. elegans. Overall, we provide experimental annotation of 26,644 putative CRMs in the embryo containing 55,890 TF footprints, as well as 15,841 putative CRMs in the L1 arrest larvae containing 32,685 TF footprints.

Stable Heritable Germline Silencing Directs Somatic Silencing at an Endogenous Locus

The importance of transgenerationally inherited epigenetic states to organismal fitness remains unknown as well-documented examples are often not amenable to mechanistic analysis or rely on artificial reporter loci. Here we describe an induced silenced state at an endogenous locus that persists, at 100% transmission without selection, for up to 13 generations. This unusually persistent silencing enables a detailed molecular genetic analysis of an inherited epigenetic state. We find that silencing is dependent on germline nuclear RNAi factors and post-transcriptional mechanisms. Consistent with these later observations, inheritance does not require the silenced locus, and we provide genetic evidence that small RNAs embody the inherited silencing signal. Notably, heritable germline silencing directs somatic epigenetic silencing. Somatic silencing does not require somatic nuclear RNAi but instead requires both maternal germline nuclear RNAi and chromatin-modifying activity. Coupling inherited germline silencing to somatic silencing may enable selection for physiologically important traits.

An Alternative STAT Signaling Pathway Acts in Viral Immunity in Caenorhabditis elegans

Across metazoans, innate immunity is vital in defending organisms against viral infection. In mammals, antiviral innate immunity is orchestrated by interferon signaling, activating the STAT transcription factors downstream of the JAK kinases to induce expression of antiviral effector genes. In the nematode Caenorhabditis elegans, which lacks the interferon system, the major antiviral response so far described is RNA interference (RNAi), but whether additional gene expression responses are employed is not known. Here we show that, despite the absence of both interferon and JAK, the C. elegans STAT homolog STA-1 orchestrates antiviral immunity. Intriguingly, mutants lacking STA-1 are less permissive to antiviral infection. Using gene expression analysis and chromatin immunoprecipitation, we show that, in contrast to the mammalian pathway, STA-1 acts mostly as a transcriptional repressor. Thus, STA-1 might act to suppress a constitutive antiviral response in the absence of infection. Additionally, using a reverse genetic screen, we identify the kinase SID-3 as a new component of the response to infection, which, along with STA-1, participates in the transcriptional regulatory network of the immune response. Our work uncovers novel physiological roles for two factors in viral infection: a SID protein acting independently of RNAi and a STAT protein acting in C. elegans antiviral immunity. Together, these results illustrate the complex evolutionary trajectory displayed by innate immune signaling pathways across metazoan organisms.IMPORTANCE Since innate immunity was discovered, a diversity of pathways has arisen as powerful first-line defense mechanisms to fight viral infection. RNA interference, reported mostly in invertebrates and plants, as well as the mammalian interferon response and JAK/STAT pathway are key in RNA virus innate immunity. We studied infection by the Orsay virus in Caenorhabditis elegans, where RNAi is known to be a potent antiviral defense. We show that, in addition to its RNAi pathway, C. elegans utilizes an alternative STAT pathway to control the levels of viral infection. We identify the transcription factor STA-1 and the kinase SID-3 as two components of this response. Our study defines C. elegans as a new example of the diversity of antiviral strategies.

Temporal regulation of epithelium formation mediated by FoxA, MKLP1, MgcRacGAP, and PAR-6

During embryo morphogenesis, minor epithelia are generated after, and then form bridges between, major epithelia (e.g., epidermis and gut). In Caenorhabditis elegans, this delay is regulated by four proteins that control production and localization of polarity proteins: the pioneer factor PHA-4/FoxA, kinesin ZEN-4/MKLP1, its partner CYK-4/MgcRacGAP, and PAR-6.

To establish the animal body plan, embryos link the external epidermis to the internal digestive tract. In Caenorhabditis elegans, this linkage is achieved by the arcade cells, which form an epithelial bridge between the foregut and epidermis, but little is known about how development of these three epithelia is coordinated temporally. The arcade cell epithelium is generated after the epidermis and digestive tract epithelia have matured, ensuring that both organs can withstand the mechanical stress of embryo elongation; mistiming of epithelium formation leads to defects in morphogenesis. Using a combination of genetic, bioinformatic, and imaging approaches, we find that temporal regulation of the arcade cell epithelium is mediated by the pioneer transcription factor and master regulator PHA-4/FoxA, followed by the cytoskeletal regulator and kinesin ZEN-4/MKLP1 and the polarity protein PAR-6. We show that PHA-4 directly activates mRNA expression of a broad cohort of epithelial genes, including junctional factor dlg-1. Accumulation of DLG-1 protein is delayed by ZEN-4, acting in concert with its binding partner CYK-4/MgcRacGAP. Our structure–function analysis suggests that nuclear and kinesin functions are dispensable, whereas binding to CYK-4 is essential, for ZEN-4 function in polarity. Finally, PAR-6 is necessary to localize polarity proteins such as DLG-1 within adherens junctions and at the apical surface, thereby generating arcade cell polarity. Our results reveal that the timing of a landmark event during embryonic morphogenesis is mediated by the concerted action of four proteins that delay the formation of an epithelial bridge until the appropriate time. In addition, we find that mammalian FoxA associates with many epithelial genes, suggesting that direct regulation of epithelial identity may be a conserved feature of FoxA factors and a contributor to FoxA function in development and cancer.

Landscape of histone modifications in a sponge reveals the origin of animal cis-regulatory complexity

Combinatorial patterns of histone modifications regulate developmental and cell type-specific gene expression and underpin animal complexity, but it is unclear when this regulatory system evolved. By analysing histone modifications in a morphologically-simple, early branching animal, the sponge Amphimedonqueenslandica, we show that the regulatory landscape used by complex bilaterians was already in place at the dawn of animal multicellularity. This includes distal enhancers, repressive chromatin and transcriptional units marked by H3K4me3 that vary with levels of developmental regulation. Strikingly, Amphimedon enhancers are enriched in metazoan-specific microsyntenic units, suggesting that their genomic location is extremely ancient and likely to place constraints on the evolution of surrounding genes. These results suggest that the regulatory foundation for spatiotemporal gene expression evolved prior to the divergence of sponges and eumetazoans, and was necessary for the evolution of animal multicellularity.

The origins of developmental gene regulation

The leap from simple unicellularity to complex multicellularity remains one of life's major enigmas. The origins of metazoan developmental gene regulatory mechanisms are sought by analyzing gene regulation in extant eumetazoans, sponges, and unicellular organisms. The main hypothesis of this manuscript is that, developmental enhancers evolved from unicellular inducible promoters that diversified the expression of regulatory genes during metazoan evolution. Promoters and enhancers are functionally similar; both can regulate the transcription of distal promoters and both direct local transcription. Additionally, enhancers have experimentally characterized structural features that reveal their origin from inducible promoters. The distal co-operative regulation among promoters identified in unicellular opisthokonts possibly represents the precursor of distal regulation of promoters by enhancers. During metazoan evolution, constitutive-type promoters of regulatory genes would have acquired novel receptivity to distal regulatory inputs from promoters of inducible genes that eventually specialized as enhancers. The novel regulatory interactions would have caused constitutively expressed genes controlling differential gene expression in unicellular organisms to become themselves differentially expressed. The consequence of the novel regulatory interactions was that regulatory pathways of unicellular organisms became interlaced and ultimately evolved into the intricate developmental gene regulatory networks (GRNs) of extant metazoans.

The gastrula transition reorganizes replication-origin selection in Caenorhabditis elegans

Although some features underlying replication-origin activation in metazoan cells have been determined, little is known about their regulation during metazoan development. Using the nascent-strand purification method, here we identified replication origins throughout Caenorhabditis elegans embryonic development and found that the origin repertoire is thoroughly reorganized after gastrulation onset. During the pluripotent embryonic stages (pregastrula), potential cruciform structures and open chromatin are determining factors that establish replication origins. The observed enrichment of replication origins in transcription factor-binding sites and their presence in promoters of highly transcribed genes, particularly operons, suggest that transcriptional activity contributes to replication initiation before gastrulation. After the gastrula transition, when embryonic differentiation programs are set, new origins are selected at enhancers, close to CpG-island-like sequences, and at noncoding genes. Our findings suggest that origin selection coordinates replication initiation with transcriptional programs during metazoan development.

SeqPlots - Interactive software for exploratory data analyses, pattern discovery and visualization in genomics

Experiments involving high-throughput sequencing are widely used for analyses of chromatin function and gene expression. Common examples are the use of chromatin immunoprecipitation for the analysis of chromatin modifications or factor binding, enzymatic digestions for chromatin structure assays, and RNA sequencing to assess gene expression changes after biological perturbations. To investigate the pattern and abundance of coverage signals across regions of interest, data are often visualized as profile plots of average signal or stacked rows of signal in the form of heatmaps. We found that available plotting software was either slow and laborious or difficult to use by investigators with little computational training, which inhibited wide data exploration. To address this need, we developed SeqPlots, a user-friendly exploratory data analysis (EDA) and visualization software for genomics. After choosing groups of signal and feature files and defining plotting parameters, users can generate profile plots of average signal or heatmaps clustered using different algorithms in a matter of seconds through the graphical user interface (GUI) controls. SeqPlots accepts all major genomic file formats as input and can also generate and plot user defined motif densities. Profile plots and heatmaps are highly configurable and batch operations can be used to generate a large number of plots at once. SeqPlots is available as a GUI application for Mac or Windows and Linux, or as an R/Bioconductor package. It can also be deployed on a server for remote and collaborative usage. The analysis features and ease of use of SeqPlots encourages wide data exploration, which should aid the discovery of novel genomic associations.