Molecular Biology of the Cell Sign up for new MBC in Press e-TOCs!

Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ostashevsky, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ostashevsky, J.

Vol. 9, Issue 11, 3031-3040, November 1998

A Polymer Model for the Structural Organization of Chromatin Loops and Minibands in Interphase Chromosomes

Joseph Ostashevsky

Department of Radiation Oncology, State University of New York, Health Science Center at Brooklyn, Brooklyn, New York 11203

Submitted May 15, 1998; Accepted August 21 1998
Monitoring Editor: Joseph Gall

    ABSTRACT
Top
Abstract
Introduction
Results & Discussion
Conclusions
References

A quantitative model of interphase chromosome higher-order structure is presented based on the isochore model of the genome and results obtained in the field of copolymer research. G1 chromosomes are approximated in the model as multiblock copolymers of the 30-nm chromatin fiber, which alternately contain two types of 0.5- to 1-Mbp blocks (R and G minibands) differing in GC content and DNA-bound proteins. A G1 chromosome forms a single-chain string of loop clusters (micelles), with each loop ~1-2 Mbp in size. The number of ~20 loops per micelle was estimated from the dependence of geometrical versus genomic distances between two points on a G1 chromosome. The greater degree of chromatin extension in R versus G minibands and a difference in the replication time for these minibands (early S phase for R versus late S phase for G) are explained in this model as a result of the location of R minibands at micelle cores and G minibands at loop apices. The estimated number of micelles per nucleus is close to the observed number of replication clusters at the onset of S phase. A relationship between chromosomal and nuclear sizes for several types of higher eukaryotic cells (insects, plants, and mammals) is well described through the micelle structure of interphase chromosomes. For yeast cells, this relationship is described by a linear coil configuration of chromosomes.

    INTRODUCTION
Top
Abstract
Introduction
Results & Discussion
Conclusions
References

The higher-order structure of interphase chromosomes is still poorly understood. Many models in the literature include a loop structure as one of the high levels of packing of a chromatin fiber in the nucleus (for review see van Holde, 1989; Wolffe, 1995). Earlier studies have suggested that an average chromatin loop contains ~50-100 kbp DNA (van Holde, 1989; Wolffe, 1995), whereas more recent studies suggest large loops of ~1-3 Mbp, which may include 50- to 100-kbp loops (Razin and Gromova, 1995; Sachs et al., 1995; Yokota et al., 1995; Johnston et al., 1997).

On the scale of 1-3 Mbp, high-resolution mapping of replication bands in S phase (Drouin et al., 1990, 1994) is similar to a quasiperiodic pattern of G (for dark in Giemsa) and R (for reverse, light in Giemsa) minibands observed in prophase chromosomes (Bak et al., 1981; Yunis, 1981). G minibands are AT rich, late replicating, and gene poor, whereas R minibands are GC and gene rich and early replicating and contain a less compact chromatin than do G minibands (Holmquist, 1992; Craig and Bickmore, 1993; Yokota et al., 1997). Isochores, long DNA segments having a size range from 0.2 to 1.3 Mbp with an excess of one type of nucleotide (e.g., AT rich or GC rich), are found in the genome of higher eukaryotes (Bernardi, 1995).

An important feature of G1 phase chromosomes is that they behave approximately as ideal Gaussian chains, which obey random-walk statistics (van den Engh et al., 1992; Sachs et al., 1995; Yokota et al., 1995). This was concluded from the proportionality between the mean square geometrical distance between two points on the chromosome and their genomic distance, on the scale up to ~1 Mbp. On a larger scale (up to 200 Mbp), this dependence has a much shallower slope than the initial one (Sachs et al., 1995; Yokota et al., 1995).

On the nuclear level, G1 chromosomes tend to occupy exclusive territories rather than overlapping extensively (Haaf and Schmid, 1991; Cremer et al., 1993; Zirbel et al., 1993; van Driel et al., 1995; Kurz et al., 1996; Ferreira et al., 1997; Zink et al., 1998). There is some contradiction between the random-walk behavior of chromatin and the discreteness of chromosomal domains: random coils do not have clear boundaries, and they are prone to overlap (de Gennes 1979; Grosberg and Khokhlov, 1994).

The size of a nucleus influences the compactness of individual chromosomes (Yokota et al., 1995, 1997; Sanchez et al., 1997). Compartmentalization of nuclear space is characteristic for chromosome functions (Spector, 1993; Strouboulis and Wolffe, 1996). In particular, DNA replication starts only in several hundred clusters per nucleus in early S phase (for review see Berezney et al., 1995a; Jackson and Cook, 1995).

Several polymer approaches to chromosomes exist in the literature (van den Engh et al., 1992; Hahnfeldt et al., 1993; Ostashevsky and Lange, 1994; Sikorav and Jannink, 1994; Duplantier et al., 1995; Sachs et al., 1995; Jannink et al., 1996; Ostashevsky, 1996, 1998; Houchmandzadeh et al., 1997; Liu and Sachs, 1997; Marko and Siggia, 1997a); however, only a few articles consider interphase chromosomes.

This study develops a model of the higher-order structure of interphase chromosomes that deals with the problems and takes into account the main facts mentioned above. In this model, based on the isochore model of the genome (Bernardi, 1995) and results obtained in the field of copolymer research (e.g., see Semenov et al., 1995, 1996), a G1 chromosome is approximated as a multiblock copolymer containing two types of blocks differing in GC content.

    RESULTS AND DISCUSSION
Top
Abstract
Introduction
Results & Discussion
Conclusions
References

The Model's Background

The presented model of interphase chromosomes is based on the following assumptions.

1) A mammalian G1 chromosome can be approximated as a multiblock copolymer alternately containing two types of polymer blocks different in GC content (Figure 1). This assumption is supported by the observation that the DNA sequence of high eukaryotes is not random but is a mosaic of isochores, which are long DNA segments (0.2-1.3 Mbp) with an excess of one type of nucleotides (AT or GC) (Bernardi, 1995). Although five families of isochores can be defined in mammalian genomes, the division of polymer blocks in two classes, R (GC rich) and G (AT rich), as made in the presented model, can be considered as a first approximation. It is argued below that the R and G blocks in the model are related to the interphase and prophase R and G minibands, which are ~1 Mbp in size (Bak et al., 1981; Yunis, 1981; Ronne et al., 1995); thus, the terms blocks and minibands will be used interchangeably in this article.


View larger version (10K):
[in this window]
[in a new window]
 
Figure 1.   A mammalian G1 chromosome is approximated as a multiblock copolymer containing two types of polymer blocks with different GC contents. Light and dark chromatin segments are R and G blocks, respectively.

2) A multiblock copolymer containing two alternately located types of blocks can form a single-chain string of loop clusters called micelles (Halperin, 1991). A micelle consists of a certain number of loops, the termini of which, formed by blocks of one type, are located in close proximity to each other (Figure 2). Micelle structures are well studied for diblock copolymers (polymer chains having only two blocks) and ionomers (polymer chains with charged groups at the ends) (e.g., see Semenov et al., 1995, 1996). Large multiblock copolymers form single-chain micelles, and small diblock copolymers form multichain micelles. Formation of loops and organization of them in micelles constitute an entropically unfavorable process, because the number of possible polymer conformations decreases, but it occurs in multiblock copolymers because of the energically favorable processes of repulsion between unlike monomer units and/or attraction between like monomer units (de Gennes 1979; Grosberg and Khokhlov, 1994). For a multiblock copolymer in aqueous solution, which contains two types of blocks with hydrophobic and hydrophilic groups, hydrophobic blocks form loop termini at the micelle cores, and hydrophilic blocks are located at loop apices.


View larger version (17K):
[in this window]
[in a new window]
 
Figure 2.   Schematic drawing of a micelle in a G1 chromosome. Dark and light chromatin blocks are G and R minibands, respectively. Segments drawn with free ends can be either chromosome ends or intermicelle links. The circle represents the micelle core (see inset), where loop termini are located. The dark dots at loop termini represent multiprotein complexes, e.g., replication complexes.

3) Incompatibility between GC- and AT-rich blocks can contribute to micelle formation. On average, R minibands are at least ~3% richer in GC content than G minibands (Saitoh and Laemmli, 1994). Approximately 80% of the known genes are found in R minibands (Craig and Bickmore, 1993). This suggests that more histones are chemically modified, e.g., acetylated, in R than in G blocks, and more transcription complexes are bound to R than to G blocks. Because of the large size of the blocks (0.5-1 Mbp), even a small difference in interaction energy per monomer between unlike versus like monomers can lead to block separation: the blocks of one type form the loop termini, and the blocks of the other type are located at loop apices. Marko and Siggia (1997a) suggested that one can determine the parameters of the GC versus AT incompatibility by mixing bacterial DNA molecules that have very different GC contents.

Another contribution to stabilization of the chromosome micelle structure could come from multiprotein complexes, which participate in many chromosome functions at various stages of the cell cycle, e.g., transcription, replication, and chromosome condensation. Multiprotein complexes may associate differently with R and G blocks, as was suggested above for transcription complexes. For DNA replication, we assume that replication complexes are located in the micelle cores at the onset of S phase (see below).

4) The average loop size in interphase chromosomes in the model is assumed to be in the range of 1-2 Mbp. This is consistent with a loop containing two isochores or two replication minibands of size 0.5-1.0 Mbp, which seems to be reasonable (Bernardi, 1995; Simon and Cedar, 1996). The range of loop sizes for a large number of mammalian cell lines was estimated to be between 1.2 and 2.2 Mbp (Johnston et al., 1997). In the nuclei of early embryos of Drosophila, contacts between chromatin and the nuclear envelope have a frequency of one per 1-2 Mbp (Marshall et al., 1996). Because a loop contacts the nuclear envelope at its apex in the presented model, this leads to a 1- to 2-Mbp loop size.

Some Properties of Chromatin and Micelles

Individual loops in micelles behave as independent Gaussian coils. On the other hand, micelles are not interpenetrating (Semenov et al., 1995). Thus, the micelle structure of interphase chromosomes reconciles the contradiction mentioned in the INTRODUCTION between the discreetness of chromosomal territories and the random-walk behavior of chromatin.

Random-walk behavior of chromatin was demonstrated (van den Engh et al., 1992; Sachs et al., 1995; Yokota et al., 1995, 1997) by the linear dependence of the mean square of the geometrical distance between two probes on the same chromosome, < hx2> (µm2), versus their genomic distance, Mx (Mbp):
<UP>⟨h</UP><SUB><UP>x</UP></SUB><SUP><UP>2</UP></SUP><UP>⟩ = BM<SUB>x</SUB></UP> (1)
where the coefficient B (µm2/Mbp) describes chromatin compactness. Equation 1 is valid on the scale up to ~1 Mbp, no matter where two probes are located on the chromosome, and deviations from the linear dependence are observed for Mx > 1 Mbp. Equation 1 describes the behavior of ideal Gaussian chains, which particularly occurs in polymer melts or theta  solvents (de Gennes 1979; Grosberg and Khokhlov, 1994). Although the conditions for chromosomes in nuclei might not be the same as in a melt or a theta  solvent, here we consider ideal Gaussian chains as a first approximation, because the values of B have been obtained using this assumption (van den Engh et al., 1992; Sachs et al., 1995; Yokota et al., 1995, 1997).

It is known (e.g., see de Gennes, 1979) that < hx2> can also be represented as < hx2>  = bLx, where b is the length of the Kuhn statistical segment, and Lx is the fiber contour length. These two quantities are interrelated through k, the mass of the Kuhn statistical segment, and Mx: b/k = Lx/Mx. Thus, B in Eq. 1 can be expressed as B = b2/k.

Because the 30-nm chromatin fiber has ~0.2 kbp per nucleosome (van Holde, 1989; Wolffe, 1995), n, the number of nucleosomes per 10 nm of chromatin fiber contour length, can be estimated from the above expressions, as:
<UP>n = A</UP>(<UP>b/B</UP>) (2)
where the coefficient A = 50 µm/Mbp (=10 nm/0.2 kbp).

It has been shown that values of B for chromatin in R minibands are ~2.5-fold greater than those in G minibands, independent of fixation technique (Yokota et al., 1997). The fixation technique strongly affects the absolute values of BG and BR (values of B for G and R minibands, respectively), in parallel with nuclear size (Sachs et al., 1995; Yokota et al., 1995, 1997). For paraformaldehyde-fixed human fibroblast nuclei, nuclear size is not changed by fixation, and BG = 0.5 µm2/Mbp and BR = 1.3 µm2/Mbp (Yokota et al., 1997). We shall use these values and their average value B = 0.9 µm2/Mbp for the calculations below.

The larger value of BR relative to BG means that chromatin in R minibands is stretched in comparison with that in G minibands. In the accordion-like structure of the chromatin fiber (Woodcock et al., 1993; Horowitz et al., 1994; Woodcock and Horowitz, 1995), angles between links increase under stretching, which leads to an increase in the ratio of chromatin contour length to chromatin mass, L/M, which = b/k, and an increase in b/k leads to an increase in the values of B (see above).

Experimental data (Castro 1994; also see Marko and Siggia, 1997b) indicate that the Kuhn segment length b ~ 60 nm. Substituting BG = 0.5 µm2/Mbp and BR = 1.3 µm2/Mbp, and b = 60 nm in Eq. 2, one obtains n = 6 and 2.3 nucleosomes per 10-nm contour length for G and R minibands, respectively. These values are consistent with experimental data for chromatin structure: n = 6-8 nucleosomes per 10 nm for a compact chromatin fiber and n = 1-2 nucleosomes per 10 nm for a stretched chromatin fiber (van Holde and Zlatanova, 1995, 1996; Woodcock and Horowitz, 1995).

It has been shown for micelles that the polymer blocks that form micelle cores are stretched (see e.g., Semenov et al., 1995, 1996), because a large block incompatibility favors an increase in micelle size, and this leads to stretching of polymer blocks in the micelle cores. Applied to R and G minibands in chromatin micelles, this suggests that R blocks, which are stretched, are located at loop termini, and G blocks, which are unstretched, are located at loop apices. This assignment of R and G minibands is consistent with their replication time patterns (see below).

Dependence of Mean-Square Geometrical Distance on Genomic Distance for G1 Phase Chromosomes

The dependence of mean-square geometrical distance, < hx2> , on genomic distance, Mx, has been obtained on the 0.1-200 Mbp scale for three chromosomes (4, 5, and 19) in fixed human fibroblasts (Sachs et al., 1995; Yokota et al., 1995). These data can be summarized as having experimental points located between two parallel lines that have a shallow slope of ~20-fold less than the slope of this dependence over a short range (<1 Mbp). The authors suggested a model of chromosome structure that includes ~3-Mbp loops containing flexible chromatin that corresponds to a steep slope and a much less flexible nonchromatin backbone that corresponds to a shallow slope (also see Liu and Sachs, 1997). However, the measurements under separation of <1 Mbp, wherever one looks in the chromosome, never reveal a shallow slope (Yokota et al. 1995, 1997); this puts in doubt the existence of a rigid backbone.

The presented model suggests that intermicelle links and micelle tails contain the same material as micelle loops, the 30-nm chromatin fiber. The < hx2> versus Mx dependence following from this model is presented schematically in Figure 3. The net increase in < hx2> inside a micelle is zero, because the loop termini are located randomly and very close to each other in the micelle core. Thus, the < hx2> versus Mx dependence in the model is due to chromosome tails and intermicelle links.


View larger version (13K):
[in this window]
[in a new window]
 
Figure 3.   Schematic drawing of the dependence of < hx2> , mean-square geometrical distance, vs. Mx, genomic distance between two points on a G1 chromosome. The < hx2> vs. Mx dependence is due to chromosome tails and intermicelle links, because the net increase in < hx2> inside a micelle is zero. Dashed lines represent the boundaries of the experimental points. The model predicts that their slope (Bapp) is f-fold shallower than the initial slope (see Eq. 6), where f is the average number of loops per micelle.

Let us consider the < hx2> versus Mx dependence quantitatively. Suppose a G1 chromosome of size Mo (megabase pairs) forms several micelles with an average loop size, Mf (megabase pairs), and an average number of loops per micelle, f. The number of micelles per chromosome, m, can be expressed as:
<UP>m ≈ M<SUB>o</SUB>/fM<SUB>f</SUB></UP> (3)
As shown in Figure 3, straight lines for the boundaries of the experimental points can be defined as those connecting a telomere and the opposite tail base. Both lines have the same slope, Bapp, which can be defined as:
<UP>⟨H<SUP>2</SUP>⟩ = B</UP><SUB><UP>app</UP></SUB>(<UP>M<SUB>o</SUB> − M</UP><SUB><UP>f</UP></SUB>)<UP> ≈ B<SUB>app</SUB>M<SUB>o</SUB></UP> (4)
On the other hand, the same mean-square distance can be expressed as caused by a chromosome tail and the m-1 intermicelle links:
<UP>⟨H<SUP>2</SUP>⟩ = mBM<SUB>f</SUB> = BM<SUB>o</SUB>/f</UP> (5)
Equating Eqs. 4 and 5, one obtains the expression for the average number of loops per micelle, f:
<UP>f = B/B<SUB>app</SUB></UP> (6)
The meaning of Eq. 6 is that because there is one linear intermicelle link per f loop, representation of micelle structure as a linear coil yields a slope = Bapp, which is f-fold shallower than that for the micelle. The estimate of f from Eq. 6 is independent of any parameter of the model and equals the ratio of two measurable values, which was found to be ~20 (Sachs et al., 1995; Yokota et al., 1995). Thus, data for the geometrical versus genomic distances for human G1 fibroblasts (Sachs et al., 1995; Yokota et al., 1995) suggest that the average number of loops per micelle, f, is ~20 in these cells.

The number of loops per micelle in G1 fibroblast chromosomes (f ~ 20) is comparable with the number of loops per micelle for ionomers, which is 5-50 (Semenov et al., 1995, 1996). The value of f is limited by the maximal number of polymer chains that can be brought together in a micelle core (Semenov et al., 1996).

A crude estimate of flim, the maximal number of loops per micelle, is as follows. The number of loop termini confined in the micelle core is approx f, the number of loops per micelle (exactly f + 1). Suppose Dc is the micelle core diameter, L is the loop terminus contour length, and h is the distance between the entrance and exit points of a loop terminus. The total volume occupied by chromatin fibers in the micelle core is approx fd2L, where d is the chromatin fiber diameter (=30 nm). The average contour length, L, can be expressed as < h2> /b (see above), where < h2> , the mean-square average cord length in a sphere, = Dc2/2, and b is the Kuhn segment length (=60 nm). Equating the total volume occupied by loop termini to the micelle core volume = pi Dc3/6, one can estimate flim as:
<UP>f<SUB>lim</SUB> ≈ </UP>(<UP>&pgr;/3</UP>)<UP>bD<SUB>c</SUB>/d<SUP>2</SUP></UP> (7)
It is argued below that chromatin replication starts at the micelle cores. Taking the diameter of replication "factories" (~0.2-0.3 µm [Hozak et al., 1993; Tomilin et al., 1995]) as a range of Dc, Eq. 7 yields flim ~ 10-20; this is comparable with the value of f estimated above.

The number of micelles per G1 nucleus, Nm, can be estimated as:
<UP>N<SUB>m</SUB> = C/fM<SUB>f</SUB></UP> (8)
where C is the G1 DNA content. Diploid human cells have C = 6.4 Gbp, and for f = 20 and Mf = 1-2 Mbp, Eq. 8 yields Nm = 160-320. This estimate is important for discussion of the relationship between micelle structure and DNA replication.

Micelle Structure and DNA Replication

In high eukaryotes, but not in yeast, chromosome structure plays an important role in replication initiation (Coverley and Laskey, 1994; Laskey and Madine, 1996; Gilbert, 1998). Let us show that the proposed model is consistent with data for early S phase replication in mammalian cells if we suggest that at the onset of S phase, chromatin replication is initiated at the loop termini in micelle cores. The following points support this suggestion.

1) If replication starts in the micelle cores, R blocks that form loop termini (see above) should replicate earlier than G blocks. This is consistent with the fact that R minibands replicate earlier than G minibands (for review see Drouin et al., 1994). Thus, in this model two features of R minibands, stretched chromatin and early replication, could be explained by their location in micelle cores.

2) In the model, chromatin fibers contact the nuclear membrane at loop apices. Thus, the observation that the nuclear periphery contains predominantly late-replicated G minibands (Ferreira et al., 1997) is consistent with the assignment of G minibands to loop apices.

3) If replication starts in the micelle cores at the beginning of S phase: a) the number of replication clusters per nucleus should be similar to the number of micelle cores (Nm) in G1 phase; b) the number of minifoci per cluster should be similar to the number of loop termini per core (f); and thus, c) the total number of minifoci per nucleus should be similar to the product fNm. At the onset of S phase, nuclei with C ~ 6 Gbp (e.g., diploid human fibroblasts and V79 cells) have 100-300 replication clusters with ~20 minifoci per cluster (Berezney et al., 1995a; Jackson and Cook, 1995); this is consistent with the estimates obtained above: Nm = 160-320 and f ~ 20. Cell lines with C = 9-10 Gbp (e.g., mouse 3T3 and human HeLa) have 600-750 replication clusters with ~10-12 minifoci per cluster (Jackson and Pombo, 1998), or a total number ~6000-9000 minifoci per nucleus, which is consistent with fNm = C/Mf = 4500-10000 for these cells. These data suggest that f is ~10-20 for various cells.

4) Replication-labeled clusters were observed through several cell cycles (Jackson and Pombo, 1998), and their number doubles in G2 versus G1 (Jackson and Pombo, 1998; Zink et al., 1998). The size of these labeled chromatid subdomains is ~0.4-0.8 µm (Zink et al., 1998). This is comparable with our estimate of the size range of G and R minibands, i.e., 0.5-1.1 µm (=[BMf/2]1/2, see Eq. 1) for B = 0.5-1.3 µm2/Mbp (Yokota et al., 1997) and Mf = 1-2 Mbp.

Thus, the suggestion that replication clusters can be considered a fundamental aspect of the higher-order structure of the genome (Berezney et al., 1995b; Jackson and Pombo, 1998; Zink et al., 1998) can have micelle cores as its basis.

Relationship Between Size of Interphase Nucleus and Chromosome Compactness

One test of the model is that the chromosome size estimated with the model must not exceed the size (length and thickness) of the corresponding nucleus, and that the total chromosomal volume or area must not exceed the nuclear volume or area. To estimate the micelle diameter, a micelle can be represented as a star-branched polymer with branch size of Mf/2. Because chromatin in loops behaves as a Gaussian chain, and because the branch ends are close to each other in the micelle core, the mean-square micelle diameter, < Dm2> , is equal to double the mean-square branch size:
<UP>⟨D</UP><SUB><UP>m</UP></SUB><SUP><UP>2</UP></SUP><UP>⟩ = BM<SUB>f</SUB></UP> (9)
where the average value B = (BG + BR)/2.

For human fibroblasts, replacing Mf = 1-2 Mbp and B = 0.9 µm2/Mbp in Eq. 9, one obtains < Dm2> 1/2 = 0.9-1.3 µm, which is close to the smallest nuclear thickness measured for cultured cells (=1.2 µm) observed for human AG1522 fibroblasts in monolayer (Raju et al., 1991). This suggests a two-dimensional (2-D) organization of micelles in a monolayer of flattened AG1522 cells.

The mean-square chromosome length, < H2> , can be presented as the end-to-end distance for a random walk of m micelles, each of length (diameter) Dm. Replacing Dm and m from Eqs. 3 and 9, one obtains:
<UP>⟨H<SUP>2</SUP>⟩ = D</UP><SUB><UP>m</UP></SUB><SUP><UP>2</UP></SUP><UP>m = BM<SUB>o</SUB>/f</UP> (10)
Equation 10 is the same as Eq. 5, where the expression for < H2> was obtained as a random walk of the intermicelle links. Thus, two approaches yield the same values of H. For human fibroblasts, replacing f = 20 and B = 0.9 µm2/Mbp in Eq. 10, one obtains < H2> 1/2 ~ 4 µm for the largest (Mo = 263 Mbp) human chromosome. For a fibroblast in monolayer, this is much smaller than its nuclear dimensions (~10-30 µm [Yokota et al., 1997]).

The 2-D chromosome territories have been measured for chromosomes 17 (=4.1 µm2; Mo = 92 Mbp) and the inactive X (=5.2 µm2; Mo = 164 Mbp), respectively, in human fibroblasts (Clemson et al., 1996). For a chromosome consisting of m micelles, the 2-D chromosome area, Ac, can be presented as a sum of micelle areas, = (pi /4)Dm2m, which is only slightly different from the expression for < H2> (see Eq. 10). It follows from Eqs. 3 and 9 that Ac = (pi /4)BMo/f. Thus, B can be estimated from Ac as:
<UP>B = </UP>(<UP>4/&pgr;</UP>)<UP>A<SUB>c</SUB>f/M<SUB>o</SUB></UP> (11)
For f = 20 and the above values of Ac, Eq. 11 yields B = 1.1 and 0.8 µm2/Mbp for chromosomes 17 and X, respectively. These values are close to the average value of B = 0.9 µm2/Mbp (Yokota et al., 1997).

Chromosome compactness seems to be adjusted to nuclear size (Yokota et al., 1995, 1997; Sanchez et al., 1997). We shall estimate B from nuclear cross-section areas and nuclear volumes. The test is whether the estimated values of B fall within the range of experimental values. We assume that the chromosomes fill all available space in the nucleus. This assumption is consistent with observations of constrained diffusional motion of chromosomes in the nucleus (Abney et al., 1997; Marshall et al., 1997).

For 2-D nuclei, the available nuclear area is alpha An, where alpha  is the occupancy factor, and An is the nuclear cross-section area. Because there is a large number (hundreds) of micelles in a nucleus (see above), the occupancy factor alpha  for tightly packed micelles in a nucleus can be taken to be alpha  ~ pi /4, as for squared circles. Equating the available nuclear area to the sum of micelle areas, = Nmpi Dm2/4, and using Eqs. 8 and 9, one obtains an estimate for B from An:
<UP>B = A<SUB>n</SUB>f/C</UP> (12)
Human HSF7 fibroblasts (C ~ 6 Gbp) and HeLa cells (C ~ 9 Gbp) form a monolayer of flattened cells, and their nuclei may be approximated as 2-D. Experimental values of their nuclear cross-section areas are An = 160, 240, and 400 µm2 for the paraformaldehyde-fixed HSF7 and MAA-fixed HeLa and HSF7 cells, respectively (Yokota et al., 1997). For f = 20, Eq. 12 yields B = 0.5-1.3 µm2/Mbp, which is close to the experimental values (Yokota et al., 1997).

For 3-D nuclei, the available nuclear volume is alpha Vn, where Vn is the nuclear volume and alpha  is the occupancy factor, which can be taken to be ~pi /6, as for cubed spheres. Equating the available nuclear volume to the sum of micelle volumes, Nmpi < Dm2> 3/2/6, and using Eqs. 8 and 9, one obtains:
<UP>B = </UP>(<UP>6&agr;/&pgr;</UP>)<SUP><UP>2/3</UP></SUP>(<UP>V<SUB>n</SUB>/C</UP>)<SUP><UP>2/3</UP></SUP> (<UP>f<SUP>2/3</SUP>/M</UP><SUB><UP>f</UP></SUB><SUP><UP>1/3</UP></SUP>) (13)
Table 1 presents estimates of B from Eq. 13 for f = 10-20, Mf = 1-2 Mbp, and alpha  = pi /6. It follows from Table 1 that the whole range of estimated B, 0.2-1.5 µm2/Mbp, is consistent with the range of experimentally determined values (Yokota et al., 1997). The values of n, the number of nucleosomes per 10 nm of chromatin fiber contour length estimated from Eq. 2, range between 2 and 13, which seems to be reasonable (van Holde and Zlatanova, 1995, 1996; Woodcock et al., 1995).

                              
View this table:
[in this window]
[in a new window]
 
Table 1.  Relationship between nuclear volume and coefficient B

Equation 13 is not applicable to cells with very small chromosomes, such as the budding yeast Saccharomyces cerevisiae, because all or most of their chromosomes are smaller than Mf = 1-2 Mbp. We suggest that yeast chromosomes have a linear coil configuration. Because replication initiation in higher eukaryotes is suggested to begin in the micelle cores (see above), the lack of micelle structure in yeast chromosomes is consistent with the observation that higher and lower eukaryotes have different patterns of replication initiation (see review in Gilbert, 1998).

The expression for B for linear coil polymers (e.g., yeast chromosomes) can be obtained from Eq. 13 by replacing f with 1 and Mf with the average chromosome size Mo = C/N, where N is the number of chromosomes per nucleus (N = 32 for S. cerevisiae):
<UP>B = </UP>(<UP>6&agr;/&pgr;</UP>)<SUP><UP>2/3</UP></SUP>(<UP>V<SUB>n</SUB>/C</UP>)<SUP><UP>2/3</UP></SUP> (<UP>N/C</UP>)<SUP><UP>1/3</UP></SUP> (14)
Equation 14 yields B = 0.6 µm2/Mbp for S. cerevisiae, which is close to the estimated values of B obtained for the micelle structure of high eukaryote chromosomes (see Table 1). This is consistent with the observation that chromatin in interphase yeast cells has the same relationship between geometrical versus genomic distances as that in mammalian cells on the 1-Mbp scale (Guacci et al., 1994).

    CONCLUSIONS
Top
Abstract
Introduction
Results & Discussion
Conclusions
References

1) A G1 phase mammalian chromosome can be approximated as a multiblock copolymer containing two alternating types (R and G) of polymer blocks, which form a string of loop clusters (micelles), with each loop ~1-2 Mbp in size. Application of the model to the experimental data (Sachs et al., 1995; Yokota et al., 1995) for the dependence of geometrical versus genomic distances between two points on the same chromosome yields an estimate of ~20 loops per micelle.

2) The number of micelles per nucleus is close to the observed number of replication clusters at the onset of S phase, and the number of loops per micelle is close to the number of replication minisites per cluster. This is consistent with loop termini being sites of initiation of DNA replication at the onset of S phase.

3) R minibands form loop termini, whereas G minibands are located at loop apices. This conclusion follows from relating the chromatin fiber being stretched in R minibands (Yokota et al., 1997) to a known feature of micelles, that polymer blocks located in micelle cores are stretched. These locations of R and G minibands are consistent with their replication pattern; the former are replicated earlier than the latter.

4) The chromosome micelle structure describes the relationship between chromosomal and nuclear sizes for several types of higher-order eukaryotic cells (insects, plants, and mammals). For yeast cells, this relationship is described by a linear coil configuration of chromosomes.

    ACKNOWLEDGMENTS

I am grateful to C.S. Lange for very helpful discussions and careful reading of the manuscript, A. Berens for drawing the figures, and the anonymous reviewers for their constructive suggestions.

    FOOTNOTES

Author's E-mail address: ostasj23{at}hscbklyn.edu.

    REFERENCES
Top
Abstract
Introduction
Results & Discussion
Conclusions
References


Molecular Biology of the Cell
Vol. 9, 3031-3040, November 1998
Copyright © 1998 by The American Society for Cell Biology



This article has been cited by other articles:


Home page
Mol. Biol. CellHome page
J. Ostashevsky
A Polymer Model for Large-scale Chromatin Organization in Lower Eukaryotes
Mol. Biol. Cell, June 1, 2002; 13(6): 2157 - 2169.
[Abstract] [Full Text] [PDF]


Home page
Mol. Biol. CellHome page
C. Carvalho, H. M. Pereira, J. Ferreira, C. Pina, D. Mendonca, A. C. Rosa, and M. Carmo-Fonseca
Chromosomal G-dark Bands Determine the Spatial Organization of Centromeric Heterochromatin in the Nucleus
Mol. Biol. Cell, November 1, 2001; 12(11): 3563 - 3572.
[Abstract] [Full Text] [PDF]


Home page
Genome Res.Home page
A. Nekrutenko and W.-H. Li
Assessment of Compositional Heterogeneity Within and Between Eukaryotic Genomes
Genome Res., December 1, 2000; 10(12): 1986 - 1995.
[Abstract] [Full Text]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ostashevsky, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ostashevsky, J.


Home Help [Feedback] [For Subscribers] [Archive] [Search] [Contents]