| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|
| ||||||||||||||||||||||||||||||||||
National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, Maryland.
Address correspondence to Bruce Howard, MD, National Institute for Child Health and Human Development, NIH, Bldg. 31, Rm. 2A25, 6 Center Dr., Bethesda, MD 20892. E-mail: howard{at}helix.nih.gov
| Abstract |
|---|
|
|
|---|
Epigenetic states that specify either reduced or elevated transcription competence are well studied in yeast and Drosophila melanogaster (1013). In these genetically tractable organisms, epigenetic control is dependent on histone modifications and histone variants, as well as interacting proteins that form heritable chromatin-based structures. The genomes of vertebrates and plants also reflect regulation by DNA methylation-dependent mechanisms (14); thus, to examine epigenome patterns fully in vertebrates, both DNA methylation-based and chromatin-based regulatory arms should be considered.
Global aspects of epigenetic control have been repeatedly investigated in the contexts of postnatal development and aging. Most early studies in mammals failed to reveal robust changes in chromatin composition or acetylation of histones (15); however, clear age-related increases are documentable in trimethylation of histone H4 at lysine 20 in rat tissues (16). Also, histone replacement variants are more abundant in the livers of adult than young mice (17). Small age-related losses of global DNA methylation can be documented in mouse and human tissues (18), and more striking decreases occur with passage in nonimmortalized cells (8). Interestingly, telomerase expression can increase total genomic 5-methyl CpG content (19).
Numerous publications describe locus-specific alterations in DNA methylation patterns with age. In some instances, experiments focus on one or a small number of CpG dinucleotide positions (20,21). More extensive age-related DNA methylation appears to occur within mouse ribosomal rRNA gene elements (proximal 5' spacer domains) (22). A series of papers by Heidmann and colleagues documents age-related activation of intracisternal A particle-associated transcripts with concomitant loss of DNA methylation in mouse liver (23). Particularly influential are observations, originally by Issa and Baylin, that age-related CpG methylation occurs at some CpG islands (24). In an excellent brief overview, Issa stresses the increasing mosaicism in DNA methylation patterns with age (25). This investigator emphasizes the stochastic nature of the process, estimating that the level of abnormality ranges from 5% to 50% (though most often near the lower value). Overall, such alterations in DNA methylation patterns appear to fall within the framework of cumulative age-related disorder.
In comparison with the body of work on DNA methylation, little is known at the molecular level concerning age-related chromatin changes in vertebrates. This contrasts with the yeast model of replicative senescence, where gene expression changes consistent with chromatin remodeling have been documented (26). Likewise, progress has been rapid in delineating differentiation-associated remodeling (2734), where genes encoding specialized differentiation functions provide a natural target for study. The task of finding development-related and age-related chromatin alterations is challenging, in part because such processes may be subtle and occur within seemingly unchanging cell statesin designing experimental approaches, it is simply not obvious in what gene-regulatory or gene-encoding regions to look. One strategy is to search for large-scale patterns by mapping microarray-based RNA data onto the genome (35). In previous work, we developed another solution to the epigenome search problem (36). A semirandom genome sampling/chromatin immunoprecipitation (gsChIP) approach was devised, validated, and applied to screen for differences in the histone H4 acetylation epigenome patterns of human fibroblasts. We looked for changes either related to in vitro passage or donor age. The searches revealed several loci subject to remodeling, including one located about 11 Mb from the 4p chromosome terminus and a second about 3 Kb from the 4q terminus.
An immediate question that arises on finding such loci is how far along the chromosome the observed acetylation changes might extend. To map a single region in two cell types is not difficult. However, if one envisages a systematic effort along these lines, it will entail the mapping of dozens to hundreds or more regions in numerous cell types, in multiple organisms, at various stages of aging, and so forth. A second major problem thus emerges, i.e., how to map flanking (or query potentially similar noncontiguous) regions efficiently in terms of time and effort. Mapping involves multiple steps, including primer design, testing of primers to verify that yields of a single product reflect template concentration, redundant amplification reactions from ChIP-derived genomic DNA fragments, and finally graphic display of the results for a substantial number of positions. We report here an integrated solution to these issues. Further, we propose from the results that new and unexpected links between development and aging may have been identified.
| METHODS |
|---|
|
|
|---|
Preparation of Nuclei and Chromatin
Nuclei and histone-H1 depleted chromatin were prepared by a procedure described previously (38). Micrococcal nuclease concentrations were titrated in order to obtain chromatin fragments consisting predominantly of mononucleosomes, dinucleosomes, and trinucleosomes. All solutions used starting from cell trypsinization through preparation of nuclei and chromatin contained 10 mM sodium butyrate to preserve histone acetylation.
Immunoprecipitation of Acetylated Chromatin
Antibody against acetylated H4 was prepared by immunizing rabbits with chemically acetylated histone H4 (39,40). Affinity purification was accomplished with acetylated BSA coupled to BrCN Sepharose and characterized by enzyme-linked immunosorbent assay (ELISA) and Western blotting as described previously (38). Affinity-purified antibody (70 µg) was incubated with magnetic beads coated with antirabbit immunoglobulin (IgG) (Dynal Biotech, Brown Deer, WI) as described by the manufacturer. Chromatin fragments (25 µg) were incubated with antibody-coated beads by rotating overnight at 4°C in 50 mM NaCl, 10 mM Tris-HCl pH 8.0, 1 mM EDTA, 10 mM Na-butyrate, and proteinase inhibitor cocktail (39,40). Supernatant from the first magnetic separation combined with the first wash was taken as the nonbound chromatin fraction. Beads were washed six times with 1 ml 150 mM NaCl, 10 mM Tris-HCl pH 8.0, 1 mM EDTA, 10 mM Na-butyrate, proteinase inhibitor cocktail (Roche Diagnostics, Indianapolis, IN). Samples were transferred to new tubes prior to a final wash. Bound DNA was recovered directly from the beads using Qiagen PCR (polymerase chain reaction) purification kits (Valencia, CA). DNA concentrations were determined fluorimetrically. The typical yield from 25 µg chromatin was 0.51 µg bound DNA.
PCR Amplification
Primer design across RepeatMasker (Smit, AFA & Green, P RepeatMasker at http://repeatmasker.org)-processed sequences was accomplished using an interface utility designed to work with Primer3 (41). Primers were obtained from Integrated DNA Technologies (Coralville, IA) in a normalized 96-well format. PCR amplification was performed in triplicate in a 20 µl reaction mix containing 0.51 ng DNA, 1X buffer (Perkin-Elmer, Wellesley, MA), 50 µM deoxynucleotides, 0.5 µM primers, 1 unit Taq polymerase (Boehringer-Mannheim/Roche). The reactions were set up using an automated liquid handling system (Multiprobe II HT, Perkin-Elmer) in 96-well plates (ABgene, Inc., Rochester, NY). The cycling parameters were: 5 minutes at 94°C, followed by 30 cycles of 30 seconds at 94°C, 30 seconds at the annealing temperature, and 30 seconds at 72°C. As a check for linearity, reactions containing 1.5, 0.5, or 0.17 ng DNA were performed initially for each primer pair. PCR products were fractionated by ion-pair reverse-phase high-pressure liquid chromatography (HPLC) on a polystyrene-divinylbenzene matrix at 50°C and detected by SYBR gold (WAVE system, Transgenomics, Inc., Omaha, NE). Gradients started at 56% buffer A (0.1 M triethylammonium acetate [TEAA]), 44% buffer B (0.1 M TEAA, 25% acetonitrile), and finished at 40.9% A, 59.1%B.
Real-time PCR was performed using SYBR Green Fast-start PCR kits and a LightCycler instrument (Roche). Each sample (0.51 ng) was analyzed in triplicate together with 0.13 ng input DNA used for standard curve calculation.
Computational Analysis
Bioinformatics software utilized included Primer3 (41), RepeatMasker, Tandem Repeats Finder (42), and TESS [TESS: Transcription Element Search Software on the WWW, Jonathan Schug and G. Christian Overton, Technical Report CBIL-TR-1997-1001-v0.0, Computational Biology and Informatics Laboratory, School of Medicine, University of Pennsylvania, 1997, URL: http://www.cbil.upenn.edu/tess]. Utilities developed in the course of this study provide a primer design interface for Primer3, generation of scalar vector graphics (SVG)-based virtual gel images, quantification of peak heights and associated standard deviations, and red-green-blue (RGB) graphic display of chromatin acetylation patterns [code available on request]. This study utilized the high-performance computational capabilities of the Biowulf PC/Linux cluster at the National Institutes of Health, Bethesda, MD (http://biowulf.nih.gov).
| RESULTS |
|---|
|
|
|---|
Rapid ChIP Analysis Using an HPLC-Based Approach
For primer design and testing, the strategy adopted was to evaluate numerous primer pairs and to discard those that failed to yield informative data without optimization. Primer design across a RepeatMasker-processed region of interest was accomplished using a console-based Unix utility that constructs primer pair lists by means of consecutive calls to a locally installed copy of Primer3. Software-generated output lists were suitably formatted for direct ordering of plate arrays containing normalized primer pairs. PCR reactions were set up using an automated liquid handling system and run in a 96-well format. Products of the reactions were sampled directly and fractionated by ion-pair reverse-phase HPLC on a poly (styrene-divinylbenzene) matrix, with detection by postfractionation staining with SYBR gold. Notably, the sensitivity of this system permits PCR amplifications to be minimized so as to remain in or close to the linear range. XML-format result files were exported and processed using a locally written utility to both quantify the results and convert them to an SVG-based virtual gel format (see Figure 1A).
|
A summary of the first-stage mapping results is provided in Figure 2A. As shown at the upper right of this figure, the region analyzed spans 2.15 Mb. A central subregion, denoted (b), is mapped over approximately 75 Kb in either direction flanking the original h9 locus, itself indicated by a black triangle. Fetal to adult ratios reach maximum values of 4.1 and 4.6 at two sites flanking h9 (gray triangles). Two more distal 1 Mb flanking regions, denoted (a) and (c), are sampled at intervals of about 40 Kb or greater. In the direction of the centromere, only minor fetal versus adult differences in acetylation are evident, whereas an extended region of acetylation difference is seen in the telomere-proximal direction.
|
Timeframe for Chromatin Change at a Position Close to the h9 Site
Several loci exhibiting relatively large acetylation changes were examined further to relate these in an approximate way to developmental stages. Chromatin was prepared from skin fibroblasts at comparable low passage levels and from the following sources: additional fetal (3 samples), neonatal foreskin (2 samples), very young [ages 1 and 2 yr; (2 samples)], and young adult [ages 24 to 30 yr; (4 samples)]. Two positions closely flanking the original h9 locus (denoted by grey triangles in Figure 2A) displayed moderate acetylation only in fetal samples, with unchanging low acetylation in other cells (not shown). One position 85 Kb upstream of the MIST geneindicated by asterisks in Figures 2A, 2B, and 3exhibited a more gradual decrease in acetylation. In this case, baseline low acetylation levels appear to be reached by or just after the period of early childhood (Figure 3). No further acetylation change was observed on comparison of samples from young and old adults (unpublished observations).
|
While further mapping remains to be done, a summary of the results is assembled in Figure 4. Near the centromere-proximal edge of the region, an area of high histone H4 acetylation is evident. As this is not associated with a CpG island or an identified gene, it could be an insulator, or have another function not yet discernible. Similarly to the 4p16.1 region, patchy subregions of acetylation change are seen over a near megabase-size chromosome segment. Young versus old adult ratio values are indicated by empty triangles at selected positions. There appears to be a trend toward higher ratio values as mapping and sampling proceed in the direction of the telomere. A histone H4 acetylation ratio of 4.3 near LOC389250, a hypothetical gene supported by expressed sequence tags (ESTs), is evident, and, like the subregion upstream from the MIST gene, will be of interest to examine in greater detail.
|
| DISCUSSION |
|---|
|
|
|---|
Mapping and sampling data for the 4p16.1 region extend over a distance of about 2 Mb. Most histone H4 acetylation changes occur within a roughly 0.7 Mb subregion delimited by clusters of 5 telomere-proximal and 12 centromere-proximal positions that exhibit little or no remodeling. Two pairs of genes are situated near or overlapping the edges of the 0.7 Mb subregion. Of particular interest is the MIST gene (LOC166522), since the most striking remodeling spans the 5' part of the gene and extends through the promoter region. MIST/Clnk and related proteins are adaptors involved in signal transduction, mostly studied in hematopoietic lineages (44,45). A murine homologue of the MIST gene is expressed in embryonic stem cells, but reverse transcriptase (RT)-PCR analysis revealed no expression in fetal fibroblasts (unpublished observations).
Yet to be determined are what sequence features in the MIST upstream region contribute to its striking age-related chromatin remodeling. As a first step to explore this issue, we generated a 160 Kb subregion map, and this revealed three more discrete areas where the magnitude of the histone H4 acetylation decrease reaches 8-fold to 11-fold. Each peak area of change still spans 8 Kb or more, so a third round submapping step should be informative. Possibly significant is the fact that the MIST promoter region is devoid of CpG islands. The nearest islands are located between this gene and KIAA1729 (near the 5' end of the latter), and close to the 5' end of LOC153027. Most gene-associated CpG islands are highly acetylated, and it is plausible that large islands might stabilize neighboring chromatin domains against development-related or age-related chromatin remodeling. Consistent with this, drift in CpG island methylation patterns is typically stochastic and occurs at low levels (25).
In mapping from the MIST gene toward the 4p telomere, a span of about 240 Kb is encountered, which exhibits minimal histone H4 acetylation change. At the most telomere-proximal position sampled, however, decreasing acetylation is again seen. Could this be the edge of another domain of chromatin remodeling? Imprinted domains with opposite maternalpaternal polarity are, in some instances, clustered into much larger regions. Quite possibly, the same may be true for the remodeling domains of the type described here. A very unusual tandem array of 4.7 Kb repeats, each copy of which contains a deubiquitinating enzyme 1A-like gene, is located approximately 0.8 Mb telomere-proximal to the edge of the presently mapped 4p16.1 domain. This large array merits comment, as tandem arrays are involved in several examples of striking age-related epigenetic change. At least two reports describe progressive DNA methylation within the tandemly arrayed ribosomal RNA genes (22,46). Another documents strong age-related silencing of transgene repeat arrays (47).
At the other end of the same chromosome, the 4q35.2 remodeling region has been subjected to preliminary mapping. While the boundaries of acetylation change are yet to be determined, the ZFP42 gene and several additional uncharacterized transcription units fall within the region already examined. Three of the transcription units are associated with minimal CpG islands. Several tandem arrays are to be found, including one located centromere-proximal at a distance of about 340 Kb, and a second located less than 50 Kb from the 4q telomere. Deletions of the latter array are clearly implicated in the etiology of facioscapulohumeral muscular dystrophy (48). The proximity of relatively large tandem arrays to both 4p16.1 and 4q35.2 remodeling regions may be coincidental; nevertheless, it is conceivable that silencing mechanisms acting on such arrays may exert clock-like influences on chromatin states of adjacent single-copy regions.
Beyond the above two domains, are there numerous other genome regions subject to age-related chromatin remodeling? And, if so, do any regions extend across and influence critical genesi.e., those already linked to important diseases and pathophysiologies of human aging? This remains to be seen, but using the techniques described here, it should be straightforward to search for and map other examples. At the current basic level of automation, mapping across regions of megabase size can be accomplished inexpensively (relative to microarrays) and in a span of several days. Simple expedients (e.g., fractionation of bands from multiple PCR amplification reactions in a single column run), as well as further HPLC development and more-sensitive staining-detection technologies should increase throughput considerably. Moreover, when full genome-based microarrays become widely accessible, this approach should be of use for fine mapping of bacterial artificial chromosomes (BAC) or larger size domains identified by hybridization.
What significance should be ascribed to age-related changes in large chromatin domains? At this point, we cannot exclude the possibility that such chromatin remodeling is an epiphenomenon, of little or no consequence in the aging process. On the other hand, the remodeling of epigenetic patterns is a topic of broad interest relevant to both differentiation-associated gene regulation and the genome reprogramming required for successful animal cloning experiments. Further, it has been hypothesized that age-related change in chromatin-based epigenetic states could serve as one of the basic driving forces for aging (49,50). The results reported in this study are predicted by, and appear to provide important new experimental support for, that hypothesis.
Finally, while this article is focused solely on pattern changes in the epigenome, it is important to emphasize that the links between epigenetic mechanisms and aging are likely to be numerous and intricate. Histone-modifying and transcription factor-modifying enzymes are typically embedded in large (and presumably multifunctional) complexes. Disturbances in the functions of these complexes or DNA methylation mechanisms modulate cell growth and senescence, mortality pathways, and organismal aging (5156). While such disturbances may prove to accelerate region-specific alterations in epigenome patterns, links can likewise be postulated to multiple damage-related pathways, metabolism, and other modulators of aging.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received August 3, 2004
Accepted October 1, 2004
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. P. de Magalhaes and G. M. Church Genomes Optimize Reproduction: Aging as a Consequence of the Developmental Program Physiology, August 1, 2005; 20(4): 252 - 259. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||
| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|