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

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.

Stable Caenorhabditis elegans chromatin domains separate broadly expressed and developmentally regulated genes

Eukaryotic genomes are organized into domains of differing structure and activity. There is evidence that the domain organization of the genome regulates its activity, yet our understanding of domain properties and the factors that influence their formation is poor. Here, we use chromatin state analyses in early embryos and third-larval stage (L3) animals to investigate genome domain organization and its regulation in Caenorhabditis elegans At both stages we find that the genome is organized into extended chromatin domains of high or low gene activity defined by different subsets of states, and enriched for H3K36me3 or H3K27me3, respectively. The border regions between domains contain large intergenic regions and a high density of transcription factor binding, suggesting a role for transcription regulation in separating chromatin domains. Despite the differences in cell types, overall domain organization is remarkably similar in early embryos and L3 larvae, with conservation of 85% of domain border positions. Most genes in high-activity domains are expressed in the germ line and broadly across cell types, whereas low-activity domains are enriched for genes that are developmentally regulated. We find that domains are regulated by the germ-line H3K36 methyltransferase MES-4 and that border regions show striking remodeling of H3K27me1, supporting roles for H3K36 and H3K27 methylation in regulating domain structure. Our analyses of C. elegans chromatin domain structure show that genes are organized by type into domains that have differing modes of regulation.

A Conserved Nuclear Cyclophilin Is Required for Both RNA Polymerase II Elongation and Co-transcriptional Splicing in Caenorhabditis elegans

The elongation phase of transcription by RNA Polymerase II (Pol II) involves numerous events that are tightly coordinated, including RNA processing, histone modification, and chromatin remodeling. RNA splicing factors are associated with elongating Pol II, and the interdependent coupling of splicing and elongation has been documented in several systems. Here we identify a conserved, multi-domain cyclophilin family member, SIG-7, as an essential factor for both normal transcription elongation and co-transcriptional splicing. In embryos depleted for SIG-7, RNA levels for over a thousand zygotically expressed genes are substantially reduced, Pol II becomes significantly reduced at the 3' end of genes, marks of transcription elongation are reduced, and unspliced mRNAs accumulate. Our findings suggest that SIG-7 plays a central role in both Pol II elongation and co-transcriptional splicing and may provide an important link for their coordination and regulation.

RNA polymerase II pausing as a context-dependent reader of the genome

The RNA polymerase II (Pol II) transcribes all mRNA genes in eukaryotes and is among the most highly regulated enzymes in the cell. The classic model of mRNA gene regulation involves recruitment of the RNA polymerase to gene promoters in response to environmental signals. Higher eukaryotes have an additional ability to generate multiple cell types. This extra level of regulation enables each cell to interpret the same genome by committing to one of the many possible transcription programs and executing it in a precise and robust manner. Whereas multiple mechanisms are implicated in cell type-specific transcriptional regulation, how one genome can give rise to distinct transcriptional programs and what mechanisms activate and maintain the appropriate program in each cell remains unclear. This review focuses on the process of promoter-proximal Pol II pausing during early transcription elongation as a key step in context-dependent interpretation of the metazoan genome. We highlight aspects of promoter-proximal Pol II pausing, including its interplay with epigenetic mechanisms, that may enable cell type-specific regulation, and emphasize some of the pertinent questions that remain unanswered and open for investigation.

DNA-guided establishment of nucleosome patterns within coding regions of a eukaryotic genome

A conserved hallmark of eukaryotic chromatin architecture is the distinctive array of well-positioned nucleosomes downstream from transcription start sites (TSS). Recent studies indicate that trans-acting factors establish this stereotypical array. Here, we present the first genome-wide in vitro and in vivo nucleosome maps for the ciliate Tetrahymena thermophila. In contrast with previous studies in yeast, we find that the stereotypical nucleosome array is preserved in the in vitro reconstituted map, which is governed only by the DNA sequence preferences of nucleosomes. Remarkably, this average in vitro pattern arises from the presence of subsets of nucleosomes, rather than the whole array, in individual Tetrahymena genes. Variation in GC content contributes to the positioning of these sequence-directed nucleosomes and affects codon usage and amino acid composition in genes. Given that the AT-rich Tetrahymena genome is intrinsically unfavorable for nucleosome formation, we propose that these “seed” nucleosomes—together with trans-acting factors—may facilitate the establishment of nucleosome arrays within genes in vivo, while minimizing changes to the underlying coding sequences.

Lessons from modENCODE

The modENCODE (Model Organism Encyclopedia of DNA Elements) Consortium aimed to map functional elements-including transcripts, chromatin marks, regulatory factor binding sites, and origins of DNA replication-in the model organisms Drosophila melanogaster and Caenorhabditis elegans. During its five-year span, the consortium conducted more than 2,000 genome-wide assays in developmentally staged animals, dissected tissues, and homogeneous cell lines. Analysis of these data sets provided foundational insights into genome, epigenome, and transcriptome structure and the evolutionary turnover of regulatory pathways. These studies facilitated a comparative analysis with similar data types produced by the ENCODE Consortium for human cells. Genome organization differs drastically in these distant species, and yet quantitative relationships among chromatin state, transcription, and cotranscriptional RNA processing are deeply conserved. Of the many biological discoveries of the modENCODE Consortium, we highlight insights that emerged from integrative studies. We focus on operational and scientific lessons that may aid future projects of similar scale or aims in other, emerging model systems.

Trimethylation of Lys36 on H3 restricts gene expression change during aging and impacts life span

Pu et al. profiled the genome-wide pattern of trimethylation of Lys36 on histone 3 (H3K36me3) in the somatic cells of young and old C. elegans. Genes with dramatic expression change during aging are marked with low or undetectable levels of H3K36me3 in their gene bodies. A similar negative correlation between H3K36me3 marking and mRNA expression change during aging is also observed in Drosophila melanogaster. Inactivation of the methyltransferase met-1 resulted in a decrease in global H3K36me3 marks, an increase in mRNA expression change with age, and a shortened life span.

Functional data indicate that specific histone modification enzymes can be key to longevity in Caenorhabditis elegans, but the molecular basis of how chromatin structure modulates longevity is not well understood. In this study, we profiled the genome-wide pattern of trimethylation of Lys36 on histone 3 (H3K36me3) in the somatic cells of young and old Caenorhabditis elegans. We revealed a new role of H3K36me3 in maintaining gene expression stability through aging with important consequences on longevity. We found that genes with dramatic expression change during aging are marked with low or even undetectable levels of H3K36me3 in their gene bodies irrespective of their corresponding mRNA abundance. Interestingly, 3′ untranslated region (UTR) length strongly correlates with H3K36me3 levels and age-dependent mRNA expression stability. A similar negative correlation between H3K36me3 marking and mRNA expression change during aging was also observed in Drosophila melanogaster, suggesting a conserved mechanism for H3K36me3 in suppressing age-dependent mRNA expression change. Importantly, inactivation of the methyltransferase met-1 resulted in a decrease in global H3K36me3 marks, an increase in mRNA expression change with age, and a shortened life span, suggesting a causative role of the H3K36me3 marking in modulating age-dependent gene expression stability and longevity.

The DREAM complex promotes gene body H2A.Z for target repression

The DREAM (DP, Retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell cycle genes. Latorre et al. show that Caenorhabditis elegans DREAM targets have an unusual pattern of high gene body HTZ-1/H2A.Z. In mutants of lin-35, the sole p130/Rb-like gene in C. elegans, DREAM targets have reduced gene body HTZ-1/H2A.Z and increased expression. Consistent with a repressive role for gene body H2A.Z, many DREAM targets are up-regulated in htz-1/H2A.Z mutants.

The DREAM (DP, Retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell cycle genes, but its mechanism of action is poorly understood. Here we show that Caenorhabditis elegans DREAM targets have an unusual pattern of high gene body HTZ-1/H2A.Z. In mutants of lin-35, the sole p130/Rb-like gene in C. elegans, DREAM targets have reduced gene body HTZ-1/H2A.Z and increased expression. Consistent with a repressive role for gene body H2A.Z, many DREAM targets are up-regulated in htz-1/H2A.Z mutants. Our results indicate that the DREAM complex facilitates high gene body HTZ-1/H2A.Z, which plays a role in target gene repression.

Comparative analysis of metazoan chromatin organization

Genome function is dynamically regulated in part by chromatin, which consists of the histones, non-histone proteins and RNA molecules that package DNA. Studies in Caenorhabditis elegans and Drosophila melanogaster have contributed substantially to our understanding of molecular mechanisms of genome function in humans, and have revealed conservation of chromatin components and mechanisms. Nevertheless, the three organisms have markedly different genome sizes, chromosome architecture and gene organization. On human and fly chromosomes, for example, pericentric heterochromatin flanks single centromeres, whereas worm chromosomes have dispersed heterochromatin-like regions enriched in the distal chromosomal 'arms', and centromeres distributed along their lengths. To systematically investigate chromatin organization and associated gene regulation across species, we generated and analysed a large collection of genome-wide chromatin data sets from cell lines and developmental stages in worm, fly and human. Here we present over 800 new data sets from our ENCODE and modENCODE consortia, bringing the total to over 1,400. Comparison of combinatorial patterns of histone modifications, nuclear lamina-associated domains, organization of large-scale topological domains, chromatin environment at promoters and enhancers, nucleosome positioning, and DNA replication patterns reveals many conserved features of chromatin organization among the three organisms. We also find notable differences in the composition and locations of repressive chromatin. These data sets and analyses provide a rich resource for comparative and species-specific investigations of chromatin composition, organization and function.

Extreme HOT regions are CpG-dense promoters in C. elegans and humans

Most vertebrate promoters lie in unmethylated CpG-dense islands, whereas methylation of the more sparsely distributed CpGs in the remainder of the genome is thought to contribute to transcriptional repression. Nonmethylated CG dinucleotides are recognized by CXXC finger protein 1 (CXXC1, also known as CFP1), which recruits SETD1A (also known as Set1) methyltransferase for trimethylation of histone H3 lysine 4, an active promoter mark. Genomic regions enriched for CpGs are thought to be either absent or irrelevant in invertebrates that lack DNA methylation, such as C. elegans; however, a CXXC1 ortholog (CFP-1) is present. Here we demonstrate that C. elegans CFP-1 targets promoters with high CpG density, and these promoters are marked by high levels of H3K4me3. Furthermore, as for mammalian promoters, high CpG content is associated with nucleosome depletion irrespective of transcriptional activity. We further show that highly occupied target (HOT) regions identified by the binding of a large number of transcription factors are CpG-rich promoters in C. elegans and human genomes, suggesting that the unusually high factor association at HOT regions may be a consequence of CpG-linked chromatin accessibility. Our results indicate that nonmethylated CpG-dense sequence is a conserved genomic signal that promotes an open chromatin state, targeting by a CXXC1 ortholog, and H3K4me3 modification in both C. elegans and human genomes.

RNA polymerase II transcription elongation and Pol II CTD Ser2 phosphorylation: A tail of two kinases

The transition between initiation and productive elongation during RNA Polymerase II (Pol II) transcription is a well-appreciated point of regulation across many eukaryotes. Elongating Pol II is modified by phosphorylation of serine 2 (Ser2) on its carboxy terminal domain (CTD) by two kinases, Bur1/Ctk1 in yeast and Cdk9/Cdk12 in metazoans. Here, we discuss the roles and regulation of these kinases and their relationship to Pol II elongation control, and focus on recent data from work in C. elegans that point out gaps in our current understand of transcription elongation.

PRDE-1 is a nuclear factor essential for the biogenesis of Ruby motif-dependent piRNAs in C. elegans

Piwi-interacting RNAs (piRNA) are small regulatory RNAs with essential roles in maintaining genome integrity in animals and protists. Most Caenorhabditis elegans piRNAs are transcribed from two genomic clusters that likely contain thousands of individual transcription units; however, their biogenesis is not understood. Here we identify and characterize prde-1 (piRNA silencing-defective) as the first essential C. elegans piRNA biogenesis gene. Analysis of prde-1 provides the first direct evidence that piRNA precursors are 28- to 29-nucleotide (nt) RNAs initiating 2 nt upstream of mature piRNAs. PRDE-1 is a nuclear germline-expressed protein that localizes to chromosome IV. PRDE-1 is required specifically for the production of piRNA precursors from genomic loci containing an 8-nt upstream motif, the Ruby motif. The expression of a second class of motif-independent piRNAs is unaffected in prde-1 mutants. We exploited this finding to determine the targets of the motif-independent class of piRNAs. Together, our data provide new insights into both the biogenesis and function of piRNAs in gene regulation.

Evolutionary conservation of the eumetazoan gene regulatory landscape

Despite considerable differences in morphology and complexity of body plans among animals, a great part of the gene set is shared among Bilateria and their basally branching sister group, the Cnidaria. This suggests that the common ancestor of eumetazoans already had a highly complex gene repertoire. At present it is therefore unclear how morphological diversification is encoded in the genome. Here we address the possibility that differences in gene regulation could contribute to the large morphological divergence between cnidarians and bilaterians. To this end, we generated the first genome-wide map of gene regulatory elements in a nonbilaterian animal, the sea anemone Nematostella vectensis. Using chromatin immunoprecipitation followed by deep sequencing of five chromatin modifications and a transcriptional cofactor, we identified over 5000 enhancers in the Nematostella genome and could validate 75% of the tested enhancers in vivo. We found that in Nematostella, but not in yeast, enhancers are characterized by the same combination of histone modifications as in bilaterians, and these enhancers preferentially target developmental regulatory genes. Surprisingly, the distribution and abundance of gene regulatory elements relative to these genes are shared between Nematostella and bilaterian model organisms. Our results suggest that complex gene regulation originated at least 600 million yr ago, predating the common ancestor of eumetazoans.

Pol II docking and pausing at growth and stress genes in C. elegans

Fluctuations in nutrient availability profoundly impact gene expression. Previous work revealed postrecruitment regulation of RNA polymerase II (Pol II) during starvation and recovery in Caenorhabditis elegans, suggesting that promoter-proximal pausing promotes rapid response to feeding. To test this hypothesis, we measured Pol II elongation genome wide by two complementary approaches and analyzed elongation in conjunction with Pol II binding and expression. We confirmed bona fide pausing during starvation and also discovered Pol II docking. Pausing occurs at active stress-response genes that become downregulated in response to feeding. In contrast, "docked" Pol II accumulates without initiating upstream of inactive growth genes that become rapidly upregulated upon feeding. Beyond differences in function and expression, these two sets of genes have different core promoter motifs, suggesting alternative transcriptional machinery. Our work suggests that growth and stress genes are both regulated postrecruitment during starvation but at initiation and elongation, respectively, coordinating gene expression with nutrient availability.

Holocentromeres are dispersed point centromeres localized at transcription factor hotspots

Centromeres vary greatly in size and sequence composition, ranging from 'point' centromeres with a single cenH3-containing nucleosome to 'regional' centromeres embedded in tandemly repeated sequences to holocentromeres that extend along the length of entire chromosomes. Point centromeres are defined by sequence, whereas regional and holocentromeres are epigenetically defined by the location of cenH3-containing nucleosomes. In this study, we show that Caenorhabditis elegans holocentromeres are organized as dispersed but discretely localized point centromeres, each forming a single cenH3-containing nucleosome. These centromeric sites co-localize with kinetochore components, and their occupancy is dependent on the cenH3 loading machinery. These sites coincide with non-specific binding sites for multiple transcription factors ('HOT' sites), which become occupied when cenH3 is lost. Our results show that the point centromere is the basic unit of holocentric organization in support of the classical polycentric model for holocentromeres, and provide a mechanistic basis for understanding how centromeric chromatin might be maintained. DOI: http://dx.doi.org/10.7554/eLife.02025.001.

Holocentromeres are dispersed point centromeres localized at transcription factor hotspots

Centromeres vary greatly in size and sequence composition, ranging from ‘point’ centromeres with a single cenH3-containing nucleosome to ‘regional’ centromeres embedded in tandemly repeated sequences to holocentromeres that extend along the length of entire chromosomes. Point centromeres are defined by sequence, whereas regional and holocentromeres are epigenetically defined by the location of cenH3-containing nucleosomes. In this study, we show that Caenorhabditis elegans holocentromeres are organized as dispersed but discretely localized point centromeres, each forming a single cenH3-containing nucleosome. These centromeric sites co-localize with kinetochore components, and their occupancy is dependent on the cenH3 loading machinery. These sites coincide with non-specific binding sites for multiple transcription factors (‘HOT’ sites), which become occupied when cenH3 is lost. Our results show that the point centromere is the basic unit of holocentric organization in support of the classical polycentric model for holocentromeres, and provide a mechanistic basis for understanding how centromeric chromatin might be maintained.

DOI:http://dx.doi.org/10.7554/eLife.02025.001

During cell division, the chromosomes in the original cell must be replicated and these ‘sister chromosomes’ must then be divided equally between the two new daughter cells. At first, the sister chromosomes are held together near a region called the centromere, which is important because the microtubules that pull the sister chromosomes apart attach themselves to the centromere. In many cases, the centromere is a small region near the middle of the chromosomes, which produces a classic X shape. However, in some organisms centromeres span the entire length of the chromosomes. There are at least 13 plant and animal lineages with such holocentromeres.

Inside the nucleus of cells, DNA is wrapped around molecules called histones. There are five major families of histones, and histones belonging to one of these families—the H3 histones—are replaced by cenH3 variant histones at both conventional centromeres and holocentromeres. There are many unanswered questions about holocentromeres. In particular, do holocentromeres truly extend along the full length of the chromosomes, or are they found at a large number of specific sites?

Now Steiner and Henikoff have studied the distribution of cenH3 in the genome of the worm C. elegans to investigate holocentromeres in greater detail. These experiments showed that the holocentromere in C. elegans is actually made of about 700 individual centromeric sites distributed along the length of the chromosomes. Each of these sites contains just one nucleosome that contains cenH3, and these sites are likely to be the sites that microtubules attach to during cell division. Surprisingly, the same sites can also act as so-called ‘HOT–sites’: these sites are bound by many proteins that are involved in regulating the process by which genes are expressed as proteins, which suggests a link between centromeres and these regulatory proteins.

The work of Steiner and Henikoff describes how centromeric nucleosomes are distributed across the genome, but why and how cenH3 ends up at these particular 700 sites remains an open question.

DOI:http://dx.doi.org/10.7554/eLife.02025.002

Chromatin and DNA sequences in defining promoters for transcription initiation

One of the key events in eukaryotic gene regulation and consequent transcription is the assembly of general transcription factors and RNA polymerase II into a functional pre-initiation complex at core promoters. An emerging view of complexity arising from a variety of promoter associated DNA motifs, their binding factors and recent discoveries in characterising promoter associated chromatin properties brings an old question back into the limelight: how is a promoter defined? In addition to position-dependent DNA sequence motifs, accumulating evidence suggests that several parallel acting mechanisms are involved in orchestrating a pattern marked by the state of chromatin and general transcription factor binding in preparation for defining transcription start sites. In this review we attempt to summarise these promoter features and discuss the available evidence pointing at their interactions in defining transcription initiation in developmental contexts. This article is part of a Special Issue entitled: Chromatin and epigenetic regulation of animal development.

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

The X-chromosome gene regulatory process called dosage compensation ensures that males (1X) and females (2X) express equal levels of X-chromosome transcripts. The mechanism in Caenorhabditis elegans has been elusive due to improperly annotated transcription start sites (TSSs). Here we define TSSs and the distribution of transcriptionally engaged RNA polymerase II (Pol II) genome-wide in wild-type and dosage-compensation-defective animals to dissect this regulatory mechanism. Our TSS-mapping strategy integrates GRO-seq, which tracks nascent transcription, with a new derivative of this method, called GRO-cap, which recovers nascent RNAs with 5' caps prior to their removal by co-transcriptional processing. Our analyses reveal that promoter-proximal pausing is rare, unlike in other metazoans, and promoters are unexpectedly far upstream from the 5' ends of mature mRNAs. We find that C. elegans equalizes X-chromosome expression between the sexes, to a level equivalent to autosomes, by reducing Pol II recruitment to promoters of hermaphrodite X-linked genes using a chromosome-restructuring condensin complex. DOI:http://dx.doi.org/10.7554/eLife.00808.001.