Molecular and Cellular Biology, November 2001, p. 7682-7695, Vol. 21, No. 22
0270-7306/01/$04.00+0 DOI: 10.1128/MCB.21.22.7682-7695.2001
Copyright © 2001, American Society for Microbiology. All rights reserved.
Department of Biochemistry and Molecular Biology,1 Center for Mammalian Genetics,2 and Division of Pediatric Genetics,3 University of Florida, Gainesville, Florida 32610
Received 8 June 2001/Returned for modification 29 July 2001/Accepted 20 August 2001
| |
ABSTRACT |
|---|
|
|
|---|
Differential chromatin structure is one of the hallmarks distinguishing active and inactive genes. For the X-linked human hypoxanthine phosphoribosyltransferase gene (HPRT), this difference in chromatin structure is evident in the differential general DNase I sensitivity and hypersensitivity of the promoter regions on active versus inactive X chromosomes. Here we characterize the nucleosomal organization responsible for the differential chromatin structure of the active and inactive HPRT promoters. The micrococcal nuclease digestion pattern of chromatin from the active allele in permeabilized cells reveals an ordered array of translationally positioned nucleosomes in the promoter region except over a 350-bp region that is either nucleosome free or contains structurally altered nucleosomes. This 350-bp region includes the entire minimal promoter and all of the multiple transcription initiation sites of the HPRT gene. It also encompasses all of the transcription factor binding sites identified by either dimethyl sulfate or DNase I in vivo footprinting of the active allele. In contrast, analysis of the inactive HPRT promoter reveals no hypersensitivity to either DNase I or a micrococcal nuclease and no translational positioning of nucleosomes. Although nucleosomes on the inactive promoter are not translationally positioned, high-resolution DNase I cleavage analysis of permeabilized cells indicates that nucleosomes are rotationally positioned over a region of at least 210 bp on the inactive promoter, which coincides with the 350-bp nuclease-hypersensitive region on the active allele, including the entire minimal promoter. This rotational positioning of nucleosomes is not observed on the active promoter. These results suggest a model in which the silencing of the HPRT promoter during X chromosome inactivation involves remodeling a transcriptionally competent, translationally positioned nucleosomal array into a transcriptionally repressed architecture consisting of rotationally but not translationally positioned nucleosomal arrays.
| |
INTRODUCTION |
|---|
|
|
|---|
Differential chromatin structure and accessibility, particularly at the promoter, have long been recognized as characteristics that distinguish active from inactive genes. Active genes are in general more accessible to regulatory factors than inactive genes, as indicated by nuclease sensitivity. In addition, the promoters of active genes often exhibit marked DNase I hypersensitivity, especially in the vicinity of transcription factor binding sites (10, 20). This hypersensitivity is postulated to be due to changes in the chromatin architecture of the promoter and may represent nucleosomal remodeling or displacement, stretches of single-stranded DNA, torsionally stressed DNA, or other distortions in chromatin structure arising from factor binding (18, 20, 62).
The functional effect of nucleosomes on transcription initiation is thought to be repressive since in vitro assembly of nucleosomal arrays on DNA templates drastically reduces the capacity of these templates to support basal transcription (25, 35, 36, 57). Furthermore, the differential accessibility and transcriptional potential of chromatin structure in active versus inactive promoters are often associated with differential nucleosomal organization (3, 5, 23, 30, 46, 62). Thus, remodeling the nucleosomal architecture of a promoter is likely to be an integral feature of mechanisms of gene activation and/or silencing, which may involve histone acetylation and chromatin-remodeling complexes such as SWI/SNF (15, 20, 60).
The organization of genomic DNA into nucleosomal arrays is defined by both the translational position of the nucleosome relative to the linear nucleotide sequence and the rotational orientation of the DNA helix relative to the surface of the histone octamer. The translational position of nucleosomes on a DNA template (i.e., the linear position of the nucleosome relative to the DNA sequence [46]) has been shown to affect the accessibility of cis-acting elements in the promoter to various transcription factors as well as the basal transcription complex. Whether a transcription factor binding site is incorporated into a nucleosomal core or is in the linker region between nucleosomes can dramatically affect its accessibility to its cognate binding protein(s) in vitro (25, 58). However, the extent to which incorporating a transcription factor binding site into a nucleosomal core reduces the accessibility of that site can vary widely among transcription factors (4, 25). The rotational orientation of the DNA helix wound around a nucleosomal core also affects the accessibility of that DNA to transcription factors (25). For instance, the rotational orientation of the TATA box, the glucocorticoid response element, and the thyroid response element within a nucleosome strongly affects the accessibility of these elements to their cognate binding factors (14, 24, 59). These findings suggest that both the translational position and rotational orientation of cis-acting regulatory elements relative to those of the histone octamer can significantly affect their accessibility to regulatory factors and therefore their function in the context of nucleosomal arrays.
In addition to their repressive effects on transcription, nucleosomal positioning and orientation appear, in some cases, to modulate transcriptional activation. In particular, translational positioning of nucleosomes can mediate the cooperative binding of transcription factors such as Gal4 (48, 54) and pho4 (53), allowing rapid transcriptional induction. Furthermore, the nucleosomal organization of chromatin may establish a sequential hierarchy of transcription factor binding, with those factors that can bind to their cognate sites in nucleosomes binding first, leading to remodeling of adjacent nucleosomes so other factors can subsequently bind and activate transcription (21, 35, 48). In addition, precisely positioned nucleosomes can juxtapose linearly distant transcription factor binding sites, allowing the factors bound to these sites to interact (29, 45, 46). Alternatively, it has been suggested that the curvature of DNA along the surface of the histone octamer may reduce steric hindrance, which would otherwise occur between adjacent factors bound to closely spaced binding sites on linear DNA (30, 52). Together, these studies strongly suggest that examination of the nucleosomal organization of endogenous promoters is vital to the understanding of their regulation.
Here, we examine the nucleosomal structure of the promoter of the X-linked human hypoxanthine phosphoribosyltransferase gene (HPRT) on the active and inactive X chromosomes. The HPRT gene is subject to X chromosome inactivation, a process that leads to the transcriptional silencing of genes on one of the two X chromosomes in each female somatic cell (9). This results in the presence of both a transcriptionally active and a transcriptionally inactive HPRT allele within each female nucleus. The HPRT promoter lies within a CpG island, lacks a TATA box, and contains a potential AP-2 binding site, a cluster of five GC boxes, a potential initiator element, and a region of multiple transcription initiation sites (13, 16, 32, 42).
Several epigenetic characteristics distinguish the active and inactive HPRT alleles, including differential DNA methylation, general DNase I sensitivity, and DNase I hypersensitivity. On the active allele, the promoter region in vivo is unmethylated (11), is relatively sensitive to DNase I (26), and contains a DNase I-hypersensitive site that maps to the 5' flanking region (13, 26, 56). Multiple transcription factor binding sites, including the potential AP-2 binding site, all five GC boxes, and the potential initiator element, have been identified in the active promoter by dimethyl sulfate (DMS) in vivo footprinting (13). In contrast, the inactive promoter is densely methylated (11), is resistant to DNase I, and does not exhibit DNase I hypersensitivity (26) or detectable transcription factor binding in vivo (13).
To examine the nucleosomal organization of the HPRT promoter, we investigated both the translational and rotational positioning of nucleosomes in the active and inactive HPRT promoter regions in permeabilized cells. Micrococcal nuclease (MNase) analysis showed that the active promoter was assembled into an ordered array of translationally positioned nucleosomes which was interrupted over a 350-bp region that showed increased accessibility to MNase and DNase I and that contained the entire functional promoter. In contrast, the inactive promoter was relatively inaccessible to both nucleases and transcription factors and was not assembled into a translationally positioned nucleosomal array. However, high-resolution DNase I cleavage analysis revealed rotational positioning of nucleosomes over a 210-bp region on the inactive promoter coincident with the 350-bp nuclease-accessible region on the active promoter. This differential organization of nucleosomes between the active and inactive HPRT promoters suggests that the transcriptional silencing of the HPRT promoter by X inactivation involves reconfiguration of the translationally positioned nucleosomal array in the active promoter to rotationally, but not translationally positioned nucleosomal arrays on the inactive promoter. Furthermore, it suggests that this rotational positioning of nucleosomes may be involved in the transcriptional repression of the HPRT promoter.
| |
MATERIALS AND METHODS |
|---|
|
|
|---|
Cell lines. Two cell lines, 8121 and 4.12, described previously (13), were used for analysis of the active and inactive human HPRT promoter regions. 8121 is a human-hamster hybrid containing an inactive human X chromosome, whereas 4.12 is a human-hamster hybrid containing an active human X chromosome. 4.12 cells were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum, 1% penicillin-streptomycin, and 1× hypoxanthine-aminopterin-thymidine (Gibco/BRL) supplement. 8121 cells were maintained in DMEM supplemented with 10% fetal bovine serum, 1% penicillin-streptomycin, and 1× 6-thioguanine (2-amino-6-mercaptopurine; final concentration, 5 µg/ml; Sigma). Cells were maintained in culture at 37°C in 5% CO2.
MNase treatment of permeabilized cells MNase treatment of permeabilized cells was based on a protocol for DNase I treatment of permeabilized cells modified from that described by Rincon Limas et al. (41). Cells were grown to confluence in T-150 culture flasks, trypsinized, combined in a 50-ml conical tube, and gently pelleted at ~500 × g in a tabletop centrifuge. They were then washed gently once in solution A (150 mM sucrose, 80 mM KCl, 35 mM HEPES [pH 7.4], 5 mM K2HPO4, 5 mM MgCl2, 0.5 mM CaCl2) and once in solution B (150 mM sucrose, 80 mM KCl, 35 mM HEPES [pH 7.4], 5 mM K2HPO4, 5 mM MgCl2, 2 mM CaCl2) and were then resuspended at 500 µl per T-150 flask in solution B. For each MNase treatment, 500 µl of cells was combined with 500 µl of solution B plus 0.4% NP-40 (Sigma) and the desired amount (0 to 100 U [see Fig. 2B] or 0 to 200 U [see Fig. 2A]) of MNase (resuspended at 20 U/µl; Sigma), gently mixed, and then incubated at room temperature for 2 min. The digestion was stopped by adding 4 ml of DNA lysis buffer (50 mM Tris [pH 8.0], 150 mM NaCl, 25 mM EDTA [pH 8.0], 0.5% sodium dodecyl sulfate [SDS], 300 µg of proteinase K/ml), and the lysate was then incubated overnight at room temperature. MNase-treated DNA was then isolated by standard phenol extraction and ethanol precipitation techniques (43), resuspended at 1 µg/µl in TE (10 mM Tris, 1 mM EDTA [pH 8.0]), and stored at 4°C.
DNase I treatment of permeabilized cells for mapping DNase I-hypersensitive sites. Cells were grown to confluence in a T-75 culture flask and then washed with phosphate-buffered saline and trypsinized. They were pelleted by centrifugation at 500 × g in a tabletop centrifuge and then washed once with solution A and resuspended in 500 µl of solution B. Then 500 µl of solution B containing 0.4% NP-40 and 20 or 40 µg of DNase I (1 mg/ml; Worthington) was added to the cell suspension, and the suspension was mixed and incubated for 2 min at 37°C. The digestion was stopped with 4 ml of lysis buffer (see above), and the suspension was incubated overnight at room temperature. DNA was subsequently isolated by standard techniques (43), resuspended at 1 µg/µl in TE, and stored at 4°C.
MNase treatment of naked DNA. Ten micrograms of purified genomic DNA was digested to completion with BclI for each MNase concentration point. The BclI-digested DNA was then extracted with phenol-chloroform, precipitated with ethanol, and resuspended in 100 µl of Ex50 buffer (10 mM HEPES [pH 7.6], 60 mM KCl, 1.5 mM MgCl2, 0.5 mM EGTA [pH 8.0], 10% glycerol, 10 mM glycerol phosphate) (1) plus 5 mM CaCl2. Then 0.125 to 1 U (see Fig. 2B) or 1 to 8 U (see Fig. 2A) of MNase was added, and the reaction mixture was allowed to incubate for 2 min at room temperature. The reaction was stopped by adding 50 µl of stop solution (2.5% Sarkosyl, 100 mM EDTA [pH 8.0]). The treated DNA was extracted with phenol-chloroform, precipitated with ethanol, resuspended at 1 µg/µl in TE, and stored at 4°C.
Southern transfer, probe preparation, and hybridization.
MNase-treated DNA samples from permeabilized cells and naked DNA were
digested to completion with BclI, size-fractionated on a
1.4% agarose gel at 60 V for 12 h, stained with ethidium bromide,
and transferred to Hybond N+ by capillary transfer. Indirect end
labeling was performed using a 400-bp hybridization probe located just
upstream of a BclI restriction site in intron I of HPRT (see Fig. 1). This probe was amplified by PCR from a
plasmid containing a BamHI/PstI restriction
fragment from intron 1 that includes a previously described single-copy
region (26) using the following primer set: MnaseBcl400
(GTTTGGGGTGCGATGGTGAGG) and MnaseBcldownstream
(CAGAACGGTTGAGGAGGGAGGCCA). The PCR product was gel purified
using the Qiagen gel extraction kit and radiolabeled by random priming.
Hybridization was performed as described by Hornstra and Yang
(12) at 65°C overnight in 4 ml of hybridization solution
(0.25 M Na2HPO4 [pH 7.2]
with o-phosphoric acid, 7% SDS, 1% bovine serum albumin
[fraction V], 1 mM EDTA [pH 8.0]). The blot was then washed at
65°C in wash solution (20 mM
Na2HPO4 [pH 7.2; with
o-phosphoric acid], 1% SDS, 1 mM EDTA) and exposed to Kodak MR film for 4 to 5 days with intensifying screens at
80°C.
Reconstitution of nucleosomes onto the HPRT promoter in vitro. The DNA template used for in vitro reconstitution of chromatin was the pBS HPRT 1.8-kb plasmid, which is a pBluescript-derived plasmid containing a 1.8-kb EcoRI-to-BamHI fragment that includes the entire HPRT promoter. DNA methylation at CpG dinucleotides in vitro was performed using HpaII, HhaI, and SssI methylases (New England Biolabs) essentially as described by the supplier. The methylated DNA was then extracted with phenol-chloroform, precipitated, and resuspended at 100 ng/µl in TE. Completeness of methylation was assayed by digestion of in vitro-methylated templates with HpaII and HhaI. Chromatin reconstitution was performed with a Drosophila melanogaster chromatin assembly extract essentially as described by Becker et al. (1). The reconstituted chromatin (total volume = 70 µl) was mixed with 100 µl of a preassembled MNase mixture (94 µl of Ex50 buffer, 5 µl of 100 mM CaCl2, 1 µl of MNase [50 U/µl; Sigma]) and incubated at room temperature. At 2 and 5 min, 40 µl of this digestion mixture was removed and mixed with 20 µl of stop solution (2.5% Sarkosyl, 100 mM EDTA [pH 8.0]) to stop the reaction. To remove RNA and proteins, 1 µl of RNase cocktail (Ambion), 8 µl of 2% SDS, and 5 µl of proteinase K (resuspended at 10 mg/ml; Gibco/BRL) were added at each time point and the mixtures were incubated at 37°C overnight. The digested DNA (approximately 150 ng per time point) was extracted with phenol-chloroform, precipitated, and resuspended in 10 µl of TE (pH 8.0). The purified DNA at each MNase digestion time point was digested to completion with BamHI. The digested DNA was size-fractionated on a 1.6% agarose gel at 100 V for 4 h and then transferred to Hybond N+ by capillary transfer. The blot was hybridized overnight at 40°C with radiolabeled oligonucleotide BamHINuc1Probe (5'-GCCCTGAGGCGCGGGATC-3'), which corresponds to the sequence immediately upstream of the BamHI site in the first intron of the HPRT promoter, and visualized by exposure to Kodak AR film at room temperature for 1 h.
DNase I treatment of permeabilized cells for high-resolution DNase I cleavage analysis. Cells were grown to 80% confluence in T-75 culture flasks. They were then washed once with phosphate-buffered saline, once with solution A, and once with solution B. One milliliter of solution B containing 0.2% NP-40 and 0 to 100 µg of DNase I (Worthington) was gently distributed over the monolayer of cells in each T-75 flask, and the cells were incubated at 37°C for 2 min. The DNase I digestion was then stopped by adding 4 ml of lysis buffer (see above), and the mixture was incubated overnight at room temperature. DNA was then isolated by standard phenol extraction and ethanol precipitation techniques (43), resuspended at 1 µg/µl in TE, and stored at 4°C.
DNase I treatment of naked DNA. Fifty micrograms of genomic DNA was first digested to completion with EcoRI for each DNase I concentration point. The EcoRI-digested DNA was then extracted with phenol-chloroform, precipitated with ethanol, and then resuspended in a solution containing 100 µl of distilled water and 200 µl of solution B. DNase I (1 mg/ml; Worthington) was added to a final concentration of 0.0125 to 0.05 µg/ml, and the DNA was digested for 2 min at room temperature. The reaction was then stopped by adding 12 µl of 0.5 M EDTA (pH 8.0) and 3 µl of 20% SDS. The DNA was extracted with phenol-chloroform, precipitated with ethanol, resuspended in TE, and stored at 4°C.
LMPCR. Amplification of DNase I-treated DNA was accomplished by ligation-mediated PCR (LMPCR) essentially as described by Hornstra and Yang (12), with an extension product capture modification described by Tormanen et al. (49) to reduce background. Primer extension was performed using Vent DNA polymerase (New England Biolabs), and PCR amplification was performed using Taq DNA polymerase (Gibco/BRL). The E1/E2 LMPCR primer set (13) was used to examine the upper strand in the minimal promoter, whereas the CA1/CA2 primer set (CCTAGTGAGCCTGCAAACTG/AAACTGGTAGGCGCCGGCGTAGG) was used to examine the lower strand. Additional upper- and lower-strand LMPCR primer sets C1/C2 and A1/A2, respectively (13), were used to analyze the region immediately upstream of the minimal promoter. Primer extension was performed essentially as described by Hornstra and Yang (12) using the same annealing temperature (45°C) and ramping parameters for all primer sets. For the PCR amplification, an annealing temperature of 64°C was used for all primers sets. Visualization of the footprint was achieved by size-fractionation of the PCR products on a DNA sequencing gel, followed by electrotransfer, hybridization, and autoradiography, essentially as described by Hornstra and Yang (12) using a radiolabeled single-strand-specific probe synthesized from single-stranded M13 templates containing the upper or lower strand of the HPRT promoter, as appropriate.
| |
RESULTS |
|---|
|
|
|---|
To determine the nucleosomal organization of the HPRT promoter region on the active and inactive X chromosomes and its relationship to transcription factor binding sites and nuclease hypersensitivity, the translational and rotational positioning of nucleosomes and the pattern of transcription factor binding in the HPRT promoter were examined. To examine translational positioning in vivo, NP-40-permeabilized cells were treated with MNase and then the MNase cleavage sites in the active and inactive HPRT promoters were mapped by Southern blotting and indirect end labeling. To identify regions of rotational nucleosomal positioning in the active and inactive promoters in vivo, NP-40-permeabilized cells were treated with DNase I and the high-resolution DNase I cleavage pattern of the promoter was determined by LMPCR. These high-resolution DNase I cleavage patterns were also used to delineate regions within the promoter that are occupied and footprinted in vivo by sequence-specific DNA-binding factors on the active and inactive HPRT alleles.
Translational positioning of nucleosomes in the HPRT
promoter region in vivo.
To examine the translational positioning
of nucleosomes in the human HPRT promoter region, 4.12 and
8121 cells (human-hamster hybrid cell lines containing either the
active or the inactive human X chromosome, respectively) were
permeabilized using detergent NP-40 and then treated with increasing
concentrations of MNase. MNase preferentially cleaves DNA within the
linker region between nucleosomal cores and, in conjunction with
Southern blotting and indirect end labeling, can detect the presence
and positions of translationally positioned nucleosomes within a given
region of chromatin. The positions of the MNase cleavages within
chromatin of the HPRT promoter region of permeabilized cells
relative to a downstream BclI site in the first intron of
the HPRT gene were mapped using a 400-bp hybridization probe
located just upstream of the BclI site (Fig.
1). To determine the positions of DNase I-hypersensitive sites relative to MNase cleavage sites,
NP-40-permeabilized cells containing the active HPRT allele
were treated with increasing concentrations of DNase I and the DNase
I-hypersensitive sites in chromatin of the HPRT promoter
relative to the same BclI site were also mapped by indirect
end labeling using the same hybridization probe.
|
|
260,
460,
660,
840,
1040,
1220, and
1420 relative to the translation initiation site (Fig. 2A, lanes 14 to 18, and B, lanes 13 to 16). This ~200-bp periodicity is consistent with the potential unit length of DNA in a mammalian nucleosome, suggesting that six
distinct translationally positioned nucleosomes are positioned immediately upstream of the potential AP-2 site (positions
264 to
272) on the active HPRT promoter (Fig. 1 and schematic
diagrams in Fig. 2A and B). Two additional translationally positioned
nucleosomes are suggested by the presence of a second ~200-bp
cleavage periodicity consisting of three weak MNase cleavage sites at
positions +90, +280, and +490 in the first intron of the
HPRT gene (Fig. 2A, lanes 14 to 18, and B, lanes 13 to 16).
While these downstream cleavage sites were relatively weak, they were
consistently reproducible in multiple independent MNase assays. This
pattern of periodic 200-bp cleavages observed in the chromatin of the
active allele was not observed on MNase-treated naked DNA through a
range of MNase concentrations that eventually achieve complete
digestion of the parental 4.3-kb BclI band (Fig. 2A, lanes
10 to 13, and B, lanes 9 to 12). This difference between the MNase
cleavage patterns of permeabilized 4.12 cells (Fig. 2A, lanes 14 to 18, and B, lanes 13 to 16) and naked DNA (Fig. 2A, lanes 10 to 13, and B,
lanes 9 to 12) suggests that the periodic cleavages observed in
permeabilized cells are due to nucleosomal organization rather than the
underlying DNA sequence.
An additional weak MNase cleavage site was observed in permeabilized
cells within the first downstream nucleosome at position +170 in the
active promoter (Fig. 2A, lanes 14 to 18). Since the MNase cleavages at
+90 and +280 suggest that a nucleosome is in fact positioned within
this region, the cleavage site at +170 may simply indicate an
intranucleosomal MNase cleavage site similar to those previously seen
at high resolution in the pS2 promoter by Sewack and Hansen
(46). Alternatively, the cleavage site at +170 could
represent (i) a modified nucleosome that is inherently more susceptible
to nuclease attack (e.g., the "split" nucleosomes observed by Lee
and Garrard [22]), (ii) the linker region of an
alternatively positioned downstream nucleosomal array in a subpopulation of cells, or (iii) the absence of a nucleosome between +90 and +280.
Between the upstream and downstream nucleosomal arrays on the active
allele in permeabilized cells was an approximately 350-bp region that
exhibited multiple MNase cleavages at nonnucleosomal intervals at
positions
140,
70, and +1 (Fig. 2A, lanes 14 to 18, and B, lanes 13 to 16). MNase cleavage at position
260 was also greatly enhanced
relative to that at all other sites on the active allele (in cells as
well as in naked DNA). The presence of cleavages at positions
140,
70, and +1 suggests that nucleosomes are not translationally
positioned over this region and, in conjunction with the strong MNase
cleavage at position
260, indicates an increased accessibility to
nucleases, consistent with the absence (or modification) of nucleosomes
in this region. This highly nuclease-accessible 350-bp region from
positions
260 to +90 contains the entire functional promoter
including the known transcription factor binding sites as well as the
multiple transcription initiation sites of the HPRT gene
(Fig. 2, diagrams). The absence of MNase cleavage at positions
140,
70, and +1 in naked DNA from 4.12 cells (Fig. 2A, lanes 10 to 13, and
B, lanes 9 to 12), even at higher MNase concentrations (Fig. 2A, lanes
10 to 13), suggests that these three cleavages in permeabilized cells
are dictated by the chromatin structure of the active allele rather
than the inherent susceptibility of the underlying DNA to MNase
cleavage. Significant cleavage within this 350-bp region was also
absent in permeabilized 8121 cells (Fig. 2A, lanes 5 to 9, and B, lanes
5 to 8) and naked DNA from 8121 cells (Fig. 2A, lanes 1 to 4, and B,
lanes 1 to 4), further supporting the notion that DNA in this region is
highly accessible in vivo on the active HPRT promoter.
DNase I hypersensitivity is characteristic of the HPRT
promoter on the active X chromosome, whereas the promoter on the
inactive X chromosome does not contain DNase I-hypersensitive sites
(26, 61). To examine the basis of the hypersensitivity of
the active promoter, we mapped DNase I-hypersensitive sites in the
promoter relative to the translationally positioned nucleosomes. Figure 2B (lanes 19 and 20) shows two major and three minor DNase
I-hypersensitive sites in chromatin of the active HPRT
promoter. All three minor hypersensitive sites, centered at positions
660,
1300, and
2700, comapped to complete or partial Alu repeat
elements upstream of the functional promoter. However, both of
the major DNase I-hypersensitive sites, centered at positions
260 and
70, map to the MNase-accessible 350-bp region corresponding to the
functional promoter on the active allele. The major DNase
I-hypersensitive site centered at position
260 mapped between the
AP-2 site and the cluster of CG boxes in the active HPRT
promoter. The other major DNase I-hypersensitive site, which spanned
the region from position
20 to
130, mapped immediately downstream
of the CG boxes and encompasses the major transcription initiation
sites as well as the potential initiator element (Fig. 1 and 2B, lanes
19 and 20). On the active promoter, the presence of both DNase I
hypersensitivity and increased MNase accessibility within the 350-bp
region suggests that this region is preferentially accessible to
nucleases and transcription factors and may be devoid of nucleosomes.
Furthermore, the locations of both the DNase I-hypersensitive sites (at
positions
260 and
70) and the MNase cleavages (at positions
260,
140, and
70) between transcription factor binding sites suggest
that transcription factor binding protects the underlying DNA from nuclease cleavage but may also distort the adjacent DNA, making it more
accessible for nuclease cleavage (20). In contrast, both
DNase I hypersensitivity (26) and MNase cleavage (Fig. 2A,
lanes 5 to 9, and B, lanes 5 to 8) are markedly absent within this
350-bp region on the inactive HPRT allele, consistent with the nuclease-inaccessible and transcriptionally silent nature of the
inactive promoter.
Overall, the data in Fig. 2 suggest that the promoter region of the
transcriptionally active HPRT allele is assembled into an
ordered array of translationally positioned nucleosomes interrupted by
a highly nuclease-accessible 350-bp interval that contains the known
functional elements of the HPRT promoter including the AP-2
site, the cluster of GC boxes, the potential initiator element, and the
region of multiple transcription initiation sites. Furthermore, three
CpG dinucleotides, at positions
97,
54, and
48, whose methylation
appears to be critical for transcriptional repression of the
HPRT promoter (6) all lie within this 350-bp
region. In contrast, nucleosomes on the inactive HPRT
promoter region do not appear to be translationally positioned, and
both hypersensitivity to DNase I and accessibility to MNase within this
350-bp region are reduced or absent.
Effects of DNA methylation on the translational positioning of
nucleosomes on the HPRT promoter in vitro.
To
examine the possibility that the differential methylation of the active
versus inactive HPRT promoters (11) might
mediate the differences in their nucleosomal organizations, a
Drosophila chromatin assembly extract (1) was
used to assemble nucleosomal arrays onto methylated and unmethylated
supercoiled human HPRT promoter templates. Templates were
methylated in vitro with HhaI, HpaII,
SssI, or HhaI plus HpaII methylase
prior to chromatin assembly. The reconstituted chromatin was then
digested with MNase, and the positions of MNase cleavage sites relative
to a BamHI site within the first intron of the
HPRT gene were mapped by Southern blotting and indirect end
labeling using a radiolabeled oligonucleotide probe
(BamHINuc1Probe; see Materials and Methods) corresponding to
the sequence immediately upstream of the BamHI site (Fig.
1). When the reconstituted chromatin was treated with MNase and not subsequently digested with BamHI (Fig.
3, lanes 1, 2, 5, 6, 9, 10, 13, 14, 17, and 18, labeled uncut), an approximately 180-bp ladder (i.e., 180, 360, 540 bp, etc.) of MNase-cleaved fragments was observed for all
templates, indicating that all templates formed regularly spaced
nucleosomal arrays. However, when these reconstituted MNase-treated
templates were subsequently digested with BamHI to map the
positions of the MNase cleavage sites relative to the intronic
BamHI site (Fig. 3, lanes 3, 4, 7, 8, 11, 12, 15, 16, 19, and 20, labeled BamHI), this 180-bp ladder disappeared in
all templates regardless of the methylation status of the template. This loss of the MNase-generated 180-bp ladder in samples digested with
BamHI suggests that nucleosomes in the reconstituted arrays are not translationally positioned (i.e., phased) in the same manner on
each template molecule. This is in direct contrast to the strong
translational positioning of nucleosomes observed on the active
HPRT promoter in vivo (i.e., in permeabilized cells; Fig.
2). Since the cleavage pattern generated by MNase treatment followed by
BamHI digestion is unaffected by the methylation status of
the template (Fig. 3, lanes labeled BamHI), DNA methylation does not appear to directly affect the translational positioning of
nucleosomes on the HPRT promoter in vitro. This suggests
that differential DNA methylation is unlikely to play a major role in
the differential translational positioning of nucleosomes on the active
versus inactive promoters in vivo. Furthermore, the absence of a
detectable 180-bp ladder of fragments in in vitro-reconstituted samples
that were MNase treated and BamHI digested suggests that the
HPRT promoter region does not contain a strong
nucleosome-positioning sequence. Interestingly, the positions of the
MNase cleavages represented by the 100-, 180-, 250-, and 340-bp
fragments seen in vitro (Fig. 3, lanes labeled BamHI)
coincide with the major hypersensitive MNase cleavage sites (at
positions +1,
70,
140, and
260) in the 350-bp nuclease-accessible
region on the active HPRT promoter in permeabilized cells
(Fig. 2A, lanes 14 to 18, and Fig. 2B, lanes 13 to 16).
|
DNase I in vivo footprinting of the human HPRT promoter. To further examine the correlation between the locations of the nuclease cleavage sites and the positions of transcription factor binding sites in the active functional HPRT promoter in permeabilized cells, DNase I in vivo footprint analysis was performed, concentrating on the region containing the DNase I-hypersensitive sites. While the pattern of DMS in vivo footprints (i.e., the pattern of guanine contacts in DNA by sequence-specific DNA-binding factors) in the HPRT promoter has been previously determined (13), DNase I in vivo footprinting allowed examination of the extent of DNA-protein contacts at each transcription factor binding site and permitted detection of the binding of additional transcription factors that do not interact with DNA at guanine residues. In vivo DNase I footprint analysis was performed by treating permeabilized 4.12 and 8121 cells with increasing concentrations of DNase I, purifying the DNase I-treated genomic DNA, and examining the DNase I cleavage pattern in the region of interest by LMPCR. The strand and region covered by each of the LMPCR primer sets used in this analysis are shown in Fig. 1.
DNase I footprint analysis of the minimal promoter in permeabilized cells identified two major footprints on the active HPRT allele, which corresponded to the potential AP-2 binding site and the cluster of GC boxes (Fig. 4). The footprint over the GC boxes extended from position
152 to
220 on
the upper strand (Fig. 4A, lanes 2 to 4) and from position
154 to
230 on the lower strand (Fig. 4B, lanes 12 to 14). The footprint over
the AP-2 binding site spanned positions
255 to
272 on the upper strand (Fig. 4C, lanes 23 to 25) and
260 to
275 on the lower strand
(Fig. 4B, lanes 12 to 14). This study complements previous DMS in vivo
footprinting studies (13), which identified footprints over the AP-2 binding site and the cluster of GC boxes, by further defining the region of DNA occupied and protected by the bound proteins. It also confirms that DNase I hypersensitivity and increased MNase cleavage appear to occur between rather than within transcription factor binding sites in this region. DNase I in vivo footprinting identified no new footprints in the functional promoter that were not
previously detected by DMS in vivo footprinting over the region analyzed (from positions
10 to
350).
|
95 to
86 that overlaps the two major transcription initiation sites
of the active promoter and that lies near a critical methylation site
at
97 (6). DNase I in vivo footprinting of the 4.12 human-hamster hybrid containing the active X chromosome did not detect
an equivalent footprint (Fig. 5, lanes 3 to 5 and 14 to 16), but an in vivo DNase footprint consisting of a
DNase I-protected region from positions
106 to
83 was observed in
the human male fibrosarcoma cell line, HT1080 (data not shown). This
observation suggests that the interactions between the hamster
trans-acting factor(s) and the human initiator element are
readily detectable by DMS but not DNase I footprinting in 4.12 cells,
while the interactions of the human factor(s) with the human sequence
in HT1080 cells are possibly more stable and therefore more readily
detectable by DNase I in vivo footprinting.
|
Rotational orientation of nucleosomes at the HPRT promoter. The high-resolution DNase I cleavage pattern generated by DNase I footprinting can also identify the rotational positioning of nucleosomes in chromatin. DNase I binds within the minor groove of DNA and can cleave DNA wrapped around the histone octamer, but the accessibility of the minor groove to DNase I varies as a function of the helical path of the DNA along the surface of the histone octamer. DNase I is thought to preferentially cleave the DNA where the minor groove faces directly away from the histone octamer and therefore is maximally exposed to the solution. When a nucleosome is rotationally positioned within a population of cells, this differential accessibility of the DNA helix wrapped around the histone octamer results in a pattern of preferential DNase I cleavages at approximately 10-base intervals corresponding to the helical pitch of DNA. This cleavage pattern can be visualized as a ladder of bands with a periodicity and spacing of approximately 10 bases in a DNA sequencing gel (23, 31, 38, 46).
Such a 10-base periodicity was, in fact, observed for DNase I cleavage of the inactive HPRT promoter of permeabilized 8121 cells (Fig. 4A, 5, and 6). This periodic 10-base cleavage pattern was the result of both enhanced DNase I cleavage at 10-base intervals and suppression of cleavage between bands of the 10-base ladder. Figure 5A shows a strong periodic 10-base cleavage ladder which extends from position
75 to
167 on the lower strand of the inactive promoter
(lanes 6 to 8) and which covers the region of multiple transcription
initiation sites and the potential initiator element. A similar pattern
of preferential cleavages at 10-base intervals also occurred on the
upper strand, extending upstream, from positions
154 to
249,
through the cluster of GC boxes on the inactive allele (Fig. 4A, lanes
5 to 7, and Fig. 6B, lanes 18 to 20), and downstream, from positions
107 to
38 (Fig. 5B, lanes 17 to 19), where the cleavages encompassed the positions of all three
critical sites of methylation located at
97,
54, and
48
(6).
|
|
284 to
387 on the
lower strand of the active promoter (Fig. 6A, lanes 4 to 6). This
rotationally positioned region occurs within a translationally
positioned nucleosome identified by MNase analysis of the active
promoter (from positions
260 to
460), suggesting that this
nucleosome is both translationally and rotationally positioned.
However, a similar cleavage pattern was also observed in the naked DNA
from 4.12 cells (Fig. 6A, lanes 1 to 3), making the
interpretation of rotational positioning in this region somewhat uncertain. The pattern of DNase I cleavage sites indicative of rotational positioning on the inactive HPRT promoter is
summarized in Fig. 7.
| |
DISCUSSION |
|---|
|
|
|---|
Analysis of both translational and rotational nucleosomal
positioning at the HPRT promoter indicates that the
nucleosomal organization of the promoter differs dramatically between
the active and inactive X chromosomes. A schematic summary of these differences is shown in Fig. 8. The
active promoter is assembled into an ordered array of six upstream and
at least two downstream translationally positioned nucleosomes which
flank a 350-bp nuclease-accessible region that includes the minimal
promoter (42) and the region of multiple transcription
initiation sites (16, 32) in the HPRT gene.
Consistent with previous DMS in vivo footprinting studies (13), DNase I footprinting analysis of permeabilized 4.12 cells (i.e., the active promoter) shows that the potential AP-2 site and the cluster of GC boxes that reside in this region are occupied on
the active promoter; no new footprints were identified by DNase I
footprinting on either the active or inactive HPRT promoter. This 350-bp region on the active allele appears to be devoid of nucleosomes since it is hypersensitive to DNase I (Fig. 2B, lanes 19 and 20), is accessible to MNase cleavage (Fig. 2A, lanes 14 to 18, and
B, lanes 13 to 16), and resembles naked DNA in its high-resolution
DNase I cleavage pattern (Fig. 4 to 6).
|
In contrast, no translational positioning of nucleosomes on the
inactive allele was observed by in vivo MNase analysis (Fig. 2A, lanes
5 to 9, and B, lanes 5 to 8). This absence of translational nucleosomal
positioning on the inactive promoter is evident both in the lack of a
200-bp periodic MNase cleavage pattern in vivo (Fig. 2A, lanes 5 to 9, and B, lanes 5 to 8) and the similarity between the patterns of MNase
cleavage of chromatin of the inactive allele (Fig. 2A, lanes 5 to 9, and B, lanes 5 to 8) and naked DNA (Fig. 2A, lanes 1 to 4, and B, lanes
1 to 4). Unlike the active allele, the inactive allele exhibits no
DNase I hypersensitivity (26), and MNase accessibility is
greatly reduced, as indicated by the significant reduction or absence
of MNase cleavages at positions
260,
140,
70, and +1 (Fig. 2A,
lanes 5 to 9, and B, lanes 5 to 8). Consistent with this finding, the
AP-2 site and the cluster of GC boxes on the inactive promoter appear
to be unoccupied by sequence-specific DNA binding proteins, as
determined by both DNase I footprinting and previous DMS in vivo
footprinting analysis (13). However, DNase I cleavage
analysis of permeabilized cells did identify rotational positioning of
nucleosomes over a region of at least 210 bp, from immediately
downstream of the potential AP-2 binding site through the cluster of GC
boxes and the region of multiple transcription initiations site to
about 40 bp upstream of the translation start site (Fig. 7 and 8). This region of rotational positioning on the inactive promoter is
particularly strong over the region of multiple transcription
initiation sites and falls completely within the 350-bp
nuclease-accessible region of the active allele. The location of this
rotational positioning and its unique association with the inactive
promoter argue that it may be involved in maintaining the
transcriptional repression of the inactive HPRT allele.
This pattern of nucleosomal organization on the active and inactive HPRT alleles is unlike those of two other endogenous promoters and one enhancer that have been examined for both rotational and translational positioning in vivo, specifically, the mouse albumin enhancer, the human pS2 promoter, and the mung bean phaseolin promoter (23, 30, 46). These promoters and the enhancer exhibit translational positioning of nucleosomes on the active allele, but only the albumin enhancer and the HPRT promoter contain a region that appears to be free of nucleosomes. On the inactive allele, the pS2, phaseolin, and HPRT promoters all exhibit rotational positioning of nucleosomes while the albumin enhancer does not. However, the inactive pS2 and phaseolin promoters also exhibit translational positioning of the rotationally positioned nucleosomes, whereas the inactive HPRT promoter does not. These differences in the nucleosomal organization of various promoters and the enhancer suggest that different nucleosomal architectures can be utilized to facilitate both transcriptional activation and transcriptional silencing.
The inactive HPRT promoter appears to be the first promoter described to demonstrate rotational positioning of nucleosomes in the absence of translational positioning (either on the active or inactive allele), indicating that translational and rotational positioning of nucleosomes are not only conceptually distinct but also physically distinct and separable modes of nucleosomal organization.
Translational positioning of nucleosomes and the active HPRT promoter. While a causal relationship between nucleosomal organization and transcription state cannot be established by these studies, it is likely that the different nucleosomal architectures of the active and inactive HPRT promoters play a significant role in promoter function. The strict translational positioning of nucleosomes over the HPRT promoter may act to preferentially expose specific transcription factor binding sites and the multiple transcription initiation sites to facilitate transcriptional activation. This hypothesis is consistent with several studies that indicate that the promoters and enhancers of expressed genes are assembled into translationally positioned nucleosomal arrays in which the nucleosomes over the transcription initiation site are absent or modified (2, 33, 46, 51). In contrast, the randomly positioned nucleosomes on the inactive promoter may serve to obstruct access to critical cis-acting elements by their cognate transcriptional activators and interfere with transcription initiation since transcriptional initiation (but not elongation) is strongly inhibited by nucleosomal assembly (28).
The mechanisms responsible for establishing and maintaining the translational positioning of nucleosomes in vivo are not well understood, but both nucleosomal positioning sequences and boundary proteins have been implicated (25). Unlike a number of promoters (which are predominantly inducible promoters) (29, 39, 40, 45, 46), the HPRT promoter region does not appear to harbor a nucleosome-positioning sequence. In vitro reconstitution of nucleosomes on the promoter region using a Drosophila chromatin assembly extract (1) does not generate a uniform translationally positioned nucleosomal array (Fig. 3). In addition, the inactive HPRT promoter region is not organized into translationally positioned nucleosomes in vivo (Fig. 2A, lanes 5 to 9, and B, lanes 5 to 8), in contrast to promoters that contain a translational positioning sequence and invariably show translational positioning of nucleosomes on both the active and inactive alleles (29, 39, 40, 45, 46). The role of boundary proteins in the translational nucleosomal positioning of the active HPRT promoter is less clear. Numerous DNA-binding proteins, when bound to their cognate site, have been shown to translationally position nucleosomes both in vitro (8, 19, 34, 36, 37, 47, 55) and in vivo (2, 30, 47). The position of the potential AP-2 binding site (from positions
86 to
94) in the HPRT promoter immediately
adjacent to the first upstream translationally positioned nucleosome in
the active promoter (between positions
260 and
460; Fig. 8)
suggests that the protein bound to the AP-2 site may act as a boundary
element to maintain the positioning of the translationally positioned
nucleosome array upstream of the AP-2 site. However, Litt et al.
(27) have shown that, during the reactivation of the
HPRT gene by 5aCdr, chromatin remodeling from a
nuclease-resistant to a nuclease-accessible conformation occurs prior
to the detectable binding of transcription factors to the promoter
region (see below). Even if the protein bound to the AP-2 site does not
mediate the initial remodeling and opening of chromatin during
5aCdr-induced reactivation, it may nevertheless be involved in
positioning the first nucleosome of the upstream array and preventing
nucleosomes from assembling on or sliding into the highly
nuclease-accessible 350-bp region in the active allele.
Because differential DNA methylation of the active and inactive
HPRT promoters is a prominent feature of the HPRT
locus and because demethylation of the promoter by 5aCdr is closely
associated with the remodeling of chromatin in the promoter region
(27, 44), it was conceivable that DNA methylation could
play a role in establishing and/or maintaining the differential
nucleosomal architecture of the active and inactive HPRT
promoter regions. However, our analysis of in vitro chromatin
assembly on methylated and unmethylated HPRT promoter
templates (Fig. 3) suggests that DNA methylation has little, if any,
role in facilitating or preventing the translational positioning of nucleosomes.
Results of in vivo analysis of the chromatin structure of the X-linked
PGK-1 promoter on the active and inactive X chromosomes appear to be somewhat inconsistent with our findings for the
HPRT promoter. Pfeifer and Riggs (38) reported
that the human PGK-1 promoter region on the inactive X
chromosome is organized into nucleosomal arrays that are
translationally positioned. This conclusion was based on the
observation that DNase I generated 10-base cleavage ladders, similar to
those found in the HPRT promoter region, on the inactive
allele in permeabilized cells. While these ladders may be indicative of
rotational positioning of nucleosomes, they do not necessarily indicate
translational positioning. No MNase analysis or other direct
examination of translational nucleosomal positioning in the
PGK-1 promoter was performed, so it is conceivable that
nucleosomes are rotationally positioned but not translationally positioned in the inactive PGK-1 promoter (as we find on the
inactive HPRT promoter).
Rotational positioning of nucleosomes and the inactive HPRT promoter. To date, the promoters of very few genes have been examined for rotational positioning of nucleosomes. Thus, the prevalence of rotational positioning and its significance for the regulation of transcription in vivo are neither well understood nor well documented. Of the four other genes examined, the inactive human PGK-1 (38), mung bean phaseolin gene (23), and human pS2 (46) promoters show rotationally positioned nucleosomes, while the inactive mouse albumin enhancer does not. However, rotational positioning is modified, but not abolished, in the active pS2 and phaseolin genes. Therefore, specific patterns of rotationally positioned nucleosomes are strongly associated with transcriptionally repressed promoters, though it is not entirely clear in each case if rotational positioning contributes to transcriptional silencing or is the result of silencing.
If rotational positioning of nucleosomes contributes to transcriptional repression, it might modulate the accessibility of transcription factors to their cognate sites in DNA since the binding of certain transcription factors appears to be sensitive to the rotational orientation of their binding sites relative to the histone octamer surface (14, 24, 59). Thus, rotational positioning of nucleosomes in the HPRT promoter may inhibit transcription factor binding to one or more of the GC boxes (particularly those at
212 and
198, which are located predominantly between periodic DNase
I cleavage sites; Fig. 7) or to a potential transcriptional initiator
element associated with the two major transcription initiation sites
(13, 16). It is also possible that the rotational orientation of specific sequences inward toward the surface of the
histone octamer hinders the binding of transiently bound factors (including chromatin-remodeling factors) that are not detectable by in
vivo footprinting.
Alternatively, the rotational positioning of nucleosomes in the
HPRT promoter may insure that a nucleosomal core, rather
than a linker region, is preferentially associated with the
transcription initiation sites of the promoter. For the HPRT
promoter, the region that exhibits the most distinct 10-base cleavage
periodicity (from positions
167 to
75) also coincides almost
perfectly with the multiple sites of transcriptional initiation for the
HPRT promoter reported by Patel et al. (32) and
Kim et al. (16), extending from positions
89 to
168.
If a nucleosomal core (rather than a linker region) is preferentially
situated over the region of transcriptional initiation, it would likely
inhibit transcription initiation (28).
The basis for rotational nucleosomal positioning is not well understood
and has not been closely examined but may involve the underlying DNA
sequence (23, 46, 50) and/or epigenetic modifications such
as DNA methylation (7, 50). For the HPRT promoter, the underlying sequence alone is unlikely to influence rotational positioning since the active promoter does not exhibit rotational positioning. However, since the active and inactive HPRT promoters are differentially methylated in vivo
(11), DNA methylation may be involved in allele-specific
rotational nucleosomal positioning (though it does not appear to affect
translational positioning; Fig. 3).
5-Azadeoxycytidine-induced chromatin remodeling of the HPRT promoter. Litt et al. (27) have shown that, during the course of HPRT gene reactivation on the inactive X chromosome by 5aCdr, chromatin remodeling of the promoter region from a nuclease-resistant to a nuclease-accessible conformation occurs prior to both transcription factor binding and the appearance of mRNA. Presumably, the 5aCdr-induced remodeling of the promoter region consists of converting the nucleosomal organization on the inactive promoter, where nucleosomes are not translationally positioned across the promoter region but rather are rotationally positioned over key regulatory regions, to a transcriptionally competent structure that contains translationally positioned nucleosomal arrays flanking a nucleosome-free region. The findings of Litt et al. suggest that initiation of this remodeling of the nucleosomal architecture in response to 5aCdr treatment does not require the stable binding of sequence-specific transcription factors and that the remodeling may in fact facilitate the subsequent transcription factor binding by exposing transcription factor binding sites in the promoter. Thus, factors responsible for remodeling the nucleosomal architecture of the HPRT promoter region during 5aCdr-induced reactivation may bind only transiently and are subsequently replaced by stably bound sequence-specific DNA-binding transcription factors that initiate HPRT gene transcription. This scenario is supported by recent findings of Kontaraki et al. (17) that indicate that changes in chromatin structure in the lysozyme promoter region during development are initiated at developmental stages well before end stage sequence-specific transcriptional activators bind to the promoter and activate transcription of the gene. Conversely, one allele of the HPRT gene also must undergo the reverse remodeling of nucleosomal architecture in the promoter region, from the active nucleosomal organization to an inactive organization, during the process of X chromosome inactivation in female embryogenesis. Little is known about the mechanisms or factors that mediate these changes in nucleosomal organization in vivo, but repression of the HPRT promoter (and, presumably, maintenance of a repressive chromatin conformation) appears to involve methylation of specific critical CpG sites in the promoter region (6).
Overall, these data suggest that the nucleosomal organization of the promoter may play an important role in activation or repression of transcription at the HPRT locus. On the active allele, the translational positioning of nucleosomes preferentially exposes the functional promoter, potentially allowing efficient transcription, whereas, on the inactive allele, rotationally positioned nucleosomes may sterically hinder transcription factor binding and formation of the preinitiation complex. The role of DNA methylation and transcription factor binding in establishing and maintaining this differential nucleosomal organization of the HPRT promoter on the active and inactive X chromosomes is currently under investigation.| |
ACKNOWLEDGMENTS |
|---|
We thank Jorg Bungert and Peter Becker for providing Drosophila chromatin assembly extracts and Jorg Bungert and Kelly Leach for helpful advice in the nucleosome reconstitution experiments.
This work was supported by NIH grant RO1 GM44286 to T.P.Y.
| |
FOOTNOTES |
|---|
* Corresponding author. Mailing address: Dept. of Biochemistry and Molecular Biology, Box 100245 JHMHC, University of Florida College of Medicine, 1600 SW Archer Rd., Gainesville, FL 32610. Phone: (352) 392-6472. Fax: (352) 392-2953. E-mail: yang{at}cmg.health.ufl.edu.
| |
REFERENCES |
|---|