High nucleosome occupancy is encoded at X-linked gene promoters in C. elegans

The role of transcript regions and amino acid choice in nucleosome positioning

Eukaryotic DNA is organized and compacted in a string of nucleosomes, DNA-wrapped protein cylinders. The positions of nucleosomes along DNA are not random but show well-known base pair sequence preferences that result from the sequence-dependent elastic and geometric properties of the DNA double helix. Here, we focus on DNA around transcription start sites, which are known to typically attract nucleosomes in multicellular life forms through their high GC content. We aim to understand how these GC signals, as observed in genome-wide averages, are produced and encoded through different genomic regions (mainly 5' UTRs, coding exons, and introns). Our study uses a bioinformatics approach to decompose the genome-wide GC signal into between-region and within-region signals. We find large differences in GC signal contributions between vertebrates and plants and, remarkably, even between closely related species. Introns contribute most to the GC signal in vertebrates, while in plants the exons dominate. Further, we find signal strengths stronger on DNA than on mRNA, suggesting a biological function of GC signals along the DNA itself, as is the case for nucleosome positioning. Finally, we make the surprising discovery that both the choice of synonymous codons and amino acids contribute to the nucleosome positioning signal.

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

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

Shortest paths through synonymous genomes

The elasticity of the DNA double helix varies with the underlying base pair sequence. This allows one to put mechanical cues into sequences that in turn influence the packaging of DNA into nucleosomes, DNA-wrapped protein cylinders. Nucleosomes dictate a broad range of biological processes, ranging from gene regulation, recombination, and replication to chromosome condensation. Here we map base pair sequences onto graphs and use shortest paths algorithms to determine which DNA stretches are easiest or hardest to bend inside a nucleosome. We further demonstrate how genetic and mechanical information can be multiplexed by studying paths through graphs of synonymous codons. Using this method we find that nucleosomes can be placed by mechanical cues nearly everywhere on the genome of baker's yeast (Saccharomyces cerevisiae).

Cooperation between a hierarchical set of recruitment sites targets the X chromosome for dosage compensation

In many organisms, it remains unclear how X chromosomes are specified for dosage compensation, since DNA sequence motifs shown to be important for dosage compensation complex (DCC) recruitment are themselves not X-specific. Here, we addressed this problem in C. elegans. We found that the DCC recruiter, SDC-2, is required to maintain open chromatin at a small number of primary DCC recruitment sites, whose sequence and genomic context are X-specific. Along the X, primary recruitment sites are interspersed with secondary sites, whose function is X-dependent. A secondary site can ectopically recruit the DCC when additional recruitment sites are inserted either in tandem or at a distance (>30 kb). Deletion of a recruitment site on the X results in reduced DCC binding across several megabases surrounded by topologically associating domain (TAD) boundaries. Our work elucidates that hierarchy and long-distance cooperativity between gene-regulatory elements target a single chromosome for regulation.

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

The DNA inside living cells is organized in structures called chromosomes. In many animals, females have two X chromosomes, whereas males have only one. To ensure that females do not end up with a double dose of the proteins encoded by the genes on the X chromosome, animals use a process called dosage compensation to correct this imbalance. The mechanisms underlying this process vary between species, but they typically involve a regulatory complex that binds to the X chromosomes of one sex to modify gene expression.

Caenorhabditis elegans, for example, is a species of nematode worm in which individuals with two X chromosomes are hermaphrodites and those with one X chromosome are males. In C. elegans, a regulatory complex, called the dosage compensation complex, attaches to both X chromosomes of a hermaphrodite, and reduces the expression of the genes on each by half to match the level seen in the males. Previous research has shown that short DNA sequences, known as motifs, recruit the dosage compensation complex to the X chromosomes. However, these sequences are also found on the other chromosomes and, until now, it was not known why the complex was only recruited to the X chromosomes.

Albritton et al. now show the X chromosomes have a ‘hierarchical’ recruitment system. A few sites on the X chromosomes contain clusters of a specific DNA motif, which initiate the process and attract the dosage compensation complex more strongly than other sites. These ‘strong’ recruitment sites are placed across the length of the X chromosomes and cooperate with several ‘weaker’ ones located in between. This way, multiple recruitment sites can cooperate over a long distance, while non-sex chromosomes, which have only one or two stronger recruitment sites, do not have thisadvantage.

Hierarchy and cooperativity may be general features of gene expression, in which proteins are targeted to chromosomes without the need for having specific motifs at every recruitment site. The way DNA sequences are distributed across the genome may give us clues about their role. Thus, knowing how genomes are structured will help us identify disrupted areas in diseases such as cancer.

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

Genomes of Multicellular Organisms Have Evolved to Attract Nucleosomes to Promoter Regions

Sequences that influence nucleosome positioning in promoter regions, and their relation to gene regulation, have been the topic of much research over the last decade. In yeast, significant nucleosome-depleted regions are found, which facilitate transcription. With the arrival of nucleosome positioning maps for the human genome, it was discovered that in our genome, unlike in that of yeast, promoters encode for high nucleosome occupancy. In this work, we look at the genomes of a range of different organisms, to provide a catalog of nucleosome positioning signals in promoters across the tree of life. We utilize a computational model of the nucleosome, based on crystallographic analyses of the structure and elasticity of the nucleosome, to predict the nucleosome positioning signals in promoter regions. To be able to apply our model to large genomic datasets, we introduce an approximative scheme that makes use of the limited range of correlations in nucleosomal sequence preferences to create a computationally efficient approximation of the full biophysical model. Our predictions show that a clear distinction between unicellular and multicellular life is visible in the intrinsically encoded nucleosome affinity. Furthermore, the strength of the nucleosome positioning signals correlates with the complexity of the organism. We conclude that encoding for high nucleosome occupancy, as in the human genome, is in fact a universal feature of multicellular life.

Nucleosome fragility is associated with future transcriptional response to developmental cues and stress in C. elegans

Nucleosomes have structural and regulatory functions in all eukaryotic DNA-templated processes. The position of nucleosomes on DNA and the stability of the underlying histone–DNA interactions affect the access of regulatory proteins to DNA. Both stability and position are regulated through DNA sequence, histone post-translational modifications, histone variants, chromatin remodelers, and transcription factors. Here, we explored the functional implications of nucleosome properties on gene expression and development in Caenorhabditis elegans embryos. We performed a time-course of micrococcal nuclease (MNase) digestion and measured the relative sensitivity or resistance of nucleosomes throughout the genome. Fragile nucleosomes were defined by nucleosomal DNA fragments that were recovered preferentially in early MNase-digestion time points. Nucleosome fragility was strongly and positively correlated with the AT content of the underlying DNA sequence. There was no correlation between promoter nucleosome fragility and the levels of histone modifications or histone variants. Genes with fragile nucleosomes in their promoters tended to be lowly expressed and expressed in a context-specific way, operating in neuronal response, the immune system, and stress response. In addition to DNA-encoded nucleosome fragility, we also found fragile nucleosomes at locations where we expected to find destabilized nucleosomes, for example, at transcription factor binding sites where nucleosomes compete with DNA-binding factors. Our data suggest that in C. elegans promoters, nucleosome fragility is in large part DNA-encoded and that it poises genes for future context-specific activation in response to environmental stress and developmental cues.

An H4K16 histone acetyltransferase mediates decondensation of the X chromosome in C. elegans males

Background

In C. elegans, in order to equalize gene expression between the sexes and balance X and autosomal expression, two steps are believed to be required. First, an unknown mechanism is hypothesized to upregulate the X chromosome in both sexes. This mechanism balances the X to autosomal expression in males, but creates X overexpression in hermaphrodites. Therefore, to restore the balance, hermaphrodites downregulate gene expression twofold on both X chromosomes. While many studies have focused on X chromosome downregulation, the mechanism of X upregulation is not known.

Results

To gain more insight into X upregulation, we studied the effects of chromatin condensation and histone acetylation on gene expression levels in male C. elegans. We have found that the H4K16 histone acetyltransferase MYS-1/Tip60 mediates dramatic decondensation of the male X chromosome as measured by FISH. However, RNA-seq analysis revealed that MYS-1 contributes only slightly to upregulation of gene expression on the X chromosome. These results suggest that the level of chromosome decondensation does not necessarily correlate with the degree of gene expression change in vivo. Furthermore, the X chromosome is more sensitive to MYS-1-mediated decondensation than the autosomes, despite similar levels of H4K16ac on all chromosomes, as measured by ChIP-seq. H4K16ac levels weakly correlate with gene expression levels on both the X and the autosomes, but highly expressed genes on the X chromosome do not contain exceptionally high levels of H4K16ac.

Conclusion

These results indicate that H4K16ac and chromosome decondensation influence regulation of the male X chromosome; however, they do not fully account for the high levels of gene expression observed on the X chromosomes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0097-x) contains supplementary material, which is available to authorized users.

The time-resolved transcriptome of C. elegans

We generated detailed RNA-seq data for the nematode Caenorhabditis elegans with high temporal resolution in the embryo as well as representative samples from post-embryonic stages across the life cycle. The data reveal that early and late embryogenesis is accompanied by large numbers of genes changing expression, whereas fewer genes are changing in mid-embryogenesis. This lull in genes changing expression correlates with a period during which histone mRNAs produce almost 40% of the RNA-seq reads. We find evidence for many more splice junctions than are annotated in WormBase, with many of these suggesting alternative splice forms, often with differential usage over the life cycle. We annotated internal promoter usage in operons using SL1 and SL2 data. We also uncovered correlated transcriptional programs that span >80 kb. These data provide detailed annotation of the C. elegans transcriptome.

Fine-Scale Crossover Rate Variation on the Caenorhabditis elegans X Chromosome

Meiotic recombination creates genotypic diversity within species. Recombination rates vary substantially across taxa, and the distribution of crossovers can differ significantly among populations and between sexes. Crossover locations within species have been found to vary by chromosome and by position within chromosomes, where most crossover events occur in small regions known as recombination hotspots. However, several species appear to lack hotspots despite significant crossover heterogeneity. The nematode Caenorhabditis elegans was previously found to have the least fine-scale variation in crossover distribution among organisms studied to date. It is unclear whether this pattern extends to the X chromosome given its unique compaction through the pachytene stage of meiotic prophase in hermaphrodites. We generated 798 recombinant nested near-isogenic lines (NILs) with crossovers in a 1.41 Mb region on the left arm of the X chromosome to determine if its recombination landscape is similar to that of the autosomes. We find that the fine-scale variation in crossover rate is lower than that of other model species, and is inconsistent with hotspots. The relationship of genomic features to crossover rate is dependent on scale, with GC content, histone modifications, and nucleosome occupancy being negatively associated with crossovers. We also find that the abundances of 4- to 6-bp DNA motifs significantly explain crossover density. These results are consistent with recombination occurring at unevenly distributed sites of open chromatin.

Nonconsensus Protein Binding to Repetitive DNA Sequence Elements Significantly Affects Eukaryotic Genomes

Recent genome-wide experiments in different eukaryotic genomes provide an unprecedented view of transcription factor (TF) binding locations and of nucleosome occupancy. These experiments revealed that a large fraction of TF binding events occur in regions where only a small number of specific TF binding sites (TFBSs) have been detected. Furthermore, in vitro protein-DNA binding measurements performed for hundreds of TFs indicate that TFs are bound with wide range of affinities to different DNA sequences that lack known consensus motifs. These observations have thus challenged the classical picture of specific protein-DNA binding and strongly suggest the existence of additional recognition mechanisms that affect protein-DNA binding preferences. We have previously demonstrated that repetitive DNA sequence elements characterized by certain symmetries statistically affect protein-DNA binding preferences. We call this binding mechanism nonconsensus protein-DNA binding in order to emphasize the point that specific consensus TFBSs do not contribute to this effect. In this paper, using the simple statistical mechanics model developed previously, we calculate the nonconsensus protein-DNA binding free energy for the entire C. elegans and D. melanogaster genomes. Using the available chromatin immunoprecipitation followed by sequencing (ChIP-seq) results on TF-DNA binding preferences for ~100 TFs, we show that DNA sequences characterized by low predicted free energy of nonconsensus binding have statistically higher experimental TF occupancy and lower nucleosome occupancy than sequences characterized by high free energy of nonconsensus binding. This is in agreement with our previous analysis performed for the yeast genome. We suggest therefore that nonconsensus protein-DNA binding assists the formation of nucleosome-free regions, as TFs outcompete nucleosomes at genomic locations with enhanced nonconsensus binding. In addition, here we perform a new, large-scale analysis using in vitro TF-DNA preferences obtained from the universal protein binding microarrays (PBM) for ~90 eukaryotic TFs belonging to 22 different DNA-binding domain types. As a result of this new analysis, we conclude that nonconsensus protein-DNA binding is a widespread phenomenon that significantly affects protein-DNA binding preferences and need not require the presence of consensus (specific) TFBSs in order to achieve genome-wide TF-DNA binding specificity.

Interactions between proteins and DNA trigger many important biological processes. Therefore, to fully understand how the information encoded on the DNA transcribes into RNA, which in turn translates into proteins in the cell, we need to unravel the molecular design principles of protein-DNA interactions. It is known that many interactions occur when a protein is attracted to a specific short segment on the DNA called a specific protein-DNA binding motif. Strikingly, recent experiments revealed that many regulatory proteins reproducibly bind to different regions on the DNA lacking such specific motifs. This suggests that fundamental molecular mechanisms responsible for protein-DNA recognition specificity are not fully understood. Here, using high-throughput protein-DNA binding data obtained by two entirely different methods for ~100 TFs in each case, we show that DNA regions possessing certain repetitive sequence elements exert the statistical attractive potential on DNA-binding proteins, and as a result, such DNA regions are enriched in bound proteins. This is in agreement with our previous analysis performed for the yeast genome. We use the term nonconsensus protein-DNA binding in order to describe protein-DNA interactions that occur in the absence of specific protein-DNA binding motifs. Here we demonstrate that the identified nonconsensus effect is highly significant for a variety of organismal genomes and it affects protein-DNA binding preferences and nucleosome occupancy at the genome-wide level.

Linking dosage compensation and X chromosome nuclear organization in C. elegans

Animal sex is determined by the number of X chromosomes in many species, creating unequal gene dosage (aneuploidy) between sexes. Dosage Compensation mechanisms equalize this dosage difference by regulating X-linked gene expression. In the nematode C. elegans the current model suggests that DC is achieved by a 2-fold transcriptional downregulation in hermaphrodites mediated by the Dosage Compensation Complex (DCC), which restricts access to RNA Polymerase II by an unknown mechanism. Taking a nuclear organization point of view, we showed that the male X chromosome resides in the pore proximal subnuclear compartment whereas the DCC bound to the X, inhibits this spatial organization in the hermaphrodites. Here we discuss our results and propose a model that reassigns the role of DCC from repression of genes to inhibition of activation.

Differential Nucleosome Occupancies across Oct4-Sox2 Binding Sites in Murine Embryonic Stem Cells

The binding sequence for any transcription factor can be found millions of times within a genome, yet only a small fraction of these sequences encode functional transcription factor binding sites. One of the reasons for this dichotomy is that many other factors, such as nucleosomes, compete for binding. To study how the competition between nucleosomes and transcription factors helps determine a functional transcription factor site from a predicted transcription factor site, we compared experimentally-generated in vitro nucleosome occupancy with in vivo nucleosome occupancy and transcription factor binding in murine embryonic stem cells. Using a solution hybridization enrichment technique, we generated a high-resolution nucleosome map from targeted regions of the genome containing predicted sites and functional sites of Oct4/Sox2 regulation. We found that at Pax6 and Nes, which are bivalently poised in stem cells, functional Oct4 and Sox2 sites show high amounts of in vivo nucleosome displacement compared to in vitro. Oct4 and Sox2, which are active, show no significant displacement of in vivo nucleosomes at functional sites, similar to nonfunctional Oct4/Sox2 binding. This study highlights a complex interplay between Oct4 and Sox2 transcription factors and nucleosomes among different target genes, which may result in distinct patterns of stem cell gene regulation.

Mechanisms of x chromosome dosage compensation

In many animals, males have one X and females have two X chromosomes. The difference in X chromosome dosage between the two sexes is compensated by mechanisms that regulate X chromosome transcription. Recent advances in genomic techniques have provided new insights into the molecular mechanisms of X chromosome dosage compensation. In this review, I summarize our current understanding of dosage imbalance in general, and then review the molecular mechanisms of X chromosome dosage compensation with an emphasis on the parallels and differences between the three well-studied model systems, M. musculus, D. melanogaster and C. elegans.

Genome-wide analysis of condensin binding in Caenorhabditis elegans

Background

Condensins are multi-subunit protein complexes that are essential for chromosome condensation during mitosis and meiosis, and play key roles in transcription regulation during interphase. Metazoans contain two condensins, I and II, which perform different functions and localize to different chromosomal regions. Caenorhabditis elegans contains a third condensin, I^DC, that is targeted to and represses transcription of the X chromosome for dosage compensation.

Results

To understand condensin binding and function, we performed ChIP-seq analysis of C. elegans condensins in mixed developmental stage embryos, which contain predominantly interphase nuclei. Condensins bind to a subset of active promoters, tRNA genes and putative enhancers. Expression analysis in kle-2-mutant larvae suggests that the primary effect of condensin II on transcription is repression. A DNA sequence motif, GCGC, is enriched at condensin II binding sites. A sequence extension of this core motif, AGGG, creates the condensin I^DC motif. In addition to differences in recruitment that result in X-enrichment of condensin I^DC and condensin II binding to all chromosomes, we provide evidence for a shared recruitment mechanism, as condensin I^DC recruiter SDC-2 also recruits condensin II to the condensin I^DC recruitment sites on the X. In addition, we found that condensin sites overlap extensively with the cohesin loader SCC-2, and that SDC-2 also recruits SCC-2 to the condensin I^DC recruitment sites.

Conclusions

Our results provide the first genome-wide view of metazoan condensin II binding in interphase, define putative recruitment motifs, and illustrate shared loading mechanisms for condensin I^DC and condensin II.

pENCODE: a plant encyclopedia of DNA elements

ENCODE projects exist for many eukaryotes, including humans, but as of yet no defined project exists for plants. A plant ENCODE would be invaluable to the research community and could be more readily produced than its metazoan equivalents by capitalizing on the preexisting infrastructure provided from similar projects. Collecting and normalizing plant epigenomic data for a range of species will facilitate hypothesis generation, cross-species comparisons, annotation of genomes, and an understanding of epigenomic functions throughout plant evolution. Here, we discuss the need for such a project, outline the challenges it faces, and suggest ways forward to build a plant ENCODE.

BRG1/SMARCA4 inactivation promotes non-small cell lung cancer aggressiveness by altering chromatin organization

SWI/SNF chromatin remodeling complexes regulate critical cellular processes, including cell-cycle control, programmed cell death, differentiation, genomic instability, and DNA repair. Inactivation of this class of chromatin remodeling complex has been associated with a variety of malignancies, including lung, ovarian, renal, liver, and pediatric cancers. In particular, approximately 10% of primary human lung non-small cell lung cancers (NSCLC) display attenuations in the BRG1 ATPase, a core factor in SWI/SNF complexes. To evaluate the role of BRG1 attenuation in NSCLC development, we examined the effect of BRG1 silencing in primary and established human NSCLC cells. BRG1 loss altered cellular morphology and increased tumorigenic potential. Gene expression analyses showed reduced expression of genes known to be associated with progression of human NSCLC. We demonstrated that BRG1 losses in NSCLC cells were associated with variations in chromatin structure, including differences in nucleosome positioning and occupancy surrounding transcriptional start sites of disease-relevant genes. Our results offer direct evidence that BRG1 attenuation contributes to NSCLC aggressiveness by altering nucleosome positioning at a wide range of genes, including key cancer-associated genes.

Crossover recombination mediated by HIM-18/SLX4-associated nucleases

Meiosis is a specialized cell division program that results in the formation of haploid gametes (i.e., sperm and eggs) from diploid parental cells, and is essential for all sexually reproducing organisms. Crossover formation, the reciprocal exchange of genetic information during recombination, is critical for accurate meiotic chromosome segregation. Misregulation of crossover formation leads to genomic instability and aneuploidy (cells with the incorrect number of chromosomes), resulting in tumorigenesis, birth defects, miscarriages, and infertility in humans. Recently, a shuriken/Swiss army knife-like multi-nuclease complex has been implicated in processing various types of DNA repair intermediates. However, how these nucleases coordinate their functions during repair remained unclear. Our studies in C. elegans revealed genetic redundancies between these nucleases for meiotic crossover formation and that they promote distinct crossover control at different chromosome regions. Specifically, XPF-1 acts redundantly with both MUS-81 and SLX-1 to resolve Holliday junction recombination intermediates into crossover products at designated future crossover sites on chromosome arms. In contrast, SLX-1 is required for suppression of crossovers at the center region of chromosomes. Altogether, our studies have shed light on the interplay between structure-specific endonucleases and uncovered their ability to exert either positive or negative meiotic crossover control on a chromosome region-specific basis.

Regulation of the X chromosomes in Caenorhabditis elegans

Dosage compensation, which regulates the expression of genes residing on the sex chromosomes, has provided valuable insights into chromatin-based mechanisms of gene regulation. The nematode Caenorhabditis elegans has adopted various strategies to down-regulate and even nearly silence the X chromosomes. This article discusses the different chromatin-based strategies used in somatic tissues and in the germline to modulate gene expression from the C. elegans X chromosomes and compares these strategies to those used by other organisms to cope with similar X-chromosome dosage differences.

Transcriptional regulation of gene expression in C. elegans

Protein coding gene sequences are converted to mRNA by the highly regulated process of transcription. The precise temporal and spatial control of transcription for many genes is an essential part of development in metazoans. Thus, understanding the molecular mechanisms underlying transcriptional control is essential to understanding cell fate determination during embryogenesis, post-embryonic development, many environmental interactions, and disease-related processes. Studies of transcriptional regulation in C. elegans exploit its genomic simplicity and physical characteristics to define regulatory events with single-cell and minute-time-scale resolution. When combined with the genetics of the system, C. elegans offers a unique and powerful vantage point from which to study how chromatin-associated proteins and their modifications interact with transcription factors and their binding sites to yield precise control of gene expression through transcriptional regulation.

Dynamic nucleosome organization at hox promoters during zebrafish embryogenesis

Nucleosome organization at promoter regions plays an important role in regulating gene activity. Genome-wide studies in yeast, flies, worms, mammalian embryonic stem cells and transformed cell lines have found well-positioned nucleosomes flanking a nucleosome depleted region (NDR) at transcription start sites. This nucleosome arrangement depends on DNA sequence (cis-elements) as well as DNA binding factors and ATP-dependent chromatin modifiers (trans-factors). However, little is understood about how the nascent embryonic genome positions nucleosomes during development. This is particularly intriguing since the embryonic genome must undergo a broad reprogramming event upon fusion of sperm and oocyte. Using four stages of early embryonic zebrafish development, we map nucleosome positions at the promoter region of 37 zebrafish hox genes. We find that nucleosome arrangement at the hox promoters is a progressive process that takes place over several stages. At stages immediately after fertilization, nucleosomes appear to be largely disordered at hox promoter regions. At stages after activation of the embryonic genome, nucleosomes are detectable at hox promoters, with positions becoming more uniform and more highly occupied. Since the genomic sequence is invariant during embryogenesis, this progressive change in nucleosome arrangement suggests that trans-factors play an important role in organizing nucleosomes during embryogenesis. Separating hox genes into expressed and non-expressed groups shows that expressed promoters have better positioned and occupied nucleosomes, as well as distinct NDRs, than non-expressed promoters. Finally, by blocking the retinoic acid-signaling pathway, we disrupt early hox gene transcription, but observe no effect on nucleosome positions, suggesting that active hox transcription is not a driving force behind the arrangement of nucleosomes at the promoters of hox genes during early development.

Global remodeling of nucleosome positions in C. elegans

Background

Eukaryotic chromatin architecture is affected by intrinsic histone-DNA sequence preferences, steric exclusion between nucleosome particles, formation of higher-order structures, and in vivo activity of chromatin remodeling enzymes.

Results

To disentangle sequence-dependent nucleosome positioning from the other factors, we have created two high-throughput maps of nucleosomes assembled in vitro on genomic DNA from the nematode worm Caenorhabditis elegans. A comparison of in vitro nucleosome positions with those observed in a mixed-stage, mixed-tissue population of C. elegans cells reveals that in vivo sequence preferences are modified on the genomic scale. Indeed, G/C dinucleotides are predicted to be most favorable for nucleosome formation in vitro but not in vivo. Nucleosome sequence read coverage in vivo is distinctly lower in chromosome arms than in central regions; the observed changes in apparent nucleosome sequence specificity, likely due to genome-wide chromatin remodeler activity, contribute to the formation of these megabase-scale chromatin domains. We also observe that the majority of well-positioned in vivo nucleosomes do not occupy thermodynamically favorable sequences observed in vitro. Finally, we find that exons are intrinsically more amenable to nucleosome formation compared to introns. Nucleosome occupancy of introns and exons consistently increases with G/C content in vitro but not in vivo, in agreement with our observation that G/C dinucleotide enrichment does not strongly promote in vivo nucleosome formation.

Conclusions

Our findings highlight the importance of both sequence specificity and active nucleosome repositioning in creating large-scale chromatin domains, and the antagonistic roles of intrinsic sequence preferences and chromatin remodelers in C. elegans.

Sequence read data has been deposited into Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra; accession number SRA050182). Additional data, software and computational predictions are available on the Nucleosome Explorer website (http://nucleosome.rutgers.edu).

DNA sequence preferences of transcriptional activators correlate more strongly than repressors with nucleosomes

Transcription factors (TFs) and histone octamers are two abundant classes of DNA binding proteins that coordinate the transcriptional program in cells. Detailed studies of individual TFs have shown that TFs bind to nucleosome-occluded DNA sequences and induce nucleosome disruption/repositioning, while recent global studies suggest this is not the only mechanism used by all TFs. We have analyzed to what extent the intrinsic DNA binding preferences of TFs and histones play a role in determining nucleosome occupancy, in addition to nonintrinsic factors such as the enzymatic activity of chromatin remodelers. The majority of TFs in budding yeast have an intrinsic sequence preference overlapping with nucleosomal histones. TFs with intrinsic DNA binding properties highly correlated with those of histones tend to be associated with gene activation and might compete with histones to bind to genomic DNA. Consistent with this, we show that activators induce more nucleosome disruption upon transcriptional activation than repressors.

H4K20me1 contributes to downregulation of X-linked genes for C. elegans dosage compensation

The Caenorhabditis elegans dosage compensation complex (DCC) equalizes X-chromosome gene dosage between XO males and XX hermaphrodites by two-fold repression of X-linked gene expression in hermaphrodites. The DCC localizes to the X chromosomes in hermaphrodites but not in males, and some subunits form a complex homologous to condensin. The mechanism by which the DCC downregulates gene expression remains unclear. Here we show that the DCC controls the methylation state of lysine 20 of histone H4, leading to higher H4K20me1 and lower H4K20me3 levels on the X chromosomes of XX hermaphrodites relative to autosomes. We identify the PR-SET7 ortholog SET-1 and the Suv4-20 ortholog SET-4 as the major histone methyltransferases for monomethylation and di/trimethylation of H4K20, respectively, and provide evidence that X-chromosome enrichment of H4K20me1 involves inhibition of SET-4 activity on the X. RNAi knockdown of set-1 results in synthetic lethality with dosage compensation mutants and upregulation of X-linked gene expression, supporting a model whereby H4K20me1 functions with the condensin-like C. elegans DCC to repress transcription of X-linked genes. H4K20me1 is important for mitotic chromosome condensation in mammals, suggesting that increased H4K20me1 on the X may restrict access of the transcription machinery to X-linked genes via chromatin compaction.

In many animals, males have one X chromosome and females have two. However, the same amount of gene expression from X chromosomes is needed in the two sexes. The process of dosage compensation (DC) globally regulates X-chromosome gene expression to make it equal between the sexes, and it occurs in different ways in different animals. In mammals, one X chromosome in females is randomly inactivated, leaving one active X chromosome. In contrast, in the nematode worm C. elegans, the two X chromosomes in hermaphrodites are repressed two-fold to match gene expression to the single X chromosome in males. Previous work in C. elegans identified proteins required for DC that bind to the X chromosome, but their mode of action is not known. Here we show that DC proteins lead to higher levels of histone H4 lysine 20 monomethylation (H4K20me1) on hermaphrodite X chromosomes and that H4K20me1 functions in repressing X-chromosome gene expression. This shows that histone modification is an important aspect of the mechanism of dosage compensation. Together with previous work linking H4K20me1 to chromatin structure regulation, our results suggest that dosage compensation might lower gene expression on hermaphrodite X chromosomes by compacting them.

Nucleosomes suppress spontaneous mutations base-specifically in eukaryotes

It is unknown how the composition and structure of DNA within the cell affect spontaneous mutations. Theory suggests that in eukaryotic genomes, nucleosomal DNA undergoes fewer C→T mutations because of suppressed cytosine hydrolytic deamination relative to nucleosome-depleted DNA. Comparative genomic analyses and a mutation accumulation experiment showed that nucleosome occupancy nearly eliminated cytosine deamination, resulting in an ~50% decrease of the C→T mutation rate in nucleosomal DNA. Furthermore, the rates of G→T and A→T mutations were also about twofold suppressed by nucleosomes. On the basis of these results, we conclude that nucleosome-dependent mutation spectra affect eukaryotic genome structure and evolution and may have implications for understanding the origin of mutations in cancers and in induced pluripotent stem cells.

Promoters recognized by forkhead proteins exist for individual 21U-RNAs

C. elegans 21U-RNAs are equivalent to the piRNAs discovered in other metazoans and have important roles in gametogenesis and transposon control. The biogenesis and molecular function of 21U-RNAs and piRNAs are poorly understood. Here, we demonstrate that transcription of each 21U-RNA is regulated separately through a conserved upstream DNA motif. We use genomic analysis to show that this motif is associated with low nucleosome occupancy, a characteristic of many promoters that drive expression of protein-coding genes, and that RNA polymerase II is localized to this nucleosome-depleted region. We establish that the most conserved 8-mer sequence in the upstream region of 21U-RNAs, CTGTTTCA, is absolutely required for their individual expression. Furthermore, we demonstrate that the 8-mer is specifically recognized by Forkhead family (FKH) transcription factors and that 21U-RNA expression is diminished in several FKH mutants. Our results suggest that thousands of small noncoding transcription units are regulated by FKH proteins.

Manipulating nucleosome disfavoring sequences allows fine-tune regulation of gene expression in yeast

Understanding how precise control of gene expression is specified within regulatory DNA sequences is a key challenge with far-reaching implications. Many studies have focused on the regulatory role of transcription factor-binding sites. Here, we explore the transcriptional effects of different elements, nucleosome-disfavoring sequences and, specifically, poly(dA:dT) tracts that are highly prevalent in eukaryotic promoters. By measuring promoter activity for a large-scale promoter library, designed with systematic manipulations to the properties and spatial arrangement of poly(dA:dT) tracts, we show that these tracts significantly and causally affect transcription. We show that manipulating these elements offers a general genetic mechanism, applicable to promoters regulated by different transcription factors, for tuning expression in a predictable manner, with resolution that can be even finer than that attained by altering transcription factor sites. Overall, our results advance the understanding of the regulatory code and suggest a potential mechanism by which promoters yielding prespecified expression patterns can be designed.

Caenorhabditis elegans dosage compensation regulates histone H4 chromatin state on X chromosomes

Dosage compensation equalizes X-linked gene expression between the sexes. This process is achieved in Caenorhabditis elegans by hermaphrodite-specific, dosage compensation complex (DCC)-mediated, 2-fold X chromosome downregulation. How the DCC downregulates gene expression is not known. By analyzing the distribution of histone modifications in nuclei using quantitative fluorescence microscopy, we found that H4K16 acetylation (H4K16ac) is underrepresented and H4K20 monomethylation (H4K20me1) is enriched on hermaphrodite X chromosomes in a DCC-dependent manner. Depletion of H4K16ac also requires the conserved histone deacetylase SIR-2.1, while enrichment of H4K20me1 requires the activities of the histone methyltransferases SET-1 and SET-4. Our data suggest that the mechanism of dosage compensation in C. elegans involves redistribution of chromatin-modifying activities, leading to a depletion of H4K16ac and an enrichment of H4K20me1 on the X chromosomes. These results support conserved roles for histone H4 chromatin modification in worm dosage compensation analogous to those seen in flies, using similar elements and opposing strategies to achieve differential 2-fold changes in X-linked gene expression.

ArchTEx: accurate extraction and visualization of next-generation sequence data

MOTIVATION

The extension of mapped sequence tags is a common step in the analysis of single-end next-generation sequencing (NGS) data from protein localization and chromatin studies. The optimal extension can vary depending on experimental and technical conditions. Improper extension of sequence tags can obscure or mislead the interpretation of NGS results. We present an algorithm, ArchTEx (Architectural Tag Extender), which identifies the optimal extension of sequence tags based on the maximum correlation between forward and reverse tags and extracts and visualizes sites of interest using the predicted extension.

AVAILABILITY AND IMPLEMENTATION

ArchTEx requires Java 1.6 or newer. Source code and the compiled program are freely available at http://sourceforge.net/projects/archtex/

CONTACT

mjbuck@buffalo.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Regulation of DCC localization by HTZ-1/H2A.Z and DPY-30 does not correlate with H3K4 methylation levels

Dosage compensation is a specialized form of gene regulation that balances sex-chromosome linked gene expression between the sexes. In C. elegans, dosage compensation is achieved by the activity of the dosage compensation complex (DCC). The DCC binds along both X chromosomes in hermaphrodites to down-regulate gene expression by half, limiting X-linked gene products to levels produced in XO males. Sequence motifs enriched on the X chromosome play an important role in targeting the DCC to the X. However, these motifs are not strictly X-specific and therefore other factors, such as the chromatin environment of the X chromosome, are likely to aid in DCC targeting. Previously, we found that loss of HTZ-1 results in partial disruption of dosage compensation localization to the X chromosomes. We wanted to know whether other chromatin components coordinated with HTZ-1 to regulate DCC localization. One candidate is DPY-30, a protein known to play a role in DCC localization. DPY-30 homologs in yeast, flies, and mammals are highly conserved members of histone H3 lysine 4 (H3K4) methyltransferase Set1/MLL complexes. Therefore, we investigated the hypothesis that the dosage compensation function of DPY-30 involves H3K4 methylation. We found that in dpy-30 animals the DCC fails to stably bind chromatin. Interestingly, of all the C. elegans homologs of Set1/MLL complex subunits, only DPY-30 is required for stable DCC binding to chromatin. Additionally, loss of H3K4 methylation does not enhance DCC mislocalization in htz-1 animals. We conclude that DPY-30 and HTZ-1 have unique functions in DCC localization, both of which are largely independent of H3K4 methylation.

Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project

We systematically generated large-scale data sets to improve genome annotation for the nematode Caenorhabditis elegans, a key model organism. These data sets include transcriptome profiling across a developmental time course, genome-wide identification of transcription factor-binding sites, and maps of chromatin organization. From this, we created more complete and accurate gene models, including alternative splice forms and candidate noncoding RNAs. We constructed hierarchical networks of transcription factor-binding and microRNA interactions and discovered chromosomal locations bound by an unusually large number of transcription factors. Different patterns of chromatin composition and histone modification were revealed between chromosome arms and centers, with similarly prominent differences between autosomes and the X chromosome. Integrating data types, we built statistical models relating chromatin, transcription factor binding, and gene expression. Overall, our analyses ascribed putative functions to most of the conserved genome.