Site-specific phasing in the chromatin of the rDNA in Dictyostelium discoideum

Analysis of chromatin organization by deep sequencing technologies

Micrococcal nuclease (MNase) is an endonuclease that cleaves native DNA at high frequency, but is blocked in chromatin by sites of intimate DNA-protein interaction, including nucleosomal regions. Protection from MNase cleavage has often been used to map transcription factor binding sites and nucleosomal positions on a single-gene basis; however, by combining MNase digestion with high--throughput, paired-end DNA sequencing, it is now possible to simultaneously map DNA-protein interaction regions across the entire genome. Biochemical and bioinformatic protocols are detailed for global mono-nucleosome positioning at ~160 bp spacing coverage, but are applicable to mapping more broadly or for site-specific binding of transcription factors at ~50 bp resolution.

DNA sequence-directed organization of chromatin: structure-based computational analysis of nucleosome-binding sequences

The folding of DNA on the nucleosome core particle governs many fundamental issues in eukaryotic molecular biology. In this study, an updated set of sequence-dependent empirical "energy" functions, derived from the structures of other protein-bound DNA molecules, is used to investigate the extent to which the architecture of nucleosomal DNA is dictated by its underlying sequence. The potentials are used to estimate the cost of deforming a collection of sequences known to bind or resist uptake in nucleosomes along various left-handed superhelical pathways and to deduce the features of sequence contributing to a particular structural form. The deformation scores reflect the choice of template, the deviations of structural parameters at each step of the nucleosome-bound DNA from their intrinsic values, and the sequence-dependent "deformability" of a given dimer. The correspondence between the computed scores and binding propensities points to a subtle interplay between DNA sequence and nucleosomal folding, e.g., sequences with periodically spaced pyrimidine-purine steps deform at low cost along a kinked template whereas sequences that resist deformation prefer a smoother spatial pathway. Successful prediction of the known settings of some of the best-resolved nucleosome-positioning sequences, however, requires a template with "kink-and-slide" steps like those found in high-resolution nucleosome structures.

Two types of telomeric chromatin in Tetrahymena thermophila

The telomeric d(GGGGTT).d(AACCCC) repeat tracts (G4T2 repeats) in Tetrahymena thermophila macronuclei were shown previously to be packaged in a non-nucleosomal DNA-protein complex. Here, we demonstrate that these telomeric repeats, together with a short region of the immediately adjacent non-telomeric sequence, exist in two distinct types of chromatin. The non-nucleosomal complex (type I complex) comprises approximately 90 to 97% of telomeric DNA, has no apparent underlying periodic nucleosomal substructure, and includes the whole telomeric tract as well as the immediately adjacent sequence. Type II chromatin, comprising the remaining approximately 3 to 10% of the total telomeric DNA, consists of tightly packed nucleosomes clustered at the inner border of the telomeric tracts, with a periodicity of 154(+/-3) bp. This packing is similar to that of telomeric nucleosomes in vertebrates. However, in contrast to the unstability of vertebrate telomeric mononucleosomes, the T. thermophila mononucleosomes were stable to micrococcal nuclease digestion. During the natural lengthening of the T. thermophila telomeric DNA tracts that occurs in vegetatively dividing cells, the overall ratio of type I and type II chromatin did not change. However, type I complex expanded with the length of the telomeric DNA repeat tract, and the number of telomeric nucleosomes increased from an average of one, up to three to four, per telomeric tract. This finding of telomeric nucleosomes in T. thermophila suggests that the difference between vertebrates and lower eukaryotes in telomeric chromatin structure is quantitative rather than qualitative. We propose that deposition of nucleosomes competes with non-nucleosomal complex formation on telomeric DNA, resulting in a sub-population of chimeric telomeres containing inner nucleosomes abutting a distal, variable length of type I complex.

The roles of telomeres and telomerase in cell life span

Telomeres cap and protect the ends of chromosomes from degradation and illegitimate recombination. The termini of a linear template cannot, however, be completely replicated by conventional DNA-dependent DNA polymerases, and thus in the absence of a mechanisms to counter this effect, telomeres of eukaryotic cells shorten every round of DNA replication. In humans and possibly other higher eukaryotes, telomere shortening may have been adopted to limit the life span of somatic cells. Human somatic cells have a finite proliferative capacity and enter a viable growth arrested state called senescence. Life span appears to be governed by cell division, not time. The regular loss of telomeric DNA could therefore serve as a mitotic clock in the senescence programme, counting cell divisions. In most eukaryotic organisms, however, telomere shortening can be countered by the de novo addition of telomeric repeats by the enzyme telomerase. Cells which are "immortal' such as the human germ line or tumour cell lines, established mouse cells, yeast and ciliates, all maintain a stable telomere length through the action of telomerase. Abolition of telomerase activity in such cells nevertheless results in telomere shortening, a process that eventually destabilizes the ends of chromosomes, leading to genomic instability and cell growth arrest or death. Therefore, loss of terminal DNA sequences may limit cell life span by two mechanisms: by acting as a mitotic clock and by denuding chromosomes of protective telomeric DNA necessary for cell viability.

Flexing and folding double helical DNA

DNA base sequence, once thought to be interesting only as a carrier of the genetic blueprint, is now recognized as playing a structural role in modulating the biological activity of genes. Primary sequences of nucleic acid bases describe real three-dimensional structures with properties reflecting those structures. Moreover, the structures are base sequence dependent with individual residues adopting characteristic spatial forms. As a consequence, the double helix can fold into tertiary arrangements, although the deformation is much more gradual and spread over a larger molecular scale than in proteins. As part of an effort to understand how local structural irregularities are translated at the macromolecular level in DNA and recognized by proteins, a series of calculations probing the structure and properties of the double helix have been performed. By combining several computational techniques, complementary information as well as a series of built-in checks and balances for assessing the significance of the findings are obtained. The known sequence dependent bending, twisting, and translation of simple dimeric fragments have been incorporated into computer models of long open DNAs of varying length and chemical composition as well as in closed double helical circles and loops. The extent to which the double helix can be forced to bend and twist is monitored with newly parameterized base sequence dependent elastic energy potentials based on the observed configurations of adjacent base pairs in the B-DNA crystallographic literature.

Nucleosomal organization of telomere-specific chromatin in rat

Rat liver interphase chromosomes have telomeres 20-100 kb in length. Micrococcal nuclease digestion of nuclei cleaves telomeres with a uniform 157 bp periodicity, producing soluble particles that sediment in sucrose gradients exactly like oligonucleosomes. The monomeric telomere particles comigrate with nucleosome core particles on nucleoprotein and DNA gels but do not bind H1. DNAase I cleaves telomere nucleoprotein into a series of bands spaced by about 10.4 bp and with the same intensity distribution as bands from bulk nucleosomes. Removal of H1 from chromatin alters the sedimentation properties of telomeres in parallel with bulk chromatin. Thus, telomeres of mammals are constructed of closely spaced nucleosomes, in contrast with the telomeres of lower eukaryotes, which show no evidence of nucleosomal structure.

Euplotes crassus has genes encoding telomere-binding proteins and telomere-binding protein homologs

We have identified two 1.6 kb macronuclear DNA molecules from Euplotes crassus that hybridize to the alpha subunit of the Oxytricha telomere protein. We have shown that one of these molecules encodes the 51 kDa Euplotes telomere protein while the other appears to encode a homolog of the telomere protein. Although this homolog clearly differs in sequence from the Euplotes telomere protein, the two proteins share extensive amino acid sequence identity with each other and with the alpha subunit of the Oxytricha telomere protein. In all three proteins 35-36% of the amino acids are identical, while 54-56% are similar. The most extended regions of sequence conservation map within the N-terminal section; this section has been shown to comprise the DNA-binding domain in the Euplotes telomere protein. Our findings suggest that some of the conserved amino acids may be involved in DNA recognition and binding. The gene encoding the telomere protein homolog contains two introns; one of these introns is only 24 bp in length. This is the smallest mRNA intron reported to date.

Saccharomyces telomeres assume a non-nucleosomal chromatin structure

The chromatin structures of the telomeric and subtelomeric regions on chromosomal DNA molecules in Saccharomyces cerevisiae were analyzed using micrococcal nuclease and DNAse I. The subtelomeric repeats X and Y' were assembled in nucleosomes. However, the terminal tracts of C1-3A repeats were protein protected in a particle larger than a nucleosome herein called a telosome. The proximal boundary of the telosome was a DNase I hypersensitive site. This boundary between the telosome and adjacent nucleosomes was completely accessible to Escherichia coli dam methylase when this enzyme was expressed in yeast, whereas a site 250 bp internal to the telomeric repeats was relatively inaccessible. Telosomes could be cleaved from chromosome ends with nuclease and solubilized as protein-DNA complexes. Immunoprecipitation of chromosomal telosomes with antiserum to the RAP1 protein indicated that RAP1 was one component of isolated telosomes. Thus, the termini of chromosomal DNA molecules in yeast are assembled in a non-nucleosomal structure encompassing the entire terminal C1-3A tract. This structure is separated from adjacent nucleosomes by a region of DNA that is highly accessible to enzymes.

The curvature vector in nucleosomal DNAs and theoretical prediction of nucleosome positioning

Using our model for predicting DNA superstructures from the sequence, the average distribution of the phases of curvature along the sequences of the set of the 177 nucleosomal DNAs investigated by Satchwell et al. (J. Mol. Biol. 191 (1986) 659) was calculated. The diagram obtained shows very significant features which allow the visualization of the intrinsic nucleosomal superstructure characterized by two quasi-parallel tracts of a flat left-handed superhelical turn connected by a left-handed inflection in a perpendicular direction; such a superstructure appears to be closely related to the nucleosome model of Travers and Klug (Phil. Trans. R. Soc. Lond. 317 (1987) 537). The nucleosomal curvature phase diagram was then adopted as a sensitive determinant for the nucleosome virtual positioning in DNAs via correlation function, obtaining a good agreement with the experimental mapping of SV40 regulatory region as recently investigated by Ambrose et al. (J. Mol. Biol. 209 (1989) 255). This analysis shows also the presence of a constant phase relation between the virtual nucleosome positions which suggests its possible implication in the nucleosome condensation in chromatin.

RAP1 protein interacts with yeast telomeres in vivo: overproduction alters telomere structure and decreases chromosome stability

The protein encoded by the RAP1 gene of S. cerevisiae binds in vitro to a consensus sequence occurring at a number of sites in the yeast genome, including the repeated sequence C2-3A(CA)1-6 found at yeast telomeres. We present two lines of evidence for the in vivo binding of RAP1 protein at telomeres: first, RAP1 is present in telomeric chromatin and second, alterations in the level of RAP1 protein affect telomere length. The length changes seen with under- and overexpression of RAP1 are consistent with the interpretation that RAP1 binding to telomeres protects them from degradation. Unexpectedly, overproduction of the RAP1 protein was also shown to decrease greatly chromosome stability, suggesting that RAP1 mediates interactions that have a more global effect on chromosome behavior than simply protecting telomeres from degradation. Such interactions may involve telomere associations both with other telomeres and/or with structural elements of the nucleus.

Position effect at S. cerevisiae telomeres: reversible repression of Pol II transcription

S. cerevisiae chromosomes end with the telomeric repeat (TG1-3)n. When any of four Pol II genes was placed immediately adjacent to the telomeric repeats, expression of the gene was reversibly repressed as demonstrated by phenotype and mRNA analyses. For example, cells bearing a telomere-linked copy of ADE2 produced predominantly red colonies (a phenotype characteristic of ade2- cells) containing white sectors (characteristic of ADE2+ cells). Repression was due to proximity to the telomere itself since an 81 bp tract of (TG1-3)n positioned downstream of URA3 when URA3 was approximately 20 kb from the end of chromosome VII did not alter expression of the gene. However, this internal tract of (TG1-3)n could spontaneously become telomeric, in which case expression of the URA3 gene was repressed. These data demonstrate that yeast telomeres exert a position effect on the transcription of nearby genes, an effect that is under epigenetic control.

The upstream limit of nuclease-sensitive chromatin in Dictyostelium rRNA genes neighbors a topoisomerase I-like cluster

In the chromatin of Dictyostelium ribosomal RNA (rRNA) genes, the coding and upstream flanking regions are sensitive to endonucleases. This sensitivity stops about 2.3 x 10(3) bases upstream from the transcription start, at a point we call the structural boundary. Upstream from the boundary an 850 base-pair region is strongly protected against micrococcal nuclease cleavage, particularly in rapidly transcribing vegetative cells, and upstream from this the pattern of nuclease protection suggests that positioned nucleosomes are present. On the gene side of the structural boundary nucleosomes are known to be absent in vegetative cells but present in differentiating slug cells where the rRNA synthesis rate is lower. We show that in slugs these nucleosomes are randomly distributed, in contrast to those upstream from the boundary. Close to the gene side of the boundary is a duplication of the putative promoter located 29 base-pairs distant from four clustered topoisomerase I recognition sequences, which are cleaved by endogenous topoisomerase I-like activity. An additional topoisomerase I recognition sequence found upstream from the structural boundary is not cleaved in chromatin. The possible significance of these sequences and structures in transcription is discussed.

Statistical positioning of nucleosomes by specific protein-binding to an upstream activating sequence in yeast

Arrays of nucleosomes were positioned with respect to the GAL1-GAL10 intergenic region inserted into Saccharomyces cerevisiae minichromosomes. Deletions of DNA flanking the upstream activation sequence left the array unaltered, showing that nucleosome positioning was not a consequence of sequence-specific histone-DNA interactions but depended on proximity to the galactose-responsive upstream activation sequence (UASG). Replacement of the upstream activation sequence by synthetic oligonucleotides with different protein-binding properties identified a short sequence within this region that is responsible for the ordered array. This sequence overlaps a binding site for GAL4 protein, a positive regulator of transcription, but exerts its effect on chromatin structure independently of GAL4, probably through binding a novel factor that is not GAL-specific.

Statistical distributions of nucleosomes: nonrandom locations by a stochastic mechanism

Expressions are derived for distributions of nucleosomes in chromatin. Nucleosomes are placed on DNA at the densities found in bulk chromatin, and their locations are allowed to vary at random. No further assumptions are required to simulate the periodic patterns of digestion obtained with various nucleases. The introduction of a boundary constraint, due for example to sequence-specific protein binding, results in an array of regularly spaced nucleosomes at nonrandom locations, similar to the arrays reported for some genes and other chromosomal regions.

Topoisomerase I cleavage sites identified and mapped in the chromatin of Dictyostelium ribosomal RNA genes

Sites of an endogenous activity that has the properties of a DNA topoisomerase I have been identified on the palindromic ribosomal RNA genes of the slime mould Dictyostelium discoideum. This was done in vitro, by treating isolated nuclei with sodium dodecyl sulphate, which denatures topoisomerase during its cycle of nicking, strand passing and resealing, and hence reveals the DNA cleavages. It was also done in vivo using the drug camptothecin, which is believed to stabilize the cleavable complex of topoisomerase I plus DNA, hence increasing the chances of cleavage when sodium dodecyl sulphate is subsequently added. The cleavages in vitro and in vivo were mapped by indirect end-labelling. Both treatments cause what appear to be strong double-stranded cleavages at 200 and 2200 base-pairs and at 17 X 10(3) base-pairs upstream from the rRNA transcription start. The cleavage at 200 base-pairs was analysed in greater detail using RNA hybridization probes specific for single DNA strands. The cleavage is in fact composed of three closely spaced nicks on each DNA strand. The DNA sequence at each of the nicks is strongly homologous across 15 base-pairs. Sodium dodecyl sulphate-induced cleavage by eukaryotic topoisomerase I is known to yield enzyme covalently attached to the 3' cut end of the DNA. We show that protein-linked DNA restriction fragments with their 3' ends at the cleavage sites are selectively retarded on denaturing gels, which provides strong evidence that the unusual cluster of cleavages is caused by a topoisomerase I. Additionally, the camptothecin results revealed cleavages not only at the specific upstream sites, but also across the transcribed region. Interestingly, the zone of camptothecin-assisted cleavage does not extend as far at the 3' end of the gene as the zone of endogenous nuclease sensitivity.

Precise nucleosome positioning in the promoter of the chicken beta A globin gene

Histone octamers were reconstituted onto 5' end-labelled DNA fragments derived from the promoter region of the chicken beta A globin gene. The location of the reconstituted histone octamer with respect to the DNA sequence of each fragment was assessed by Exonuclease III digestion of purified nucleosome monomers. By this approach we have found a strong preference for histone octamers to be positioned over nucleotides -206 to -62 relative to the gene cap site. This stretch of DNA contains all those 5' beta globin sequences which, by DNase footprinting, bind specific protein factors and incorporates three promoter consensus sequence motifs. The upstream terminal 32 base pairs of this DNA segment contains the binding sites for the erythrocyte specific G-string binding protein and transcription factor Spl and appears to be relatively weakly bound to the histone octamer.

Telomere proteins: specific recognition and protection of the natural termini of Oxytricha macronuclear DNA

The macronuclear DNA in the ciliated protozoan O. nova consists of integral of 10(7) gene-sized DNA molecules, all of which terminate with 20 bp of C4A4 repeats followed by a 3' (G4T4)2 single-stranded tail. Two immunologically distinct proteins of 55 and 26 kd, which are tenaciously, but noncovalently associated with Oxytricha macronuclear DNA termini, have been purified. These proteins protect DNA termini from degradation by the exonuclease Bal31. They also facilitate retention of natural and synthetic telomeric DNAs onto nitrocellulose. The Oxytricha proteins are not simply C4A4-binding proteins. Rather, their efficient binding requires both the 3' single-stranded (G4T4)2 tail and the adjacent duplex region. Thus, these proteins require both the sequence and the structure of natural DNA termini for efficient binding. As such they represent the first described example of telomeric-specific proteins.

Mapping of endogenous nuclease-sensitive regions and of putative topoisomerase sites of action along the chromatin of Dictyostelium ribosomal RNA genes

Indirect end-labelling and the digestion patterns of endogenous and exogenous nucleases were used to analyse chromatin organization along the ribosomal RNA genes of Dictyostelium discoideum cells. A zone just upstream from the 5' end of the coding region was particularly sensitive to endogenous nucleases. In exponentially growing cells, this hypersensitive zone extended from -350 to -1600 bp relative to the transcription start. In sharp contrast, the DNA between 0 and -350 bp was strongly protected. In differentiating cells, in which the ribosomal RNA transcription rate is low, the 5' hypersensitive zone was more diffuse than in exponentially growing cells, and the protected region at the 5' end of the transcribed region was less pronounced. It is known that where DNA topoisomerase is acting on DNA, the addition of sodium dodecyl sulphate will result in cleavage of the DNA and covalent attachment of the enzyme to the cut DNA end. Treatment of nuclei from both exponentially growing cells and differentiating cells with SDS caused double-stranded cleavages at -200 (i.e. within the protected region), at -2200, and at two sites at about -17 kb. A fraction of the cleavage products appeared to be strongly associated with protein. Novobiocin, a DNA topoisomerase II inhibitor, did not inhibit the SDS-induced cleavages in vegetative cells. However, it significantly reduced the extent of nuclease cleavage within the -350 to -1600 bp hypersensitive zone. The possibility is discussed that there are two DNA topoisomerase-like activities on the ribosomal genes. One is site-specific and novobiocin-insensitive. We speculate that the other is responsible for maintaining DNA at the 5' end of the gene in a torsionally strained, nuclease-hypersensitive state.

Methidiumpropyl-EDTA-iron(II) cleavage of ribosomal DNA chromatin from Dictyostelium discoideum

We have used methidiumpropyl-EDTA-iron(II) [MPE.Fe(II)] in parallel with micrococcal nuclease to investigate the chromatin structure of the extrachromosomal palindrome ribosomal RNA genes of Dictyostelium. Confirming our earlier results with micrococcal nuclease (1,2), MPE.Fe(II) digested the coding region of rapidly transcribing rRNA genes as a smear, indicating the absence or severe disruption of nucleosomes, whereas in slowly transcribing rRNA genes, a nucleosomal ladder was produced. In the central non-transcribed spacer region of the palindrome, MPE.Fe(II) digestion resulted in a normal nucleosomal repeat, whereas micrococcal nuclease gave a complex banding pattern. The difference is attributed to the lower sequence specificity of MPE.Fe(II) compared to micrococcal nuclease. In the terminal region of the palindrome, however, both substances gave a complex chromatin digestion pattern. In this region the DNA appears to be packaged in structures strongly positioned with respect to the underlying DNA sequence.

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.

DNA bending and its relation to nucleosome positioning

X-ray and solution studies have shown that the conformation of a DNA double helix depends strongly on its base sequence. Here we show that certain sequence-dependent modulations in structure appear to determine the rotational positioning of DNA about the nucleosome. Three different experiments are described. First, a piece of DNA of defined sequence (169 base-pairs long) is closed into a circle, and its structure examined by digestion with DNAase I: the helix adopts a highly preferred configuration, with short runs of (A, T) facing in and runs of (G, C) facing out. Secondly, the same sequence is reconstituted with a histone octamer: the angular orientation around the histone core remains conserved, apart from a small uniform increase in helix twist. Finally, it is shown that the average sequence content of DNA molecules isolated from chicken nucleosome cores is non-random, as in a reconstituted nucleosome: short runs of (A, T) are preferentially positioned with minor grooves facing in, while runs of (G, C) tend to have their minor grooves facing out. The periodicity of this modulation in sequence content (10.17 base-pairs) corresponds to the helix twist in a local frame of reference (a result that bears on the change in linking number upon nucleosome formation). The determinants of translational positioning have not been identified, but one possibility is that long runs of homopolymer (dA) X (dT) or (dG) X (dC) will be excluded from the central region of the supercoil on account of their resistance to curvature.