DNAase I, DNAase II and staphylococcal nuclease cut at different, yet symmetrically located, sites in the nucleosome core

A practical guide for DNase-seq data analysis: from data management to common applications

Deoxyribonuclease I (DNase I)-hypersensitive site sequencing (DNase-seq) has been widely used to determine chromatin accessibility and its underlying regulatory lexicon. However, exploring DNase-seq data requires sophisticated downstream bioinformatics analyses. In this study, we first review computational methods for all of the major steps in DNase-seq data analysis, including experimental design, quality control, read alignment, peak calling, annotation of cis-regulatory elements, genomic footprinting and visualization. The challenges associated with each step are highlighted. Next, we provide a practical guideline and a computational pipeline for DNase-seq data analysis by integrating some of these tools. We also discuss the competing techniques and the potential applications of this pipeline for the analysis of analogous experimental data. Finally, we discuss the integration of DNase-seq with other functional genomics techniques.

Nick-seq for single-nucleotide resolution genomic maps of DNA modifications and damage

DNA damage and epigenetic marks are well established to have profound influences on genome stability and cell phenotype, yet there are few technologies to obtain high-resolution genomic maps of the many types of chemical modifications of DNA. Here we present Nick-seq for quantitative, sensitive, and accurate mapping of DNA modifications at single-nucleotide resolution across genomes. Pre-existing breaks are first blocked and DNA modifications are then converted enzymatically or chemically to strand-breaks for both 3'-extension by nick-translation to produce nuclease-resistant oligonucleotides and 3'-terminal transferase tailing. Following library preparation and next generation sequencing, the complementary datasets are mined with a custom workflow to increase sensitivity, specificity and accuracy of the map. The utility of Nick-seq is demonstrated with genomic maps of site-specific endonuclease strand-breaks in purified DNA from Eschericia coli, phosphorothioate epigenetics in Salmonella enterica Cerro 87, and oxidation-induced abasic sites in DNA from E. coli treated with a sublethal dose of hydrogen peroxide. Nick-seq applicability is demonstrated with strategies for >25 types of DNA modification and damage.

Genetic and epigenetic determinants establish a continuum of Hsf1 occupancy and activity across the yeast genome

Heat shock factor 1 is the master transcriptional regulator of molecular chaperones and binds to the same cis-acting heat shock element (HSE) across the eukaryotic lineage. In budding yeast, Hsf1 drives the transcription of ∼20 genes essential to maintain proteostasis under basal conditions, yet its specific targets and extent of inducible binding during heat shock remain unclear. Here we combine Hsf1 chromatin immunoprecipitation sequencing (seq), nascent RNA-seq, and Hsf1 nuclear depletion to quantify Hsf1 binding and transcription across the yeast genome. We find that Hsf1 binds 74 loci during acute heat shock, and these are linked to 46 genes with strong Hsf1-dependent expression. Notably, Hsf1's induced DNA binding leads to a disproportionate (∼7.5-fold) increase in nascent transcription. Promoters with high basal Hsf1 occupancy have nucleosome-depleted regions due to the presence of "pioneer factors." These accessible sites are likely critical for Hsf1 occupancy as the activator is incapable of binding HSEs within a stably positioned, reconstituted nucleosome. In response to heat shock, however, Hsf1 accesses nucleosomal sites and promotes chromatin disassembly in concert with the Remodels Structure of Chromatin (RSC) complex. Our data suggest that the interplay between nucleosome positioning, HSE strength, and active Hsf1 levels allows cells to precisely tune expression of the proteostasis network.

DNase-seq predicts regions of rotational nucleosome stability across diverse human cell types

DNase-seq is primarily used to identify nucleosome-depleted DNase I hypersensitive (DHS) sites genome-wide that correspond to active regulatory elements. However, ∼40 yr ago it was demonstrated that DNase I also digests with a ∼10-bp periodicity around nucleosomes matching the exposure of the DNA minor groove as it wraps around histones. Here, we use DNase-seq data from 49 samples representing diverse cell types to reveal this digestion pattern at individual loci and predict genomic locations where nucleosome rotational positioning, the orientation of DNA with respect to the histone surface, is stably maintained. We call these regions DNase I annotated regions of nucleosome stability (DARNS). Compared to MNase-seq experiments, we show DARNS correspond well to annotated nucleosomes. Interestingly, many DARNS are positioned over only one side of annotated nucleosomes, suggesting that the periodic digestion pattern attenuates over the nucleosome dyad. DARNS reproduce the arrangement of nucleosomes around transcription start sites and are depleted at ubiquitous DHS sites. We also generated DARNS from multiple lymphoblast cell line (LCL) samples. We found that LCL DARNS were enriched at DHS sites present in most of the original 49 samples but absent in LCLs, while multi-cell-type DARNS were enriched at LCL-specific DHS sites. This indicates that variably open DHS sites are often occupied by rotationally stable nucleosomes in cell types where the DHS site is closed. DARNS provide additional information about precise DNA orientation within individual nucleosomes not available from other nucleosome positioning assays and contribute to understanding the role of chromatin in gene regulation.

Nucleosome depletion at yeast terminators is not intrinsic and can occur by a transcriptional mechanism linked to 3'-end formation

Genome-wide mapping of nucleosomes generated by micrococcal nuclease (MNase) suggests that yeast promoter and terminator regions are very depleted of nucleosomes, predominantly because their DNA sequences intrinsically disfavor nucleosome formation. However, MNase has strong DNA sequence specificity that favors cleavage at promoters and terminators and accounts for some of the correlation between occupancy patterns of nucleosomes assembled in vivo and in vitro. Using an improved method for measuring nucleosome occupancy in vivo that does not involve MNase, we confirm that promoter regions are strongly depleted of nucleosomes, but find that terminator regions are much less depleted than expected. Unlike at promoter regions, nucleosome occupancy at terminators is strongly correlated with the orientation of and distance to adjacent genes. In addition, nucleosome occupancy at terminators is strongly affected by growth conditions, indicating that it is not primarily determined by intrinsic histone-DNA interactions. Rapid removal of RNA polymerase II (pol II) causes increased nucleosome occupancy at terminators, strongly suggesting a transcription-based mechanism of nucleosome depletion. However, the distinct behavior of terminator regions and their corresponding coding regions suggests that nucleosome depletion at terminators is not simply associated with passage of pol II, but rather involves a distinct mechanism linked to 3'-end formation.

Small DNA pieces in C. elegans are intermediates of DNA fragmentation during apoptosis

While studying small noncoding RNA in C. elegans, we discovered that protocols used for isolation of RNA are contaminated with small DNA pieces. After electrophoresis on a denaturing gel, the DNA fragments appear as a ladder of bands, ∼10 nucleotides apart, mimicking the pattern of nuclease digestion of DNA wrapped around a nucleosome. Here we show that the small DNA pieces are products of the DNA fragmentation that occurs during apoptosis, and correspondingly, are absent in mutant strains incapable of apoptosis. In contrast, the small DNA pieces are present in strains defective for the engulfment process of apoptosis, suggesting they are produced in the dying cell prior to engulfment. While the small DNA pieces are also present in a number of strains with mutations in predicted nucleases, they are undetectable in strains containing mutations in nuc-1, which encodes a DNase II endonuclease. We find that the small DNA pieces can be labeled with terminal deoxynucleotidyltransferase only after phosphatase treatment, as expected if they are products of DNase II cleavage, which generates a 3′ phosphate. Our studies reveal a previously unknown intermediate in the process of apoptotic DNA fragmentation and thus bring us closer to defining this important pathway.

Role of chromatin states in transcriptional memory

Establishment of cellular memory and its faithful propagation is critical for successful development of multicellular organisms. As pluripotent cells differentiate, choices in cell fate are inherited and maintained by their progeny throughout the lifetime of the organism. A major factor in this process is the epigenetic inheritance of specific transcriptional states or transcriptional memory. In this review, we discuss chromatin transitions and mechanisms by which they are inherited by subsequent generations. We also discuss illuminating cases of cellular memory in budding yeast and evaluate whether transcriptional memory in yeast is nuclear or cytoplasmically inherited.

High-resolution mapping and characterization of open chromatin across the genome

Mapping DNase I hypersensitive (HS) sites is an accurate method of identifying the location of genetic regulatory elements, including promoters, enhancers, silencers, insulators, and locus control regions. We employed high-throughput sequencing and whole-genome tiled array strategies to identify DNase I HS sites within human primary CD4+ T cells. Combining these two technologies, we have created a comprehensive and accurate genome-wide open chromatin map. Surprisingly, only 16%-21% of the identified 94,925 DNase I HS sites are found in promoters or first exons of known genes, but nearly half of the most open sites are in these regions. In conjunction with expression, motif, and chromatin immunoprecipitation data, we find evidence of cell-type-specific characteristics, including the ability to identify transcription start sites and locations of different chromatin marks utilized in these cells. In addition, and unexpectedly, our analyses have uncovered detailed features of nucleosome structure.

Gene splice sites correlate with nucleosome positions

Gene sequences in the vicinity of splice sites are found to possess dinucleotide periodicities, especially RR and YY, with the period close to the pitch of nucleosome DNA. This confirms previously reported findings about preferential positioning of splice junctions within the nucleosomes. The RR and YY dinucleotides oscillate counter-phase, i.e., their respective preferred positions are shifted about half-period from one another, as it was observed earlier for AA and TT dinucleotides. Species specificity of nucleosome positioning DNA pattern is indicated by the predominant use of the periodical GG(CC) dinucleotides in human and mouse genes, as opposed to predominant AA(TT) dinucleotides in Arabidopsis and C. elegans.

Redefinition of the cleavage sites of DNase I on the nucleosome core particle

DNase I has been widely used for the footprinting of DNA-protein interactions including analyses of nucleosome core particle (NCP) structure. Our understanding of the relationship between the footprint and the structure of the nucleosome complex comes mainly from digestion studies of NCPs, since they have a well-defined quasi-symmetrical structure and have been widely investigated. However, several recent results suggest that the established consensus of opinion regarding the mode of digestion of NCPs by DNase I may be based on erroneous interpretation of results concerning the relationship between the NCP ends and the dyad axis. Here, we have used reconstituted NCPs with defined ends, bulk NCPs prepared with micrococcal nuclease and molecular modelling to reassess the mode of DNase I digestion. Our results indicate that DNase I cuts the two strands of the nucleosomal DNA independently with an average stagger of 4 nt with the 3'-ends protruding. The previously accepted value of 2 nt stagger is explained by the finding that micrococcal nuclease produces NCPs not with flush ends, but with approximately 1 nt 5'-recessed ends. Furthermore we explain why the DNA stagger is an even and not an odd number of nucleotides. These results are important for studies using DNase I to probe nucleosome structure in complex with other proteins or any DNA-protein complex containing B-form DNA. We also determine the origin of the 10n +/- 5 nt periodicity found in the internucleosomal ladder of DNase I digests of chromatin from various species. The explanation of the 10n +/- 5 nt ladder may have implications for the structure of the 30 nm fibre.

Early necrotic DNA degradation: presence of blunt-ended DNA breaks, 3' and 5' overhangs in apoptosis, but only 5' overhangs in early necrosis

The structure of DNA breaks in early necrosis was analyzed and compared with apoptotic DNA degradation using in vivo and cell culture models. Early necrosis (1 hour after cell death) was produced in vivo by the freezing-thawing of rat thymus and in cell culture of Jurkat cells. Apoptosis was induced in the same cell types using dexamethasone for thymus and staurosporine for Jurkat cells. Selective detection of double-strand DNA breaks with blunt ends was performed by in situ ligation. Blunt-ended breaks bearing 5' phosphates were detected in apoptotic but not in early necrotic cells. Pretreatment of apoptotic and necrotic tissue with Klenow enzyme with or without added dNTPs reduced all 3' or 5' overhangs to blunt ends. Subsequent in situ ligation with blunt-ended probes revealed no 3' overhangs in necrotic cells. However double-strand cuts with 5' overhangs were abundant in necrotic DNA. 5' Overhangs were also detected in apoptotic cells. Presence of exclusively 5' overhangs in early necrosis with absence of a variety of possible DNA ends, suggests the existence of a specific orderly mechanism of DNA degradation.

Nuclear apoptosis detection by flow cytometry: influence of endogenous endonucleases

Nuclear apoptosis is characterized by chromatin condensation and progressive DNA cleavage into high-molecular-weight fragments and oligonucleosomes. These complex phenomena can be mediated by the activation of a multiplicity of enzymes, characterized by specific patterns of cation dependance, pH requirement, and mode of activation. The significance of this multiplicity of enzymes that cleave genomic DNA has been attributed to the need of death effector pathways specific for cell types/tissues, the level of cell differenciation, and the nature of the apoptotic stimuli. The activation of these factors contributes to the development of alterations that can be detected specifically by flow cytometric assays, namely, propidium iodide assays, acridine orange/ethidium bromide double staining, the TUNEL and ISNT techniques, and the assays of DNA sensitivity to denaturation. Although applicable to a wide spectrum of cell types, an increasing body of literature indicates that these techniques cannot be universally applied to all cell lines and apoptotic conditions: The requirement of a particular mediator(s) of nuclear apoptosis or the absence of endonuclease activity can limit the relevance of certain techniques. Finally, endonucleases recruited during primary necrosis can introduce nuclear alterations detected by some assays and raise the problem of their specificity. This review underlines the need for strategies to accurately detect and quantify nuclear apoptosis by flow cytometry when new cell systems and apoptotic conditions are considered.

Cerebral ischemia produces laddered DNA fragments distinct from cardiac ischemia and archetypal apoptosis

The electrophoretic pattern of laddered DNA fragments which has been observed after cerebral ischemia is considered to indicate that neurons are dying by apoptosis. Herein the authors directly demonstrate using ligation-mediated polymerase chain reaction methods that 99% of the DNA fragments produced after either global or focal ischemia in adult rats, or produced after hypoxia-ischemia in neonatal rats, have staggered ends with a 3' recess of approximately 8 to 10 nucleotides. This is in contrast to archetypal apoptosis in which the DNA fragments are blunt ended as seen during developmental programmed cell death in dying cortical neurons, neuroblastoma, or thymic lymphocytes. It is not simply ischemia that results in staggered ends in DNA fragments because ischemic myocardium is similar to archetypal apoptosis with a vast majority of blunt-ended fragments. It is concluded that the endonucleases that produce this staggered fragmentation of the DNA backbone in ischemic brain must be different than those of classic or type I apoptosis.

Evidence for involvement of tumor necrosis factor-alpha in apoptotic death of bone marrow cells in myelodysplastic syndromes

We previously reported excessive apoptosis and high levels of tumor necrosis factor-alpha (TNF-alpha) in the bone marrows of patients with myelodysplastic syndromes (MDS), using histochemical techniques. The present studies provide further circumstantial evidence for the involvement of TNF-alpha in apoptotic death of the marrow cells in MDS. Using our newly developed in situ double-labeling technique that sequentially employs DNA polymerase (DNA Pol) followed by terminal deoxynucleotidyl transferase (TdT) to label cells undergoing apoptosis, we have characterized DNA fragmentation patterns during spontaneous apoptosis in MDS bone marrow and in HL60 cells treated with TNF-alpha or etoposide (VP16). Clear DNA laddering detected by gel electrophoresis in MDS samples confirmed the unique length of apoptotic DNA fragments (180-200 bp). Surprisingly, however, phenotypically heterogeneous population of MDS cells as well as the homogenous population of HL60 cells showed three distinct labeling patterns after double labeling--only DNA-Pol reaction, only TdT reaction, and a combined DNA Pol + TdT reaction, albeit in different cohorts of cells. Each labeling pattern was found at all morphological stages of apoptosis. MDS mononuclear cells, during spontaneous apoptosis in 4 hr cultures, showed highest increase in double-labeled cells (DNA Pol + TdT reaction). Interestingly, this was paralleled by TNF-alpha-induced apoptosis in HL60 cells. In contrast, VP16 treatment of HL60 cells led to increased apoptosis in cells showing only TdT reaction. The double-labeling technique was applied to normal bone marrow and peripheral blood mononuclear cells after treatment with known endonucleases that specifically cause 3' recessed (BamHI), 5' recessed (PstI), or blunt ended (DraI) double-stranded DNA breaks. It was found that the DNA-Pol reaction in MDS and HL60 cells corresponds to 3' recessed DNA fragments, the TdT reaction to 5' recessed and/or blunt ended fragments, and a combined "DNA Pol + TdT reaction" corresponds to a copresence of 3' recessed with 5' recessed and/or blunt ended fragments. Clearly, therefore, apoptotic DNA fragments, in spite of a unique length, may have differently staggered ends that could be cell (or tissue) specific and be selectively triggered by different inducers of apoptosis. The presence of TNF-alpha-inducible apoptotic DNA fragmentation pattern in MDS supports its involvement in these disorders and suggests that anti-TNF-alpha (or anticytokine) therapy may be of special benefit to MDS patients, where no definitive treatment is yet available.

Protective nucleosome centering at splice sites as suggested by sequence-directed mapping of the nucleosomes

The characteristic AA(TT) sequence pattern of the nucleosome DNA derived earlier is used for prediction of nucleosome positions around splice junctions of eukaryotic genes. Two large datasets (2000 sequences each) were collected consisting of DNA segments with the exon/intron and intron/exon splice junctions, from various eukaryotic species. Positions of predicted nucleosomes near the junction sites were calculated. Those junctions which are found to belong to the nucleosomes, are located preferentially within a few base pairs from the midpoint of the nucleosome DNA. That is, obligatory GT- and AG-ends of the introns are more frequently located near the nucleosome dyad axis, within the best protected middle 10-15 base pairs of the nucleosome DNA. In addition, a tendency is observed for the strongest nucleosomes to form more often in the introns, in accordance with the hypothesis on the chromatin-organizing role of introns.

Presence of double-strand breaks with single-base 3' overhangs in cells undergoing apoptosis but not necrosis

Apoptotic cells in rat thymus were labeled in situ in paraffin-embedded and frozen tissue sections by ligation of double-stranded DNA fragments containing digoxigenin or Texas red. Two forms of double-stranded DNA fragments were prepared using the polymerase chain reaction: one was synthesized using Taq polymerase, which yields products with single-base 3' overhangs, and one using Pfu polymerase, which produces blunt-ended products. Both types of fragment could be ligated to apoptotic nuclei in thymus, indicating the presence in such nuclei of DNA double-strand breaks with single-base 3' overhangs as well as blunt ends. However, in nuclei with DNA damage resulting from a variety of nonapoptotic processes (necrosis, in vitro autolysis, peroxide damage, and heating) single-base 3' overhangs were either nondetectable or present at much lower concentrations than in apoptotic cells. Blunt DNA ends were present in such tissues, but at lower concentrations than in apoptotic cells. In contrast, in all of these forms of DNA damage, nuclei contained abundant 3'-hydroxyls accessible to labeling with terminal deoxynucleotidyl transferase. Thus, although single-base 3' overhangs and blunt ends are present in apoptotic nuclei, the specificity of the in situ ligation of 3'-overhang fragments to apoptotic nuclei indicates that apoptotic cells labeled in this way can readily be distinguished from cells with nonapoptotic DNA damage. These data are consistent with the involvement of an endonuclease similar to DNase I in apoptosis, which is predicted to leave short 3' overhangs as well as blunt ends in digestion of chromatin.

Multitude of inverted repeats characterizes a class of anchorage sites of chromatin loops to the nuclear matrix

In order to understand the nature of DNA sequences that organize chromatin into domains or loops, we have cloned the nuclear matrix DNA (1.7% of the total DNA) from human myelogenous leukemia cells in culture. Nuclear matrix is formed by interactions between specific stretches of DNA of about 0.1 to 5.0 kb with protein transcription factors, nuclear enzymes, and structural proteins. Nuclear matrix is believed to be the exclusive nuclear microenvironment in which initiation of DNA replication, transcription, and repair take place. The matrix attachment regions (MARs) of DNA have transcriptional enhancer activity, harbor the origins of replication of the human genome, and define the borders between neighboring chromatin loops. In this study we report the sequence of the human MAR fragment 19.2 of a size of 542 bp. Hum. MAR 19.2 is composed of TG-, CA-, CT-, and GA-rich blocks and shows 8 perfect and imperfect inverted repeats. Thus, we have identified a novel class of MARs with sequence characteristics divergent from the AT-rich class of MARs. The inverted repeats of the 19.2 sequence might be stabilized into their cruciform configuration by torsional strain and by specific transcription/replication protein factors. This MAR might function in the initiation of replication of the flanking chromatin domain and in the regulation of the transcriptional activity of the gene(s) that reside in this domain.

Rearrangements of the nucleosome structure in chromatin by poly(ADP-ribose)

In order to approach and clarify the effect of poly(ADP-ribose) on the nucleosomal structure, polynucleosomes from calf thymus were incubated with long poly(ADP-ribose) chains prepared in vitro and examined by ELISA with antibodies directed against the five individual histones H1, H2A, H2B, H3 and H4 as well as against two synthetic peptides in residues 1-25 of H2B and 130-135 of H3. The results showed that: (i) free ADP-ribose polymers did indeed interact with the nucleosomes; (ii) the accessibility of epitopes recognized by any of the different antibodies was altered, the binding of antibodies being increased or decreased depending on the quantity of poly(ADP-ribose) added thereby suggesting a modulation in nucleosome structure; (iii) for any ADP-ribose polymer concentration, core histones as well as histone H1 were always recognized by their respective antibodies, thus suggesting that poly(ADP-ribose) does not seem to cause complete stripping of histones from nucleosomal DNA.

Mobility of positioned nucleosomes on 5 S rDNA

We report on a dynamic aspect of nucleosome positioning, in the absence of transcription-related events, on sea urchin 5 S rDNA. On tandem repeats of nucleosome length DNA of this strongly positioning sequence, histone octamers assemble in one dominant position surrounded by minor positions, ten base-pairs apart and therefore with identical rotational setting of the DNA coil. The existence of this cluster of positions, determined using micrococcal nuclease is confirmed by the results from DNase I footprinting and restriction enzyme analysis. The results from these techniques and from two-dimensional nucleoprotein polyacrylamide gel analysis indicate that the cluster of octamer positions is in dynamic equilibrium, in low ionic conditions, suggesting that the minor positions reflect fluctuations around the major nucleosome site. Histone octamer mobility appears to be temperature dependent and is reversibly inhibited by Mg2+.

Alternative DNA-protein interactions in variable-length internucleosomal regions associated with Drosophila Adh distal promoter expression

Chromatin at the Drosophila Adh distal promoter displays an ordered but different conformation in different cell types as detected by a modified exonuclease protection assay and accessibility to endonucleases. In cells not transcribing Adh (ADH-) sequences between -40 to +30 of the distal RNA initiation site exist as a DNA linker between positioned nucleosomes, and appear to interact with a specific DNA-binding protein. In contrast, a longer linker DNA, from -140 to +30, is bound in a multi-protein transcription initiation complex in cells that specifically transcribe the distal (adult) ADH RNA (ADH+A). These DNA-protein interactions can account for a localized open chromatin structure at the distal promoter in ADH+A cells. The observed mutually exclusive patterns of DNA-protein interactions in the linkers of different ADH cell types between -40 to +30 suggest a model for organizing alternative chromatin structure associated with gene regulation. Two DNA binding proteins, one being a TATA box binding factor, compete for overlapping sites to allow either assembly of a transcription initiation complex and transcription, or positioning of nucleosomes for stable repression.

DNA strand breaks alter histone ADP-ribosylation

Histone ADP-ribosylation was studied using two-dimensional gel electrophoresis after cleavage of the nuclear DNA with nucleases. Modified histones carrying different numbers of ADP-ribose groups form a ladder of bands above each variant histone. Cellular lysates containing unfragmented DNA mainly synthesize mono(ADP-ribosylated) histones. Cleavage of the DNA with either DNase I or micrococcal nuclease to fragments of an average size of 10-20 kilobases (kb) dramatically induces the formation of poly(ADP-ribosylated) species of histones in nuclei. As the number of DNA strand breaks produced by either DNase I or micrococcal nuclease increases and a great number of DNA cuts is introduced (fragments of 0.4-0.2 kb), the size of the poly(ADP-ribose) chains on the histones decreases. Finally, in the presence of 10 mM cAMP as an inhibitor of poly(ADP-ribose) glycohydrolase, human lymphoid nuclei synthesize hyper(ADP-ribosylated) histone H2B with at least 40 ADP-ribose groups attached to it. Lateral ladders emanating at precise points of the linear ladder on hypermodified H2B can arise from branching of poly(ADP-ribose) or from multiple monomodifications of glutamic (or aspartic) acid residues. Branching or de novo monomodifications occur after a precise number of ADP-ribose groups have been added to a histone molecule. Poly(ADP-ribosylated) histones thus appear to be intermediates in nuclear processes involving DNA strand breaks.

Protein/DNA architecture of the DNase I hypersensitive region of the Drosophila hsp26 promoter

Genomic footprinting on the Drosophila hsp26 promoter in isolated nuclei has shown that a TATA box binding factor is present before and after induction by heat shock, while three of the seven heat shock consensus sequences 5' of the gene are occupied (presumably by heat shock factor, HSF) specifically on heat shock. The sites of HSF interaction are separated by greater than 200 bp of which approximately 150 bp are bound to the surface of a nucleosome. The juxtaposition of these various macromolecules on the DNA suggests a basis for the major DNase I hypersensitive site 5' of hsp26 and a novel tertiary structure for the promoter complex.

Sequence-specific positioning of core histones on an 860 base-pair DNA. Experiment and theory

Previous experiments have shown that the locations of the histone octamer on DNA molecules of 140 to 240 base-pairs (bp) are influenced strongly by the nucleotide sequence. Here we have studied the locations of the histone octamer on a relatively long DNA molecule of 860 bp, using two different nucleases, micrococcal and DNAase I. Data were obtained from both the protein--DNA complexes and from the naked DNA at single-bond resolution, and then were analyzed by densitometry to yield plots of differential cleavage, which show clearly the changes in cutting due to the addition of protein. Our results show that the placement of core histones on the 860 bp molecule is definitely non-random. The digestion data provide evidence for five nucleosome cores, the centers of which lie in defined locations. In all but one of these protein--DNA complexes, the DNA adopts a unique, highly preferred rotational setting with respect to the protein surface. Another protein--DNA complex is unusual in that it protects 200 bp from digestion, yet is cut in its very center as if it were split into two parts. The apparent average twist of the DNA within all of these protein--DNA complexes is 10.2(+/- 0.1) bp, as measured by the periodicity of DNAase I digestion. This value is in excellent agreement with the twist of 10.21(+/- 0.05) bp deduced from the periodicity of sequence content in chicken nucleosome core DNA. In addition, we observe a discontinuity in the periodic cutting by DNAase I of about -1 to -3 bonds in going from any nucleosome core to the next. The most plausible interpretation of this discontinuity is that it reflects the angle by which adjacent protein--DNA complexes are aligned. Thus, any nucleosome may be related to its neighbor by a left-handed rotation in space of -1/10.2 to -3/10.2 helix turns, or -35 degrees to -105 degrees. Repeated many times, this operation would build a long, left-handed helix of nucleosomes similar to that described by many workers for the packing of nucleosomes in chromatin. In order to look for any long-range influences on the positioning of the histone octamer in the 860 bp molecule (as would be expected if the nucleosomes have to fit into some higher-order structure), we have examined the locations of the histone octamer on five different isolated short fragments of the 860-mer, all of nucleosomal length.(ABSTRACT TRUNCATED AT 400 WORDS)

Ion-induced DNA structure change in nucleosomes

Physical methods have been used to study calcium binding to the nucleosome core particle. Equilibrium dialysis of Ca2+ and spectroscopic analysis of a Ca2+ analogue show that the ion binds tightly to the particles, resulting in a significant change of DNA circular dichroism. This suggests that base stacking may be altered as a result of Ca2+ binding. In the presence of Ca2+, the absorbance and fluorescence properties of methylene blue (MB), a DNA-specific intercalator, confirm that the dye binds tightly to nucleosomes by intercalation. However, secondary changes occur which suggest that the MB binding site is altered as a result of Ca2+ binding. Triplet state anisotropy decay and triplet lifetime quenching both show that in the Ca2+-nucleosome complex, methylene blue is capable of wobbling over a substantial angular range at its binding site. To explain these data, it is proposed that Ca2+ binding to nucleosomes causes DNA to fold by means of a series of sharp bends (kinks). The properties of bound MB are best explained if it is presumed that the intercalator binds tightly to such kinked sites in the nucleosome. On the basis of these observations, we discuss the possibility that multivalent ion concentration in the nucleus is high enough that the smooth to kinked helix equilibrium may be near to its midpoint. Near such a midpoint, the secondary structure of DNA in the nucleosome might prove to be sensitive to effector molecule binding and to site-specific variation of DNA or histone composition within genes.

Histones and their modifications

Histones constitute the protein core around which DNA is coiled to form the basic structural unit of the chromosome known as the nucleosome. Because of the large amount of new histone needed during chromosome replication, the synthesis of histone and DNA is regulated in a complex manner. During RNA transcription and DNA replication, the basic nucleosomal structure as well as interactions between nucleosomes must be greatly altered to allow access to the appropriate enzymes and factors. The presence of extensive and varied post-translational modifications to the otherwise highly conserved histone primary sequences provides obvious opportunities for such structural alterations, but despite concentrated and sustained effort, causal connections between histone modifications and nucleosomal functions are not yet elucidated.

Crystallographic structure of the octameric histone core of the nucleosome at a resolution of 3.3 A

The structure of the (H2A-H2B-H3-H4)2 histone octamer has been determined by means of x-ray crystallographic techniques at a resolution of 3.3 angstroms. The octamer is a prolate ellipsoid 110 angstroms long and 65 to 70 angstroms in diameter, and its general shape is that of a rugby ball. The size and shape are radically different from those determined in earlier studies. The most striking feature of the histone octamer is its tripartite organization, that is, a central (H3-H4)2 tetramer flanked by two H2A-H2B dimers. The DNA helix, placed around the octamer in a path suggested by the features on the surface of the protein, appears like a spring holding the H2A-H2B dimers at either end of the (H3-H4)2 tetramer.

DNase I sensitivity of nuclear DNA measured by flow cytometry

The DNase I digestion kinetics of DNA in isolated nuclei (from HeLa or murine mammary carcinoma, 67 cells) were assayed flow cytometrically by measuring the changes in ethidium bromide (EtBr) fluorescence following various digestion time intervals. The DNase I digestion curve was characterized by an initial 25-30% increase in fluorescence upon addition of the enzyme, a rapid reduction in fluorescence to approximately 50-55% in 30 minutes, and a limit digest of 45-50% beyond 45 minutes. Throughout digestion, the DNA histogram retained its characteristic bimodal shape, showing that histogram rearrangement was not responsible for the changes in EtBr fluorescence. Irradiation with 5 X 10(6) rads (137Cs-gamma-rays) or exposure to 50 mM EDTA caused an increase in EtBr fluorescence similar to that caused by DNase I, suggesting that DNA nicking and/or chromatin loosening were responsible for this increase. Residual DNA assayed by the solubilization of 14C-TdR (thymidine)-labeled DNA indicated a similar kinetic pattern without the initial increase. However, at the limit digest, the fraction of DNA remaining trichloroacetic acid (TCA) insoluble (10%) was smaller than that measured by loss of EtBr fluorescence (50% of initial, 40% of maximum). Part of this difference was due to the presence of TCA soluble DNA trapped within the nuclear matrix (15-20%). This trapped DNA was released when the digested nuclei were exposed to 0.5-1.0 M NaCl just prior to EtBr staining. Exposure of HeLa cells to three agents that are believed to cause changes in chromatin structure resulted in alterations in the DNase I digestion kinetics measured flow cytometrically.(ABSTRACT TRUNCATED AT 250 WORDS)

Chromatin assembly in Xenopus oocytes: in vivo studies

We have followed the time course of chromatin assembly, DNA supercoiling, and transcription on a Xenopus 5S RNA gene clone injected into germinal vesicles of Xenopus oocytes. During the first 2 hr after DNA injection, there is a gradual enhancement in transcription that correlates with the increase in superhelical density of the DNA template; on the other hand, nucleosome assembly is already completed by 10-30 min after DNA injection. To probe further the DNA structure in the assembled minichromosomes, we injected enzymes and chemicals into the germinal vesicle. DNAase I and topoisomerase I injections reveal that the circular DNA has been assembled into two discrete and equally abundant types of chromatin: one type, which we call "dynamic" chromatin, is torsionally strained and is thus fully relaxed by those two enzymes. The other type, which we call "static" chromatin, still yields supercoiled DNA molecules after deproteinization. The dynamic chromatin is also relaxed by injection of novobiocin, and simultaneously, 5S RNA transcription is turned off. The results of our in vivo experiments suggest that the dynamic chromatin is the one that is transcriptionally active. We discuss the biological relevance of these findings.

Analysis of the chromatin assembled in germinal vesicles of Xenopus oocytes

We have injected circular DNA, labeled with 32P at a single restriction site, into germinal vesicles of Xenopus laevis oocytes in order to study the nucleosome arrangement on the assembled minichromosomes. Two types of genes were used in these studies, the somatic 5 S RNA gene unit of Xenopus borealis and the histone gene unit of Drosophila melanogaster. We find that injections of labeled DNA alone, at 1 ng DNA per oocyte, results in irregularly spaced nucleosomes and partially supercoiled DNA molecules. However, perfectly spaced nucleosomes are assembled and fully supercoiled DNA is recovered if 5 to 20 nanograms of cold vector DNA is coinjected with the labeled DNA. At the optimum chromatin assembly conditions, the nucleosomes are perfectly spaced with a 180 base-pair periodicity, but they are randomly positioned on the DNA. The assembly of a periodic chromatin structure is accompanied by a dramatic enhancement in the expression of the injected 5 S RNA gene.

Location of the primary sites of micrococcal nuclease cleavage on the nucleosome core

The positions and relative frequencies of the primary cleavages made by micrococcal nuclease on the DNA of nucleosome core particles have been found by fractionating the double-stranded products of digestion and examining their single-stranded compositions. This approach overcomes the problems caused by secondary events such as the exonucleolytic and pseudo-double-stranded actions of the nuclease and, combined with the use of high resolution gel electrophoresis, enables the cutting site positions to be determined with a higher precision than has been achieved hitherto. The micrococcal nuclease primary cleavage sites lie close (on average, within 0.5 nucleotide) to those previously determined by Lutter (1981) for the nucleases DNase I and DNase II. These similarities show that the accessible regions are the same for all three nucleases, the cleavage sites being dictated by the structure of the nucleosome core. The differences in the final products of the digestion are explained in terms of secondary cleavage events of micrococcal nuclease. While the strongly protected regions of the nucleosome core DNA are common to all three nucleases, there are differences in the relative degrees of cutting at the more exposed sites characteristic of the particular enzyme. In particular, micrococcal nuclease shows a marked polarity in the 3'-5' direction in the cutting rates as plotted along a single strand of the nucleosomal DNA. This is explained in terms of the three-dimensional structure of the nucleosome where, in any accessible region of the double helix, the innermost strand is shielded by the outermost strand on the one side and the histone core on the other. The final part of the paper is concerned with the preference of micrococcal nuclease to cleave at (A,T) sequences in chromatin.

The chromatin structure of an actively expressed, single copy yeast gene

When the yeast galactokinase gene is not active (repressed, not expressed, quiescent), there is an exceptionally regular nucleosome array on coding sequence galactokinase chromatin, as shown by both denaturing and non-denaturing gel analysis of staphylococcal nuclease digests. Expression of the gene results in a limited smearing of the nucleosome repeat peaks and an increase in interpeak DNA, appearing as a regular ladder of DNA bands on denaturing gels. On non-denaturing gels the pattern is more complex and molecular weight dependent. These data suggest an increase in intracore particle DNA accessibility, allowing staphylococcal nuclease to digest throughout the nucleosome in expressed chromatin. Comparison to bulk chromatin and to an operationally inactive gene (35S rDNA) show that the alteration is specific to expressed chromatin. In contrast, DNase I shows no differences in the digestion of the gene specific chromatin in expressed or inactive states.

A simple procedure for parallel sequence analysis of both strands of 5'-labeled DNA

Ligation of a 5'-labeled DNA restriction fragment results in a circular DNA molecule carrying the two 32Ps at the reformed restriction site. Double digestions of the circular DNA with the original enzyme and a second restriction enzyme cleavage near the labeled site allows direct chemical sequencing of one 5'-labeled DNA strand. Similar double digestions, using an isoschizomer that cleaves differently at the 32P-labeled site, allows direct sequencing of the now 3'-labeled complementary DNA strand. It is possible to directly sequence both strands of cloned DNA inserts by using the above protocol and a multiple cloning site vector that provides the necessary restriction sites. The simultaneous and parallel visualization of both DNA strands eliminates sequence ambiguities. In addition, the labeled circular molecules are particularly useful for single-hit DNA cleavage studies and DNA footprint analysis. As an example, we show here an analysis of the micrococcal nuclease-induced breaks on the two strands of the somatic 5S RNA gene of Xenopus borealis, which suggests that the enzyme may recognize and cleave small AT-containing palindromes along the DNA helix.

Another potential artifact in the study of nucleosome phasing by chromatin digestion with micrococcal nuclease

We show that, contrary to expectations, restriction enzyme cleavage of chicken erythrocyte nucleosome core particle DNA generates a series of distinct subnucleosome fragments. These fragments do not result from bulk nucleosome phasing in vivo, but arise from micrococcal nuclease cleavages internal to the core particle, at roughly 10-base pair intervals and at AT-rich sequences. Those 145-base pair DNA fragments remaining intact are a biased population in which the guanine content can fluctuate by as much as 10%, with a 10-base pair period. We suggest that these same considerations, when applied to a unique DNA sequence, are the true explanation for several previous claims for nucleosome phasing.

Simultaneous analysis of conformation and transcription of A and B groups of vitellogenin genes in male and female Xenopus during primary and secondary activation by estrogen

In male Xenopus, primary estradiol administration results in noncoordinate activation in the liver of the A and B groups of vitellogenin genes, both as judged by transcription and DNase I sensitivity in isolated nuclei, B group genes being activated preferentially in the first 20 hr. Secondary induction in males or "primary" induction in females results in a coordinate and equal transcription of these two groups of genes. The elevated transcriptional activity following primary estrogen stimulation returns to low levels rapidly but the high DNase I sensitivity of these genes persists for 2-3 months. A non-coordinate activation of the A and B groups of vitellogenin genes is however re-established in response to a second administration of estradiol 8 months after primary stimulation of male Xenopus.

Crystallization of the tetramer of histones H3 and H4

Crystals of the histone tetramer (H3-H4)2 from calf thymus have been grown. The crystals yield x-ray diffraction patterns with Bragg spacings as small as 3.5 angstroms. Crystals grown from two types of preparations have the symmetry of the space group P61 (or P65). The best crystals were grown from histones that had the amino terminal arms removed by mild trypsinization.

Correlation of nitrosourea murine bone marrow toxicity with deoxyribonucleic acid alkylation and chromatin binding sites

All of the clinically available nitrosourea antitumor agents produce serious treatment-limiting bone marrow toxicity. A reduction in this toxicity can be achieved by attaching the chloroethylnitrosourea cytotoxic group to C2 (chlorozotocin) or C1 (1-(2-chloroethyl)-3-(beta-D-glucopyranosyl)-1-nitrosourea, GANU) of glucose. Both glucose analogs are less myelotoxic in mice than 1-(2-chloroethyl)-3-cyclohepyl-1-nitrosourea (CCNU) or 1-(4-amino-2-methylpyrimidin-5-yl)methyl-3-(2-chloroethyl)-3-nitrosourea (ACNU), while retaining comparable antitumor activity against the murine L1210 leukemia. To define the nuclear mechanisms for this reduced myelotoxicity, alkylation of L1210 and murine bone marrow DNA was quantitated. With the use of the endonuclease micrococcal nuclease and DNase I, the sites of alkylation within the chromatin substructure were determined. Experiments were performed on L1210 leukemia or bone marrow cells that had been incubated in vitro for 2 hr with 0.1 mM [14C]chloroethyl drug. The quantitative alkylation of DNA by GANU was 1.3-fold greater in L1210, as compared to bone marrow, cells. This ratio of DNA alkylation is comparable to the 1.3 ratio we previously reported for chlorozotocin [L. C. Panasci, D. Green and P. S. Schein, J. clin. Invest. 64, 1103 (1979)]. In contrast, the ratio of alkylation (L1210:bone marrow DNA) for the myelotoxic ACNU was 0.66, similar to 0.59 for CCNU. Nuclease digestion experiments demonstrated that chlorozotocin and GANU preferentially alkylated internucleosomal linker regions of bone marrow chromatin, while nucleosome core particles were the preferred targets of CCNU and ACNU. The reduced myelotoxicity of chlorozotocin and GANU may be correlated with the advantageous ratio of L1210:bone marrow DNA alkylation and preferential alkylation of internucleosomal regions of bone marrow chromatin.

Effect of the B--Z transition in poly(dG-m5dC) . poly(dG-m5dC) on nucleosome formation

We have studied the properties of complexes formed between histones and the methylated synthetic polydeoxynucleotide poly(dG-m5dC). poly(dG-m5dC). This polymer undergoes the transition from B DNA to left-handed Z DNA at moderate ionic strength. When the polymer is in the Z form it will bind histones, but nucleosomes are not detected. When the polymer in the B form is combined with equimolar quantities of the four core histones and digested with micrococcal nuclease, particles are formed which behave in all respects as normal nucleosome cores. When these core particles are placed in solvents that would result in conversion of the protein-free polymer to the Z form, no transition is observed. The formation of a nucleosome core particle thus stabilizes the B form, whereas the presence of the Z form prevents nucleosome formation. The results suggest that if Z DNA is present in eukaryotic nuclei, it will serve to disrupt the normal chromatin structure.

Responses to androgens of rat ventral prostate nuclear androgen-binding sites sensitive and resistant to micrococcal nuclease

Rat ventral prostate nuclei were separated into three major fractions by mild digestion with micrococcal nuclease and two fractions by extensive digestion. All fractions contained androgen-binding sites. Almost 50% of nuclear binding sites were resistant to enzymic digestion when only 5-15% of total DNA was resistant. Under milder digestion conditions, 21% of nuclear binding sites were associated with an intermediate fraction, representing 16% of total nuclear DNA, which was enriched in specific androgen-regulated gene sequences. This fraction was rapidly degraded by more extensive digestion. The nuclease sensitivity of these particular genes was markedly influenced by castration and the administration of dihydrotestosterone to castrated animals. The nuclear content of both nuclease-resistant and -sensitive androgen-binding sites was decreased by castration. Whereas the administration of androgen to animals castrated 1 day previously preferentially replenished nuclease-resistant sites, nuclease-sensitive sites, including those associated with transcriptionally active regions, had apparent priority when androgen was supplied to animals castrated 7 days previously. The significance of these observations to the regulation of nuclear processes and the possible interrelationships of nuclease-sensitive and -insensitive sites are discussed.

Involvement of histone H1 in the structure of the linker DNA in nucleosomes as revealed by nucleases

This report describes experiments designed to study the organization of the linker DNA in nucleosomes. When rat liver nucleosomes (145 to 188 base pairs in length) were digested by Exonuclease III and then by nuclease S1 a series of bands on sizes 90 - 102 - 112 - 125 - 135 - 142 - 154 - 166 - 172 - 181 - bases was observed in denaturing electrophoretic gels. Digestion of H1-depleted nucleosomes under the same conditions results in a series of products of sizes (10.4) xn in base (n less than or equal to 14) only. This result is interpreted as reflecting a particular arrangement of linker DNA under the influence of histone H1.

Core particle stability critically depends upon a small number of terminal nucleotides

An apparently homogeneous population of core particles is in fact composed of three subpopulations which behave differently when exposed to a high concentration of ethidium bromide or to 0.6 M NaCl. These subspecies have been identified by the use of several techniques, viz., electron microscopy, sedimentation velocity and circular dichroism. The electrophoretic analysis of their DNA leads to the conclusion that core particle stability critically depends upon a small number of terminal nucleotides.

DNase II digestion of the nucleosome core: precise locations and relative exposures of sites

The precise locations and relative exposures of the DNase II-accessible sites in the nucleosome core DNA are determined using techniques previously employed for the enzyme DNase I. It is found that there are a number of similarities between the site exposure patterns for the two enzymes but that in general the DNase II seems to discriminate less among adjacent sites' accessibilities than does DNase I. The two enzymes attack essentially the same positions in the DNA, the average difference between the precise location of the site being less than one-half base for the two enzymes. Such close similarities in the digestion patterns of two enzymes with such different mechanisms of scission show that the patterns reflect the structure of the nucleosome core and not merely the properties of the particular enzyme used.

The helical periodicity of DNA on the nucleosome

The precise number of base pairs per turn of the DNA double helix in the nucleosome core particle has been the subject of controversy. In this paper the positions of nuclease cutting sites are analysed in three dimensions. Using this midpoint of the DNA on the nucleosome dyad as origin, the cutting site locations measured along a strand of DNA are mapped onto models of the nucleosome core containing DNA of different helical periodicities. It is found that a helical periodicity of 10.5 base pairs per turn leads to cutting site positions which are sterically inaccessible. In contrast, a periodicity of 10.0 base pairs per turn leads to cutting site positions which are not only sterically sound, but which fall into a pattern such as would be expected when the access of the nuclease to the DNA is restricted by the presence of the histone core on one side and of the adjacent superhelical turn of DNA on the other. As proposed earlier by us (1), a value for the helical periodicity close to 10 base pairs per turn on the nucleosome, taken together with a periodicity close to 10.5 for DNA in solution - a value now established - resolves the so-called linkage number paradox.

Contacts between DNA gyrase and its binding site on DNA: features of symmetry and asymmetry revealed by protection from nucleases

DNA gyrase supercoils DNA by passing one DNA segment through another by means of a reversible double-strand break at specific DNA sites. We determined the nucleotide sequence of two highly preferred gyrase binding sites and analyzed the grip of gyrase on the DNA by using protection from nuclease attack. The DNA-breakage site of gyrase was centered in about 50 base pairs (bp) of DNA that was completely protected from DNase I and flanked by DNA regions cut at average intervals of 9.9 bases. The same pattern of protection from DNase I was observed with topoisomerase II', an enzyme that shares structural homology with gyrase. The gyrase site of DNA breakage was off-center in the 140 bp of DNA protected from exonuclease III digestion. ATP or inhibitors of gyrase had little specific effect on DNase I protection. On addition of a nonhydrolyzable analogue of ATP, previously stable barriers to exonuclease III were invaded and new barriers appeared. We discuss a detailed model uniting these results with previous data on gyrase structure and mechanism.

Spatial organization of histones and DNA in the nucleosome core particle: a model

We present here an attempt to build up a space-filling model of the nucleosome core particle based on the chemical crosslinking data of Mirzabekov and co-workers (23). It is shown that the models proposed earlier are inconsistent with the results of these authors. The main characteristics of our model are as follows: a) the DNA superhelix contains at least 90 base pairs (bp) per turn; b) the particle has a dyad axis of symmetry; c) the histone octamer may be regarded as consisting of two heterotypic tetramers. The possible shape and function of core histones are discussed in the light of the model.