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

Chromatin and epigenetics in aging biology

This book chapter will focus on modifications to chromatin itself, how chromatin modifications are regulated, and how these modifications are deciphered by the cell to impact aging. In this chapter, we will review how chromatin modifications change with age, examine how chromatin-modifying enzymes have been shown to regulate aging and healthspan, discuss how some of these epigenetic changes are triggered and how they can regulate the lifespan of the individual and its naïve descendants, and speculate on future directions for the field.

Cohesin organizes 3D DNA contacts surrounding active enhancers in C. elegans

In mammals, cohesin and CTCF organize the 3D genome into topologically associating domains (TADs) to regulate communication between cis-regulatory elements. Many organisms, including S. cerevisiae, C. elegans, and A. thaliana contain cohesin but lack CTCF. Here, we used C. elegans to investigate the function of cohesin in 3D genome organization in the absence of CTCF. Using Hi-C data, we observe cohesin-dependent features called "fountains," which have also been reported in zebrafish and mice. These are population average reflections of DNA loops originating from distinct genomic regions and are ∼20-40 kb in C. elegans Hi-C analysis upon cohesin and WAPL-1 depletion supports the idea that cohesin is preferentially loaded at sites bound by the C. elegans ortholog of NIPBL and loop extrudes in an effectively two-sided manner. ChIP-seq analyses show that cohesin translocation along the fountain trajectory depends on a fully intact complex and is extended upon WAPL-1 depletion. Hi-C contact patterns at individual fountains suggest that cohesin processivity is unequal on each side, possibly owing to collision with cohesin loaded from surrounding sites. The putative cohesin loading sites are closest to active enhancers, and fountain strength is associated with transcription. Compared with mammals, the average processivity of C. elegans cohesin is about 10-fold shorter, and the binding of NIPBL ortholog does not depend on cohesin. We propose that preferential loading and loop extrusion by cohesin is an evolutionarily conserved mechanism that regulates the 3D interactions of enhancers in animal genomes.

CelEst: a unified gene regulatory network for estimating transcription factor activities in C. elegans

Transcription factors (TFs) play a pivotal role in orchestrating critical intricate patterns of gene regulation. Although gene expression is complex, differential expression of hundreds of genes is often due to regulation by just a handful of TFs. Despite extensive efforts to elucidate TF-target regulatory relationships in Caenorhabditis elegans, existing experimental datasets cover distinct subsets of TFs and leave data integration challenging. Here, I introduce CelEst, a unified gene regulatory network designed to estimate the activity of 487 distinct C. elegans TFs-∼58% of the total-from gene expression data. To integrate data from ChIP-seq, DNA-binding motifs, and eY1H screens, optimal processing of each data type was benchmarked against a set of TF perturbation RNA-seq experiments. Moreover, I showcase how leveraging TF motif conservation in target promoters across genomes of related species can distinguish highly informative interactions, a strategy which can be applied to many model organisms. Integrated analyses of data from commonly studied conditions including heat shock, bacterial infection, and sex differences validates CelEst's performance and highlights overlooked TFs that likely play major roles in coordinating the transcriptional response to these conditions. CelEst can infer TF activity on a standard laptop computer within minutes. Furthermore, an R Shiny app with a step-by-step guide is provided for the community to perform rapid analysis with minimal coding required. I anticipate that widespread adoption of CelEsT will significantly enhance the interpretive power of transcriptomic experiments, both present and retrospective, thereby advancing our understanding of gene regulation in C. elegans and beyond.

Inferring the Selective History of CNVs Using a Maximum Likelihood Model

Copy number variations (CNVs)-structural variations generated by deletion and/or duplication that result in a change in DNA dosage-are prevalent in nature. CNVs can drastically affect the phenotype of an organism and have been shown to be both involved in genetic disorders and be used as raw material in adaptive evolution. Unlike single-nucleotide variations, the often large and varied effects of CNVs on phenotype hinders our ability to infer their selective advantage based on the population genetics data. Here, we present a likelihood-based approach, dubbed PoMoCNV (POlymorphism-aware phylogenetic MOdel for CNVs), that estimates the evolutionary parameters such as mutation rates among different copy numbers and relative fitness loss per copy deletion at a genomic locus based on population genetics data. As a case study, we analyze the genomics data of 40 strains of Caenorhabditis elegans, representing four different populations. We take advantage of the data on chromatin accessibility to interpret the mutation rate and fitness of copy numbers, as inferred by PoMoCNV, specifically in open or closed chromatin loci. We further test the reliability of PoMoCNV by estimating the evolutionary parameters of CNVs for mutation-accumulation experiments in C. elegans with varying levels of genetic drift.

Gene expression variation across genetically identical individuals predicts reproductive traits

In recent decades, genome-wide association studies (GWAS) have been the major approach to understand the biological basis of individual differences in traits and diseases. However, GWAS approaches have limited predictive power to explain individual differences, particularly for complex traits and diseases in which environmental factors play a substantial role in their etiology. Indeed, individual differences persist even in genetically identical individuals, although fully separating genetic and environmental causation is difficult in most organisms. To understand the basis of individual differences in the absence of genetic differences, we measured two quantitative reproductive traits in 180 genetically identical young adult Caenorhabditis elegans roundworms in a shared environment and performed single-individual transcriptomics on each worm. We identified hundreds of genes for which expression variation was strongly associated with reproductive traits, some of which depended on individuals' historical environments and some of which was random. Multiple small sets of genes together were highly predictive of reproductive traits, explaining on average over half and over a quarter of variation in the two traits. We manipulated mRNA levels of predictive genes to identify a set of causal genes, demonstrating the utility of this approach for both prediction and understanding underlying biology. Finally, we found that the chromatin environment of predictive genes was enriched for H3K27 trimethylation, suggesting that gene expression variation may be driven in part by chromatin structure. Together, this work shows that individual, non-genetic differences in gene expression are both highly predictive and causal in shaping reproductive traits.

Cohesin organizes 3D DNA contacts surrounding active enhancers in C. elegans

In mammals, cohesin and CTCF organize the 3D genome into topologically associated domains (TADs) to regulate communication between cis-regulatory elements. Many organisms, including S. cerevisiae, C. elegans, and A. thaliana contain cohesin but lack CTCF. Here, we used C. elegans to investigate the function of cohesin in 3D genome organization in the absence of CTCF. Using Hi-C data, we observe cohesin-dependent features called "fountains", which are also reported in zebrafish and mice. These are population average reflections of DNA loops originating from distinct genomic regions and are ~20-40 kb in C. elegans. Hi-C analysis upon cohesin and WAPL depletion support the idea that cohesin is preferentially loaded at NIPBL occupied sites and loop extrudes in an effectively two-sided manner. ChIP-seq analyses show that cohesin translocation along the fountain trajectory depends on a fully intact complex and is extended upon WAPL-1 depletion. Hi-C contact patterns at individual fountains suggest that cohesin processivity is unequal on each side, possibly due to collision with cohesin loaded from surrounding sites. The putative cohesin loading sites are closest to active enhancers and fountain strength is associated with transcription. Compared to mammals, average processivity of C. elegans cohesin is ~10-fold shorter and NIPBL binding does not depend on cohesin. We propose that preferential loading and loop extrusion by cohesin is an evolutionarily conserved mechanism that regulates the 3D interactions of enhancers in animal genomes.

Condensin I folds the Caenorhabditis elegans genome

The structural maintenance of chromosome (SMC) complexes-cohesin and condensins-are crucial for chromosome separation and compaction during cell division. During the interphase, mammalian cohesins additionally fold the genome into loops and domains. Here we show that, in Caenorhabditis elegans, a species with holocentric chromosomes, condensin I is the primary, long-range loop extruder. The loss of condensin I and its X-specific variant, condensin IDC, leads to genome-wide decompaction, chromosome mixing and disappearance of X-specific topologically associating domains, while reinforcing fine-scale epigenomic compartments. In addition, condensin I/IDC inactivation led to the upregulation of X-linked genes and unveiled nuclear bodies grouping together binding sites for the X-targeting loading complex of condensin IDC. C. elegans condensin I/IDC thus uniquely organizes holocentric interphase chromosomes, akin to cohesin in mammals, as well as regulates X-chromosome gene expression.

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 eXcision Repair-sequencing (XR-seq) in the model organism Caenorhabditis elegans to generate genome-wide repair maps in a wild-type strain with normal excision repair, a strain lacking TCR (csb-1), and a strain 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.

Clock gene homologs lin-42 and kin-20 regulate circadian rhythms in C. elegans

Circadian rhythms are endogenous oscillations in nearly all organisms, from prokaryotes to humans, allowing them to adapt to cyclical environments for close to 24 h. Circadian rhythms are regulated by a central clock, based on a transcription-translation feedback loop. One important protein in the central loop in metazoan clocks is PERIOD, which is regulated in part by Casein kinase 1ε/δ (CK1ε/δ) phosphorylation. In the nematode Caenorhabditis elegans, period and casein kinase 1ε/δ are conserved as lin-42 and kin-20, respectively. Here, we studied the involvement of lin-42 and kin-20 in the circadian rhythms of the adult nematode using a bioluminescence-based circadian transcriptional reporter. We show that mutations of lin-42 and kin-20 generate a significantly longer endogenous period, suggesting a role for both genes in the nematode circadian clock, as in other organisms. These phenotypes can be partially rescued by overexpression of either gene under their native promoter. Both proteins are expressed in neurons and epidermal seam cells, as well as in other cells. Depletion of LIN-42 and KIN-20, specifically in neuronal cells after development, was sufficient to lengthen the period of oscillating sur-5 expression. Therefore, we conclude that LIN-42 and KIN-20 are critical regulators of the adult nematode circadian clock through neuronal cells.

Regulation of the circadian clock in C. elegans by clock gene homologs kin-20 and lin-42

Circadian rhythms are endogenous oscillations present in nearly all organisms from prokaryotes to humans, allowing them to adapt to cyclical environments close to 24 hours. Circadian rhythms are regulated by a central clock, which is based on a transcription-translation feedback loop. One important protein in the central loop in metazoan clocks is PERIOD, which is regulated in part by Casein kinase 1 ε/δ (CK1 ε/δ ) phosphorylation. In the nematode Caenorhabditis elegans , period and casein kinase 1ε/δ are conserved as lin-42 and kin-20 , respectively. Here we studied the involvement of lin-42 and kin-20 in circadian rhythms of the adult nematode using a bioluminescence-based circadian transcriptional reporter. We show that mutations of lin-42 and kin-20 generate a significantly longer endogenous period, suggesting a role for both genes in the nematode circadian clock, as in other organisms. These phenotypes can be partially rescued by overexpression of either gene under their native promoter. Both proteins are expressed in neurons and seam cells, a population of epidermal stem cells in C. elegans that undergo multiple divisions during development. Depletion of LIN-42 and KIN-20 specifically in neuronal cells after development was sufficient to lengthen the period of oscillating sur-5 expression. Therefore, we conclude that LIN-42 and KIN-20 are critical regulators of the adult nematode circadian clock through neuronal cells.

Cisplatin exposure alters tRNA-derived small RNAs but does not affect epimutations in C. elegans

BACKGROUND

The individual lifestyle and environment of an organism can influence its phenotype and potentially the phenotype of its offspring. The different genetic and non-genetic components of the inheritance system and their mutual interactions are key mechanisms to generate inherited phenotypic changes. Epigenetic changes can be transmitted between generations independently from changes in DNA sequence. In Caenorhabditis elegans, epigenetic differences, i.e. epimutations, mediated by small non-coding RNAs, particularly 22G-RNAs, as well as chromatin have been identified, and their average persistence is three to five generations. In addition, previous research showed that some epimutations had a longer duration and concerned genes that were enriched for multiple components of xenobiotic response pathways. These results raise the possibility that environmental stresses might change the rate at which epimutations occur, with potential significance for adaptation.

RESULTS

In this work, we explore this question by propagating C. elegans lines either in control conditions or in moderate or high doses of cisplatin, which introduces genotoxic stress by damaging DNA. Our results show that cisplatin has a limited effect on global small non-coding RNA epimutations and epimutations in gene expression levels. However, cisplatin exposure leads to increased fluctuations in the levels of small non-coding RNAs derived from tRNA cleavage. We show that changes in tRNA-derived small RNAs may be associated with gene expression changes.

CONCLUSIONS

Our work shows that epimutations are not substantially altered by cisplatin exposure but identifies transient changes in tRNA-derived small RNAs as a potential source of variation induced by genotoxic stress.

Chromatin in 3D distinguishes dMes-4/NSD and Hypb/dSet2 in protecting genes from H3K27me3 silencing

Cell type-specific barcoding of genomes requires the establishment of hundreds of heterochromatin domains where heterochromatin-associated repressive complexes hinder chromatin accessibility thereby silencing genes. At heterochromatin-euchromatin borders, regulation of accessibility not only depends on the delimitation of heterochromatin but may also involve interplays with nearby genes and their transcriptional activity, or alternatively on histone modifiers, chromatin barrier insulators, and more global demarcation of chromosomes into 3D compartmentalized domains and topological-associating domain (TADs). Here, we show that depletion of H3K36 di- or tri-methyl histone methyltransferases dMes-4/NSD or Hypb/dSet2 induces reproducible increasing levels of H3K27me3 at heterochromatin borders including in nearby promoters, thereby repressing hundreds of genes. Furthermore, dMes-4/NSD influences genes demarcated by insulators and TAD borders, within chromatin hubs, unlike transcription-coupled action of Hypb/dSet2 that protects genes independently of TADs. Insulator mutants recapitulate the increase of H3K27me3 upon dMes-4/NSD depletion unlike Hypb/dSet2. Hi-C data demonstrate how dMes-4/NSD blocks propagation of long-range interactions onto active regions. Our data highlight distinct mechanisms protecting genes from H3K27me3 silencing, highlighting a direct influence of H3K36me on repressive TADs.

Tissue-specific chromatin-binding patterns of Caenorhabditis elegans heterochromatin proteins HPL-1 and HPL-2 reveal differential roles in the regulation of gene expression

Heterochromatin is characterized by an enrichment of repetitive elements and low gene density and is often maintained in a repressed state across cell division and differentiation. The silencing is mainly regulated by repressive histone marks such as H3K9 and H3K27 methylated forms and the heterochromatin protein 1 (HP1) family. Here, we analyzed in a tissue-specific manner the binding profile of the two HP1 homologs in Caenorhabditis elegans, HPL-1 and HPL-2, at the L4 developmental stage. We identified the genome-wide binding profile of intestinal and hypodermal HPL-2 and intestinal HPL-1 and compared them with heterochromatin marks and other features. HPL-2 associated preferentially to the distal arms of autosomes and correlated positively with the methylated forms of H3K9 and H3K27. HPL-1 was also enriched in regions containing H3K9me3 and H3K27me3 but exhibited a more even distribution between autosome arms and centers. HPL-2 showed a differential tissue-specific enrichment for repetitive elements conversely with HPL-1, which exhibited a poor association. Finally, we found a significant intersection of genomic regions bound by the BLMP-1/PRDM1 transcription factor and intestinal HPL-1, suggesting a corepressive role during cell differentiation. Our study uncovers both shared and singular properties of conserved HP1 proteins, providing information about genomic binding preferences in relation to their role as heterochromatic markers.

Distinct roles for canonical and variant histone H3 lysine-36 in Polycomb silencing

Polycomb complexes regulate cell type-specific gene expression programs through heritable silencing of target genes. Trimethylation of histone H3 lysine 27 (H3K27me3) is essential for this process. Perturbation of H3K36 is thought to interfere with H3K27me3. We show that mutants of Drosophila replication-dependent (H3.2K36R) or replication-independent (H3.3K36R) histone H3 genes generally maintain Polycomb silencing and reach later stages of development. In contrast, combined (H3.3K36RH3.2K36R) mutants display widespread Hox gene misexpression and fail to develop past the first larval stage. Chromatin profiling revealed that the H3.2K36R mutation disrupts H3K27me3 levels broadly throughout silenced domains, whereas these regions are mostly unaffected in H3.3K36R animals. Analysis of H3.3 distributions showed that this histone is enriched at presumptive Polycomb response elements located outside of silenced domains but relatively depleted from those inside. We conclude that H3.2 and H3.3 K36 residues collaborate to repress Hox genes using different mechanisms.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

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

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

Maternal H3K36 and H3K27 HMTs protect germline development via regulation of the transcription factor LIN-15B

Maternally synthesized products play critical roles in the development of offspring. A premier example is the Caenorhabditis elegans H3K36 methyltransferase MES-4, which is essential for germline survival and development in offspring. How maternal MES-4 protects the germline is not well understood, but its role in H3K36 methylation hinted that it may regulate gene expression in primordial germ cells (PGCs). We tested this hypothesis by profiling transcripts from nascent germlines (PGCs and their descendants) dissected from wild-type and mes-4 mutant (lacking maternal and zygotic MES-4) larvae. mes-4 nascent germlines displayed downregulation of some germline genes, upregulation of some somatic genes, and dramatic upregulation of hundreds of genes on the X chromosome. We demonstrated that upregulation of one or more genes on the X is the cause of germline death by generating and analyzing mes-4 mutants that inherited different endowments of X chromosome(s). Intriguingly, removal of the THAP transcription factor LIN-15B from mes-4 mutants reduced X misexpression and prevented germline death. lin-15B is X-linked and misexpressed in mes-4 PGCs, identifying it as a critical target for MES-4 repression. The above findings extend to the H3K27 methyltransferase MES-2/3/6, the C. elegans version of polycomb repressive complex 2. We propose that maternal MES-4 and PRC2 cooperate to protect germline survival by preventing synthesis of germline-toxic products encoded by genes on the X chromosome, including the key transcription factor LIN-15B.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Accessible Region Conformation Capture (ARC-C) gives high-resolution insights into genome architecture and regulation

Nuclear organization and chromatin interactions are important for genome function, yet determining chromatin connections at high resolution remains a major challenge. To address this, we developed Accessible Region Conformation Capture (ARC-C), which profiles interactions between regulatory elements genome-wide without a capture step. Applied to Caenorhabditis elegans, ARC-C identifies approximately 15,000 significant interactions between regulatory elements at 500-bp resolution. Of 105 TFs or chromatin regulators tested, we find that the binding sites of 60 are enriched for interacting with each other, making them candidates for mediating interactions. These include cohesin and condensin II. Applying ARC-C to a mutant of transcription factor BLMP-1 detected changes in interactions between its targets. ARC-C simultaneously profiles domain-level architecture, and we observe that C. elegans chromatin domains defined by either active or repressive modifications form topologically associating domains (TADs) that interact with A/B (active/inactive) compartment-like structure. Furthermore, we discover that inactive compartment interactions are dependent on H3K9 methylation. ARC-C is a powerful new tool to interrogate genome architecture and regulatory interactions at high resolution.

Mutation rate and spectrum in obligately outcrossing Caenorhabditis elegans mutation accumulation lines subjected to RNAi-induced knockdown of the mismatch repair gene msh-2

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutSα homolog, msh-2, in Caenorhabditis elegans male–female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male–female background substantially increased the frequency of nuclear base substitutions (∼23×) and small indels (∼328×) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

Coordinated maintenance of H3K36/K27 methylation by histone demethylases preserves germ cell identity and immortality

Germ cells have evolved unique mechanisms to ensure the transmission of genetically and nongenetically encoded information, whose alteration compromises germ cell immortality. Chromatin factors play fundamental roles in these mechanisms. H3K36 and H3K27 methyltransferases shape and propagate a pattern of histone methylation essential for C. elegans germ cell maintenance, but the role of respective histone demethylases remains unexplored. Here, we show that jmjd-5 regulates H3K36me2 and H3K27me3 levels, preserves germline immortality, and protects germ cell identity by controlling gene expression. The transcriptional and biological effects of jmjd-5 loss can be hindered by the removal of H3K27demethylases, indicating that H3K36/K27 demethylases act in a transcriptional framework and promote the balance between H3K36 and H3K27 methylation required for germ cell immortality. Furthermore, we find that in wild-type, but not in jmjd-5 mutants, alterations of H3K36 methylation and transcription occur at high temperature, suggesting a role for jmjd-5 in adaptation to environmental changes.

•

jmjd-5 is required for germ cell immortality at high temperature

•

jmjd-5 sustains the expression of germline genes and represses somatic fate

•

Mutations in jmjd-5 result in a global increase of H3K36me2 and H3K27me3

•

Ablation of H3K27 demethylases counteracts the effects of jmjd-5 mutations

Transmission of chromatin states across generations in C. elegans

Epigenetic inheritance refers to the transmission of phenotypes across generations without affecting the genomic DNA sequence. Even though it has been documented in many species in fungi, animals and plants, the mechanisms underlying epigenetic inheritance are not fully uncovered. Epialleles, the heritable units of epigenetic information, can take the form of several biomolecules, including histones and their post-translational modifications (PTMs). Here, we review the recent advances in the understanding of the transmission of histone variants and histone PTM patterns across generations in C. elegans. We provide a general overview of the intergenerational and transgenerational inheritance of histone PTMs and their modifiers and discuss the interplay among different histone PTMs. We also evaluate soma-germ line communication and its impact on the inheritance of epigenetic traits.

The zinc-finger protein OEF-1 stabilizes histone modification patterns and promotes efficient splicing in the C. elegans germ line

To ensure stable transmission of genetic information to the next generation, germ cells frequently silence sex chromosomes, as well as autosomal loci that promote inappropriate differentiation programs. In Caenorhabditis elegans, silenced and active genomic domains are established in germ cells by the histone modification complexes MES-2/3/6 and MES-4, which promote silent and active chromatin states, respectively. These states are generally mutually exclusive and modulation of one state influences the pattern of the other. Here, we identify the zinc-finger protein OEF-1 as a novel modifier of this epigenetic balance in the C. elegans germline. Loss of oef-1 genetically enhances mes mutant phenotypes. Moreover, OEF-1 binding correlates with the active modification H3K36me3 and sustains H3K36me3 levels in the absence of MES-4 activity. OEF-1 also promotes efficient mRNA splicing activity, a process that is influenced by H3K36me3 levels. Finally, OEF-1 limits deposition of the silencing modification H3K27me3 on the X chromosome and at repressed autosomal loci. We propose that OEF-1 might act as an intermediary to mediate the downstream effects of H3K36me3 that promote transcript integrity, and indirectly affect gene silencing as a consequence.

Harnessing model organism genomics to underpin the machine learning-based prediction of essential genes in eukaryotes - Biotechnological implications

The availability of high-quality genomes and advances in functional genomics have enabled large-scale studies of essential genes in model eukaryotes, including the 'elegant worm' (Caenorhabditis elegans; Nematoda) and the 'vinegar fly' (Drosophila melanogaster; Arthropoda). However, this is not the case for other, much less-studied organisms, such as socioeconomically important parasites, for which functional genomic platforms usually do not exist. Thus, there is a need to develop innovative techniques or approaches for the prediction, identification and investigation of essential genes. A key approach that could enable the prediction of such genes is machine learning (ML). Here, we undertake an historical review of experimental and computational approaches employed for the characterisation of essential genes in eukaryotes, with a particular focus on model ecdysozoans (C. elegans and D. melanogaster), and discuss the possible applicability of ML-approaches to organisms such as socioeconomically important parasites. We highlight some recent results showing that high-performance ML, combined with feature engineering, allows a reliable prediction of essential genes from extensive, publicly available 'omic data sets, with major potential to prioritise such genes (with statistical confidence) for subsequent functional genomic validation. These findings could 'open the door' to fundamental and applied research areas. Evidence of some commonality in the essential gene-complement between these two organisms indicates that an ML-engineering approach could find broader applicability to ecdysozoans such as parasitic nematodes or arthropods, provided that suitably large and informative data sets become/are available for proper feature engineering, and for the robust training and validation of algorithms. This area warrants detailed exploration to, for example, facilitate the identification and characterisation of essential molecules as novel targets for drugs and vaccines against parasitic diseases. This focus is particularly important, given the substantial impact that such diseases have worldwide, and the current challenges associated with their prevention and control and with drug resistance in parasite populations.

Single-cell dynamics of chromatin activity during cell lineage differentiation in Caenorhabditis elegans embryos

Elucidating the chromatin dynamics that orchestrate embryogenesis is a fundamental question in developmental biology. Here, we exploit position effects on expression as an indicator of chromatin activity and infer the chromatin activity landscape in every lineaged cell during Caenorhabditis elegans early embryogenesis. Systems-level analyses reveal that chromatin activity distinguishes cellular states and correlates with fate patterning in the early embryos. As cell lineage unfolds, chromatin activity diversifies in a lineage-dependent manner, with switch-like changes accompanying anterior-posterior fate asymmetry and characteristic landscapes being established in different cell lineages. Upon tissue differentiation, cellular chromatin from distinct lineages converges according to tissue types but retains stable memories of lineage history, contributing to intra-tissue cell heterogeneity. However, the chromatin landscapes of cells organized in a left-right symmetric pattern are predetermined to be analogous in early progenitors so as to pre-set equivalent states. Finally, genome-wide analysis identifies many regions exhibiting concordant chromatin activity changes that mediate the co-regulation of functionally related genes during differentiation. Collectively, our study reveals the developmental and genomic dynamics of chromatin activity at the single-cell level.

Long-term experimental evolution reveals purifying selection on piRNA-mediated control of transposable element expression

Background

Transposable elements (TEs) are an almost universal constituent of eukaryotic genomes. In animals, Piwi-interacting small RNAs (piRNAs) and repressive chromatin often play crucial roles in preventing TE transcription and thus restricting TE activity. Nevertheless, TE content varies widely across eukaryotes and the dynamics of TE activity and TE silencing across evolutionary time is poorly understood.

Results

Here, we used experimentally evolved populations of C. elegans to study the dynamics of TE expression over 409 generations. The experimental populations were evolved at population sizes of 1, 10 and 100 individuals to manipulate the efficiency of natural selection versus genetic drift. We demonstrate increased TE expression relative to the ancestral population, with the largest increases occurring in the smallest populations. We show that the transcriptional activation of TEs within active regions of the genome is associated with failure of piRNA-mediated silencing, whilst desilenced TEs in repressed chromatin domains retain small RNAs. Additionally, we find that the sequence context of the surrounding region influences the propensity of TEs to lose silencing through failure of small RNA-mediated silencing.

Conclusions

Our results show that natural selection in C. elegans is responsible for maintaining low levels of TE expression, and provide new insights into the epigenomic features responsible.

Broad genic repression domains signify enhanced silencing of oncogenes

Cancers result from a set of genetic and epigenetic alterations. Most known oncogenes were identified by gain-of-function mutations in cancer, yet little is known about their epigenetic features. Through integrative analysis of 11,596 epigenomic profiles and mutations from >8200 tumor-normal pairs, we discover broad genic repression domains (BGRD) on chromatin as an epigenetic signature for oncogenes. A BGRD is a widespread enrichment domain of the repressive histone modification H3K27me3 and is further enriched with multiple other repressive marks including H3K9me3, H3K9me2, and H3K27me2. Further, BGRD displays widespread enrichment of repressed cis-regulatory elements. Shortening of BGRDs is linked to derepression of transcription. BGRDs at oncogenes tend to be conserved across normal cell types. Putative tumor-promoting genes and lncRNAs defined using BGRDs are experimentally verified as required for cancer phenotypes. Therefore, BGRDs play key roles in epigenetic regulation of cancer and provide a direction for mutation-independent discovery of oncogenes.

Epigenetically altered genes can have a key role in cancer pathobiology but epigenetic signatures that distinguish oncogenes are not yet known. Here, the authors identify broad genic repression domains, defined by widespread H3K27me3 modification, as an epigenetic signature to provide mutation-independent information for discovery of potential oncogenes.

Alanine tRNAs Translate Environment Into Behavior in Caenorhabditis elegans

Caenorhabditis elegans nematodes produce and maintain imprints of attractive chemosensory cues to which they are exposed early in life. Early odor-exposure increases adult chemo-attraction to the same cues. Imprinting is transiently or stably inherited, depending on the number of exposed generations. We show here that the Alanine tRNA (UGC) plays a central role in regulating C. elegans chemo-attraction. Naive worms fed on tRNAAla (UGC) purified from odor-experienced worms, acquire odor-specific imprints. Chemo-attractive responses require the tRNA-modifying Elongator complex sub-units 1 (elpc-1) and 3 (elpc-3) genes. elpc-3 deletions impair chemo-attraction, which is fully restored by wild-type tRNAAla (UGC) feeding. A stably inherited decrease of odor-specific responses ensues from early odor-exposition of elpc-1 deletion mutants. tRNAAla (UGC) may adopt various chemical forms to mediate the cross-talk between innately-programmed and environment-directed chemo-attractive behavior.

Targeting of the Dosage-Compensated Male X-Chromosome during Early Drosophila Development

Dosage compensation, which corrects for the imbalance in X-linked gene expression between XX females and XY males, represents a model for how genes are targeted for coordinated regulation. However, the mechanism by which dosage compensation complexes identify the X chromosome during early development remains unknown because of the difficulty of sexing embryos before zygotic transcription using X- or Y-linked reporter transgenes. We used meiotic drive to sex Drosophila embryos before zygotic transcription and ChIP-seq to measure the dynamics of dosage compensation factor targeting. The Drosophila male-specific lethal dosage compensation complex (MSLc) requires the ubiquitous zinc-finger protein chromatin-linked adaptor for MSL proteins (CLAMP) to identify the X chromosome. We observe a multi-stage process in which MSLc first identifies CLAMP binding sites throughout the genome, followed by concentration at the strongest X-linked MSLc sites. We provide insight into the dynamics of binding site recognition by a large transcription complex during early development.

Rieder et al. establish a meiotic drive system to study Drosophila X chromosome dosage compensation before the maternal-zygotic transition. This study uncovers another step in the process during which the dosage compensation complex identifies binding sites genome-wide before becoming enriched on the X chromosome.

Epimutations driven by small RNAs arise frequently but most have limited duration in Caenorhabditis elegans

Epigenetic regulation involves changes in gene expression independent of DNA sequence variation that are inherited through cell division. In addition to a fundamental role in cell differentiation, some epigenetic changes can also be transmitted transgenerationally through meiosis. Epigenetic alterations (epimutations) could thus contribute to heritable variation within populations and be subject to evolutionary processes such as natural selection and drift. However, the rate at which epimutations arise and their typical persistence are unknown, making it difficult to evaluate their potential for evolutionary adaptation. Here, we perform a genome-wide study of epimutations in a metazoan organism. We use experimental evolution to characterize the rate, spectrum and stability of epimutations driven by small silencing RNAs in the model nematode Caenorhabditis elegans. We show that epimutations arise spontaneously at a rate approximately 25 times greater than DNA sequence changes and typically have short half-lives of two to three generations. Nevertheless, some epimutations last at least ten generations. Epimutations mediated by small RNAs may thus contribute to evolutionary processes over a short timescale but are unlikely to bring about long-term divergence in the absence of selection.

Caenorhabditis elegans Deficient in DOT-1.1 Exhibit Increases in H3K9me2 at Enhancer and Certain RNAi-Regulated Regions

The methylation of histone H3 at lysine 79 is a feature of open chromatin. It is deposited by the conserved histone methyltransferase DOT1. Recently, DOT1 localization and H3K79 methylation (H3K79me) have been correlated with enhancers in C. elegans and mammalian cells. Since earlier research implicated H3K79me in preventing heterochromatin formation both in yeast and leukemic cells, we sought to inquire whether a H3K79me deficiency would lead to higher levels of heterochromatic histone modifications, specifically H3K9me2, at developmental enhancers in C. elegans. Therefore, we used H3K9me2 ChIP-seq to compare its abundance in control and dot-1.1 loss-of-function mutant worms, as well as in rde-4; dot-1.1 and rde-1; dot-1.1 double mutants. The rde-1 and rde-4 genes are components of the RNAi pathway in C. elegans, and RNAi is known to initiate H3K9 methylation in many organisms, including C. elegans. We have previously shown that dot-1.1(-) lethality is rescued by rde-1 and rde-4 loss-of-function. Here we found that H3K9me2 was elevated in enhancer, but not promoter, regions bound by the DOT-1.1/ZFP-1 complex in dot-1.1(-) worms. We also found increased H3K9me2 at genes targeted by the ALG-3/4-dependent small RNAs and repeat regions. Our results suggest that ectopic H3K9me2 in dot-1.1(-) could, in some cases, be induced by small RNAs.

Predicting gene essentiality in Caenorhabditis elegans by feature engineering and machine-learning

Defining genes that are essential for life has major implications for understanding critical biological processes and mechanisms. Although essential genes have been identified and characterised experimentally using functional genomic tools, it is challenging to predict with confidence such genes from molecular and phenomic data sets using computational methods. Using extensive data sets available for the model organism Caenorhabditis elegans, we constructed here a machine-learning (ML)-based workflow for the prediction of essential genes on a genome-wide scale. We identified strong predictors for such genes and showed that trained ML models consistently achieve highly-accurate classifications. Complementary analyses revealed an association between essential genes and chromosomal location. Our findings reveal that essential genes in C. elegans tend to be located in or near the centre of autosomal chromosomes; are positively correlated with low single nucleotide polymorphim (SNP) densities and epigenetic markers in promoter regions; are involved in protein and nucleotide processing; are transcribed in most cells; are enriched in reproductive tissues or are targets for small RNAs bound to the argonaut CSR-1. Based on these results, we hypothesise an interplay between epigenetic markers and small RNA pathways in the germline, with transcription-based memory; this hypothesis warrants testing. From a technical perspective, further work is needed to evaluate whether the present ML-based approach will be applicable to other metazoans (including Drosophila melanogaster) for which comprehensive data sets (i.e. genomic, transcriptomic, proteomic, variomic, epigenetic and phenomic) are available.

Transcription-independent TFIIIC-bound sites cluster near heterochromatin boundaries within lamina-associated domains in C. elegans

BACKGROUND

Chromatin organization is central to precise control of gene expression. In various eukaryotic species, domains of pervasive cis-chromatin interactions demarcate functional domains of the genomes. In nematode Caenorhabditis elegans, however, pervasive chromatin contact domains are limited to the dosage-compensated sex chromosome, leaving the principle of C. elegans chromatin organization unclear. Transcription factor III C (TFIIIC) is a basal transcription factor complex for RNA polymerase III, and is implicated in chromatin organization. TFIIIC binding without RNA polymerase III co-occupancy, referred to as extra-TFIIIC binding, has been implicated in insulating active and inactive chromatin domains in yeasts, flies, and mammalian cells. Whether extra-TFIIIC sites are present and contribute to chromatin organization in C. elegans remains unknown.

RESULTS

We identified 504 TFIIIC-bound sites absent of RNA polymerase III and TATA-binding protein co-occupancy characteristic of extra-TFIIIC sites in C. elegans embryos. Extra-TFIIIC sites constituted half of all identified TFIIIC binding sites in the genome. Extra-TFIIIC sites formed dense clusters in cis. The clusters of extra-TFIIIC sites were highly over-represented within the distal arm domains of the autosomes that presented a high level of heterochromatin-associated histone H3K9 trimethylation (H3K9me3). Furthermore, extra-TFIIIC clusters were embedded in the lamina-associated domains. Despite the heterochromatin environment of extra-TFIIIC sites, the individual clusters of extra-TFIIIC sites were devoid of and resided near the individual H3K9me3-marked regions.

CONCLUSION

Clusters of extra-TFIIIC sites were pervasive in the arm domains of C. elegans autosomes, near the outer boundaries of H3K9me3-marked regions. Given the reported activity of extra-TFIIIC sites in heterochromatin insulation in yeasts, our observation raised the possibility that TFIIIC may also demarcate heterochromatin in C. elegans.

C. elegans synMuv B proteins regulate spatial and temporal chromatin compaction during development

Tissue-specific establishment of repressive chromatin through creation of compact chromatin domains during development is necessary to ensure proper gene expression and cell fate. Caenorhabditis elegans synMuv B proteins are important for the soma/germline fate decision and mutants demonstrate ectopic germline gene expression in somatic tissue, especially at high temperature. We show that C. elegans synMuv B proteins regulate developmental chromatin compaction and that the timing of chromatin compaction is temperature sensitive in both wild type and synMuv B mutants. Chromatin compaction in mutants is delayed into developmental time periods when zygotic gene expression is upregulated and demonstrates an anterior-to-posterior pattern. Loss of this patterned compaction coincides with the developmental time period of ectopic germline gene expression, which leads to a developmental arrest in synMuv B mutants. Finally, accelerated cell division rates at elevated temperature may contribute to a lack of coordination between expression of tissue specific transcription programs and chromatin compaction at high temperature. Thus, chromatin organization during development is regulated both spatially and temporally by synMuv B proteins to establish repressive chromatin in a tissue-specific manner to ensure proper gene expression.

DOT1L complex suppresses transcription from enhancer elements and ectopic RNAi in Caenorhabditis elegans

Methylation of histone H3 on lysine 79 (H3K79) by DOT1L is associated with actively transcribed genes. Earlier, we described that DOT-1.1, the Caenorhabditis elegans homolog of mammalian DOT1L, cooperates with the chromatin-binding protein ZFP-1 (AF10 homolog) to negatively modulate transcription of highly and widely expressed target genes. Also, the reduction of ZFP-1 levels has consistently been associated with lower efficiency of RNA interference (RNAi) triggered by exogenous double-stranded RNA (dsRNA), but the reason for this is not clear. Here, we demonstrate that the DOT1L complex suppresses transcription originating from enhancer elements and antisense transcription, thus potentiating the expression of enhancer-regulated genes. We also show that worms lacking H3K79 methylation do not survive, and this lethality is suppressed by a loss of caspase-3 or Dicer complex components that initiate gene silencing response to exogenous dsRNA. Our results suggest that ectopic elevation of endogenous dsRNA directly or indirectly resulting from global misregulation of transcription in DOT1L complex mutants may engage the Dicer complex and, therefore, limit the efficiency of exogenous RNAi.

Epigenetic, transcriptional and phenotypic responses in two generations of Daphnia magna exposed to the DNA methylation inhibitor 5-azacytidine

The water flea Daphnia magna is a keystone species in freshwater ecosystems and has been widely used as a model organism in environmental ecotoxicology. This aquatic crustacean is sensitive to environmental stressors and displays considerable plasticity in adapting to changing environmental conditions. Part of this plasticity may be due to epigenetic regulation of gene expression, including changes to DNA methylation and histone modifications. Because of the generally hypomethylated genome of this species, we hypothesized that the histone code may have an essential role in the epigenetic control and that histone modifications might be an early marker for stress. This study aims to characterize the epigenetic, transcriptional and phenotypic responses and their causal linkages in directly exposed adult (F0) Daphnia and peritoneal exposed neonates (F1) after a chronic (7-day) exposure to a sublethal concentration (10 mg/l) of 5-azacytidine, a well-studied vertebrate DNA methylation inhibitor. Exposure of the F0 generation significantly reduced the cumulative fecundity, accompanied with differential expression of genes in the one-carbon-cycle metabolic pathway. In the epigenome of the F0 generation, a decrease in global DNA methylation, but no significant changes on H3K4me3 or H3K27me3, were observed. In the F1 offspring generation, changes in gene expression, a significant reduction in global DNA methylation and changes in histone modifications were identified. The results indicate that exposure during adulthood may result in more pronounced effects on early development in the offspring generation, though interpretation of the data should be carefully done since both the exposure regime and developmental period is different in the two generations examined. The obtained results improve our understanding of crustacean epigenetics and the tools developed may promote use of epigenetic markers in hazard assessment of environmental stressors.

Interplay between small RNA pathways shapes chromatin landscapes in C. elegans

The nematode Caenorhabditis elegans contains several types of endogenous small interfering RNAs (endo-siRNAs) produced by RNA-dependent RNA polymerase (RdRP) complexes. Both 'silencing' siRNAs bound by Worm-specific Argonautes (WAGO) and 'activating' siRNAs bound by the CSR-1 Argonaute require the DRH-3 helicase, an RdRP component. Here, we show that, in the drh-3(ne4253) mutant deficient in RdRP-produced secondary endo-siRNAs, the silencing histone mark H3K9me3 is largely depleted, whereas in the csr-1 partially rescued null mutant strain (WM193), this mark is ectopically deposited on CSR-1 target genes. Moreover, we observe ectopic H3K9me3 at enhancer elements and an increased number of small RNAs that match enhancers in both drh-3 and csr-1 mutants. Finally, we detect accumulation of H3K27me3 at highly expressed genes in the drh-3(ne4253) mutant, which correlates with their reduced transcription. Our study shows that when abundant RdRP-produced siRNAs are depleted, there is ectopic elevation of noncoding RNAs linked to sites with increased silencing chromatin marks. Moreover, our results suggest that enhancer small RNAs may guide local H3K9 methylation.

Repression of Germline Genes in Caenorhabditis elegans Somatic Tissues by H3K9 Dimethylation of Their Promoters

Repression of germline-promoting genes in somatic cells is critical for somatic development and function. To study how germline genes are repressed in somatic tissues, we analyzed key histone modifications in three Caenorhabditis elegans synMuv B mutants, lin-15B, lin-35, and lin-37-all of which display ectopic expression of germline genes in the soma. LIN-35 and LIN-37 are members of the conserved DREAM complex. LIN-15B has been proposed to work with the DREAM complex but has not been shown biochemically to be a member of the complex. We found that, in wild-type worms, synMuv B target genes and germline genes are enriched for the repressive histone modification dimethylation of histone H3 on lysine 9 (H3K9me2) at their promoters. Genes with H3K9me2 promoter localization are evenly distributed across the autosomes, not biased toward autosomal arms, as are the broad H3K9me2 domains. Both synMuv B targets and germline genes display a dramatic reduction of H3K9me2 promoter localization in lin-15B mutants, but much weaker reduction in lin-35 and lin-37 mutants. This difference between lin-15B and DREAM complex mutants likely represents a difference in molecular function for these synMuv B proteins. In support of the pivotal role of H3K9me2 in regulation of germline genes by LIN-15B, global loss of H3K9me2 but not H3K9me3 results in phenotypes similar to synMuv B mutants, high-temperature larval arrest, and ectopic expression of germline genes in the soma. We propose that LIN-15B-driven enrichment of H3K9me2 at promoters of germline genes contributes to repression of those genes in somatic tissues.

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

Inverted repeat structures are associated with essential and highly expressed genes on C. elegans autosome distal arms

Complementary sequences in cellular transcripts base-pair to form double-stranded RNA (dsRNA) structures. Because transposon-derived repeats often give rise to self-complementary sequences, dsRNA structures are prevalent in eukaryotic genomes, typically occurring in gene introns and untranslated regions (UTRs). However, the regulatory impact of double-stranded structures within genes is not fully understood. We used three independent methods to define loci in Caenorhabditis elegans predicted to form dsRNA and correlated these structures with patterns of gene expression, gene essentiality, and genome organization. As previously observed, dsRNA loci are enriched on distal arms of C. elegans autosomes, where genes typically show less conservation and lower overall expression. In contrast, we find that dsRNAs are associated with essential genes on autosome arms, and dsRNA-associated genes exhibit higher-than-expected expression and histone modification patterns associated with transcriptional elongation. Genes with significant repetitive sequence content are also highly expressed, and, thus, observed gene expression trends may relate either to dsRNA structures or to repeat content. Our results raise the possibility that as-yet-undescribed mechanisms promote expression of loci that produce dsRNAs, despite their well-characterized roles in gene silencing.

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.

Programmed DNA Breaks Activate the Germline Genome in Caenorhabditis elegans

In Caenorhabditis elegans, the primordial germ cells Z2 and Z3 are born during early embryogenesis and then held in a transcriptionally quiescent state where the genome is highly compacted. When hatched L1s feed, the germline genome decompacts, and RNAPII is abruptly and globally activated. A previously documented yet unexplained feature of germline genome activation in the worm is the appearance of numerous DNA breaks coincident with RNAPII transcription. Here, we show that the DNA breaks are induced by topoisomerase II and that they function to recruit the RUVB complex to chromosomes so that RUVB can decompact the chromatin. DNA break- and RUVB-mediated decompaction is required for zygotic genome activation. This work highlights the importance of global chromatin decompaction in the rapid induction of gene expression and shows that one way cells achieve global decompaction is through programmed DNA breaks.

The Germline-Specific Factor OEF-1 Facilitates Coordinated Progression Through Germ Cell Development in Caenorhabditis elegans

The purpose of germ cells is to ensure the faithful transmission of genetic material to the next generation. To develop into mature gametes, germ cells must pass through cell cycle checkpoints while maintaining totipotency and genomic integrity. How germ cells coordinate developmental events while simultaneously protecting their unique fate is not well understood. Here, we characterize a novel nuclear protein, Oocyte-Excluded Factor-1 (OEF-1), with highly specific germline expression in Caenorhabditis elegans OEF-1 is initially detected early in embryogenesis and is expressed in the nuclei of all germ cells during larval stages. In adults, OEF-1 expression abruptly decreases just prior to oocyte differentiation. In oef-1 mutants, the developmental progression of germ cells is accelerated, resulting in subtle defects at multiple stages of germ cell development. Lastly, OEF-1 is primarily associated with the bodies of germline-expressed genes, and as such is excluded from the X chromosome. We hypothesize that OEF-1 may regulate the rate of progression through germ cell development, providing insight into how these critical maturation events are coordinated.

Repressive Chromatin in Caenorhabditis elegans: Establishment, Composition, and Function

Chromatin is organized and compacted in the nucleus through the association of histones and other proteins, which together control genomic activity. Two broad types of chromatin can be distinguished: euchromatin, which is generally transcriptionally active, and heterochromatin, which is repressed. Here we examine the current state of our understanding of repressed chromatin in Caenorhabditis elegans, focusing on roles of histone modifications associated with repression, such as methylation of histone H3 lysine 9 (H3K9me2/3) or the Polycomb Repressive Complex 2 (MES-2/3/6)-deposited modification H3K27me3, and on proteins that recognize these modifications. Proteins involved in chromatin repression are important for development, and have demonstrated roles in nuclear organization, repetitive element silencing, genome integrity, and the regulation of euchromatin. Additionally, chromatin factors participate in repression with small RNA pathways. Recent findings shed light on heterochromatin function and regulation in C. elegans, and should inform our understanding of repressed chromatin in other animals.

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.

Broad Chromatin Domains: An Important Facet of Genome Regulation

Chromatin composition differs across the genome, with distinct compositions characterizing regions associated with different properties and functions. Whereas many histone modifications show local enrichment over genes or regulatory elements, marking can also span large genomic intervals defining broad chromatin domains. Here we highlight structural and functional features of chromatin domains marked by histone modifications, with a particular emphasis on the potential roles of H3K27 methylation domains in the organization and regulation of genome activity in metazoans.

Histone gene replacement reveals a post-transcriptional role for H3K36 in maintaining metazoan transcriptome fidelity

Histone H3 lysine 36 methylation (H3K36me) is thought to participate in a host of co-transcriptional regulatory events. To study the function of this residue independent from the enzymes that modify it, we used a ‘histone replacement’ system in Drosophila to generate a non-modifiable H3K36 lysine-to-arginine (H3K36R) mutant. We observed global dysregulation of mRNA levels in H3K36R animals that correlates with the incidence of H3K36me3. Similar to previous studies, we found that mutation of H3K36 also resulted in H4 hyperacetylation. However, neither cryptic transcription initiation, nor alternative pre-mRNA splicing, contributed to the observed changes in expression, in contrast with previously reported roles for H3K36me. Interestingly, knockdown of the RNA surveillance nuclease, Xrn1, and members of the CCR4-Not deadenylase complex, restored mRNA levels for a class of downregulated, H3K36me3-rich genes. We propose a post-transcriptional role for modification of replication-dependent H3K36 in the control of metazoan gene expression.

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

In a single human cell there is enough DNA to stretch over a meter if laid end to end. To fit this DNA inside the cell – which is less than 20 micrometers in diameter – the DNA is tightly wrapped around millions of proteins known as histones, which look like “beads” along a “string” of DNA. These histones can prevent other proteins from binding to DNA and activating specific genes. Therefore, cells use enzymes to chemically modify histones to allow particular stretches of DNA to be unwrapped at specific times.

Proteins are made up of building blocks called amino acids. A specific amino acid on histones known as H3K36 is modified in certain sections of DNA that suggest it affects the activities of many genes. However, the precise role of this amino acid remains unclear. Previous studies have tried to investigate this by removing the enzymes that modify it, but these enzymes can also modify many other proteins, making it difficult to know what exactly causes the changes in gene activity.

Fruit flies are often used in experiments as models of how genetic processes work in humans and other animals. Like us, fruit flies also package their DNA using histones. To investigate the role of H3K36, Meers et al. generated a mutant fruit fly that has a version of the amino acid that cannot be chemically modified by the normal enzymes. Unexpectedly, the experiments suggest that some changes in gene activity that have been previously reported to be caused by modifying H3K36 might actually be due to other factors. Meers et al. found that H3K36 modifications may instead “mark” certain genes to be more active than they otherwise would be.

These findings provide a starting point for understanding exactly how H3K36 regulates gene activity. The next challenge is to refine our understanding of how H3K36 modification affects genes in cancer and other diseases, which may aid the development of new therapies to treat these conditions.

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