Genome-wide chromatin footprinting reveals changes in replication origin architecture induced by pre-RC assembly.

Start sites of DNA replication are marked by the origin recognition complex (ORC), which coordinates Mcm2-7 helicase loading to form the prereplicative complex (pre-RC). Although pre-RC assembly is well characterized in vitro, the process is poorly understood within the local chromatin environment surrounding replication origins. To reveal how the chromatin architecture modulates origin selection and activation, we "footprinted" nucleosomes, transcription factors, and replication proteins at multiple points during the Saccharomyces cerevisiae cell cycle. Our nucleotide-resolution protein occupancy profiles resolved a precise ORC-dependent footprint at 269 origins in G2. A separate class of inefficient origins exhibited protein occupancy only in G1, suggesting that stable ORC chromatin association in G2 is a determinant of origin efficiency. G1 nucleosome remodeling concomitant with pre-RC assembly expanded the origin nucleosome-free region and enhanced activation efficiency. Finally, the local chromatin environment restricts the loading of the Mcm2-7 double hexamer either upstream of or downstream from the ARS consensus sequence (ACS).

Comparative analysis of metazoan chromatin organization.

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

Dynamic loading and redistribution of the Mcm2-7 helicase complex through the cell cycle.

Eukaryotic replication origins are defined by the ORC-dependent loading of the Mcm2-7 helicase complex onto chromatin in G1. Paradoxically, there is a vast excess of Mcm2-7 relative to ORC assembled onto chromatin in G1. These excess Mcm2-7 complexes exhibit little co-localization with ORC or replication foci and can function as dormant origins. We dissected the mechanisms regulating the assembly and distribution of the Mcm2-7 complex in the Drosophila genome. We found that in the absence of cyclin E/Cdk2 activity, there was a 10-fold decrease in chromatin-associated Mcm2-7 relative to the levels found at the G1/S transition. The minimal amounts of Mcm2-7 loaded in the absence of cyclin E/Cdk2 activity were strictly localized to ORC binding sites. In contrast, cyclin E/Cdk2 activity was required for maximal loading of Mcm2-7 and a dramatic genome-wide reorganization of the distribution of Mcm2-7 that is shaped by active transcription. Thus, increasing cyclin E/Cdk2 activity over the course of G1 is not only critical for Mcm2-7 loading, but also for the distribution of the Mcm2-7 helicase prior to S-phase entry.

DNA replication and transcription programs respond to the same chromatin cues.

DNA replication is a dynamic process that occurs in a temporal order along each of the chromosomes. A consequence of the temporally coordinated activation of replication origins is the establishment of broad domains (>100 kb) that replicate either early or late in S phase. This partitioning of the genome into early and late replication domains is important for maintaining genome stability, gene dosage, and epigenetic inheritance; however, the molecular mechanisms that define and establish these domains are poorly understood. The modENCODE Project provided an opportunity to investigate the chromatin features that define the Drosophila replication timing program in multiple cell lines. The majority of early and late replicating domains in the Drosophila genome were static across all cell lines; however, a small subset of domains was dynamic and exhibited differences in replication timing between the cell lines. Both origin selection and activation contribute to defining the DNA replication program. Our results suggest that static early and late replicating domains were defined at the level of origin selection (ORC binding) and likely mediated by chromatin accessibility. In contrast, dynamic domains exhibited low ORC densities in both cell types, suggesting that origin activation and not origin selection governs the plasticity of the DNA replication program. Finally, we show that the male-specific early replication of the X chromosome is dependent on the dosage compensation complex (DCC), suggesting that the transcription and replication programs respond to the same chromatin cues. Specifically, MOF-mediated hyperacetylation of H4K16 on the X chromosome promotes both the up-regulation of male-specific transcription and origin activation.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans.

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.

Epigenome characterization at single base-pair resolution.

We have combined standard micrococcal nuclease (MNase) digestion of nuclei with a modified protocol for constructing paired-end DNA sequencing libraries to map both nucleosomes and subnucleosome-sized particles at single base-pair resolution throughout the budding yeast genome. We found that partially unwrapped nucleosomes and subnucleosome-sized particles can occupy the same position within a cell population, suggesting dynamic behavior. By varying the time of MNase digestion, we have been able to observe changes that reflect differential sensitivity of particles, including the eviction of nucleosomes. To characterize DNA-binding features of transcription factors, we plotted the length of each fragment versus its position in the genome, which defined the minimal protected region of each factor. This process led to the precise mapping of protected and exposed regions at and around binding sites, and also determination of the degree to which they are flanked by phased nucleosomes and subnucleosome-sized particles. Our protocol and mapping method provide a general strategy for epigenome characterization, including nucleosome phasing and dynamics, ATP-dependent nucleosome remodelers, and transcription factors, from a single-sequenced sample.

Nucleosomes influence multiple steps during replication initiation.

Eukaryotic replication origin licensing, activation and timing are influenced by chromatin but a mechanistic understanding is lacking. Using reconstituted nucleosomal DNA replication assays, we assessed the impact of nucleosomes on replication initiation. To generate distinct nucleosomal landscapes, different chromatin-remodeling enzymes (CREs) were used to remodel nucleosomes on origin-DNA templates. Nucleosomal organization influenced two steps of replication initiation: origin licensing and helicase activation. Origin licensing assays showed that local nucleosome positioning enhanced origin specificity and modulated helicase loading by influencing ORC DNA binding. Interestingly, SWI/SNF- and RSC-remodeled nucleosomes were permissive for origin licensing but showed reduced helicase activation. Specific CREs rescued replication of these templates if added prior to helicase activation, indicating a permissive chromatin state must be established during origin licensing to allow efficient origin activation. Our studies show nucleosomes directly modulate origin licensing and activation through distinct mechanisms and provide insights into the regulation of replication initiation by chromatin.

Noncoding Transcription Is a Driving Force for Nucleosome Instability in spt16 Mutant Cells.

FACT (facilitates chromatin transcription) consists of two essential subunits, Spt16 and Pob3, and functions as a histone chaperone. Mutation of spt16 results in a global loss of nucleosomes as well as aberrant transcription. Here, we show that the majority of nucleosome changes upon Spt16 depletion are alterations in nucleosome fuzziness and position shift. Most nucleosomal changes are suppressed by the inhibition of RNA polymerase II (Pol II) activity. Surprisingly, a small subgroup of nucleosome changes is resistant to transcriptional inhibition. Notably, Spt16 and distinct histone modifications are enriched at this subgroup of nucleosomes. We also report 1,037 Spt16-suppressed noncoding transcripts (SNTs) and found that the SNT start sites are enriched with the subgroup of nucleosomes resistant to Pol II inhibition. Finally, the nucleosomes at genes overlapping SNTs are more susceptible to changes upon Spt16 depletion than those without SNTs. Taken together, our results support a model in which Spt16 has a role in maintaining local nucleosome stability to inhibit initiation of SNT transcription, which once initiated drives additional nucleosome loss upon Spt16 depletion.

Rare codons capacitate Kras-driven de novo tumorigenesis.

The KRAS gene is commonly mutated in human cancers, rendering the encoded small GTPase constitutively active and oncogenic. This gene has the unusual feature of being enriched for rare codons, which limit protein expression. Here, to determine the effect of the rare codon bias of the KRAS gene on de novo tumorigenesis, we introduced synonymous mutations that converted rare codons into common codons in exon 3 of the Kras gene in mice. Compared with control animals, mice with at least 1 copy of this Kras(ex3op) allele had fewer tumors following carcinogen exposure, and this allele was mutated less often, with weaker oncogenic mutations in these tumors. This reduction in tumorigenesis was attributable to higher expression of the Kras(ex3op) allele, which induced growth arrest when oncogenic and exhibited tumor-suppressive activity when not mutated. Together, our data indicate that the inherent rare codon bias of KRAS plays an integral role in tumorigenesis.