|Home | About | Journals | Submit | Contact Us | Français|
Immunolabeling of metaphase chromosome spreads can map components of the human epigenome at the single cell level. Previously, there has been no systematic attempt to explore the potential of this approach for epigenomic mapping and thereby to complement approaches based on chromatin immunoprecipitation (ChIP) and sequencing technologies.
By immunostaining and immunofluorescence microscopy, we have defined the distribution of selected histone modifications across metaphase chromosomes from normal human lymphoblastoid cells and constructed immunostained karyotypes. Histone modifications H3K9ac, H3K27ac and H3K4me3 are all located in the same set of sharply defined immunofluorescent bands, corresponding to 10- to 50-Mb genomic segments. Primary fibroblasts gave broadly the same banding pattern. Bands co-localize with regions relatively rich in genes and CpG islands. Staining intensity usually correlates with gene/CpG island content, but occasional exceptions suggest that other factors, such as transcription or SINE density, also contribute. H3K27me3, a mark associated with gene silencing, defines a set of bands that only occasionally overlap with gene-rich regions. Comparison of metaphase bands with histone modification levels across the interphase genome (ENCODE, ChIP-seq) shows a close correspondence for H3K4me3 and H3K27ac, but major differences for H3K27me3.
At metaphase the human genome is packaged as chromatin in which combinations of histone modifications distinguish distinct regions along the euchromatic chromosome arms. These regions reflect the high-level interphase distributions of some histone modifications, and may be involved in heritability of epigenetic states, but we also find evidence for extensive remodeling of the epigenome at mitosis.
Large scale projects are underway to map the epigenomes of various eukaryotes, including humans. The objective is usually to define the distribution across the genome of modified histones, various non-histone proteins or methylated cytosines, and then link these modifications to genomic functions [1-3]. Genome-wide analyses have been made possible by coupling the long-established technique of chromatin immunoprecipitation (ChIP) with either high density DNA microarrays (ChIP-chip) or next-generation DNA sequencing (ChIP-seq) . These powerful technologies require material from large numbers of cells and the data generated inevitably represent a mean value derived from cells with differing patterns of expression from a significant subset of genes. Differences can arise through intrinsic transcriptional noise or because cells are in different phases of the cell cycle. Such cell to cell heterogeneity inevitably limits the precision with which histone modifications can be linked to chromatin function. In principle, this issue can be addressed by using immunomicroscopy to examine the distribution of histone modifications at the single cell level. Metaphase chromosome spreads provide a source of material in which individual chromosomes can be identified and in which the entire human epigenome can be scanned in a single cell. This approach has several additional advantages: there is little or no transcription at metaphase, removing a major source of variability between cells, consistency from cell to cell can be monitored, fluorescent probes are extremely sensitive (offering detection at the single gene level if required) and the procedure is quick (once experimental conditions are established) and relatively cheap. It should also be noted that immunostaining, if properly controlled, can detect modified histones and other proteins across the entire genome, including repeat-rich regions that are inaccessible to sequencing-based approaches . While microscopy cannot match the ultimate resolving power of ChIP-seq, it has the potential to provide a valuable complementary approach to epigenomic mapping.
Immunolabeling of metaphase chromosomes is a well established technique and has revealed dramatic regional differences in the distribution of specific histone modifications, particularly the distinctive pattern of modifications present on centric (constitutive) heterochromatin in plants and animals [5-7] and the facultative heterochromatin of the inactive X chromosome in female mammals [8,9]. Immunolabeling of meiotic (pachytene) chromosomes in maize has shown regional variation in levels of various methylated histone isoforms, with distinctive differences between heterochromatin and euchromatin .
Surprisingly, there has been only limited use of metaphase chromosome immunostaining to map histone modifications across individual chromosomes [11,12], and no systematic attempt to explore the genome-wide distribution of multiple histone modifications.
Here we describe a systematic analysis of the distribution of selected histone modifications across metaphase chromosomes from normal human cells. Antibodies to histone modifications previously linked to active transcription (H3K9ac, H3K27ac and H3K4me3, described collectively as active modifications) all highlight the same 10- to 50-Mb genomic regions, giving a characteristic and consistent banding pattern. Bands closely correspond to regions rich in genes and CpG islands (CGIs). In contrast, H3K27me3, a mark associated with gene silencing, shows a preference for telomeric regions and defines bands that only occasionally overlap with gene rich regions. At 10-Mb resolution, active modifications have similar, though not identical, distributions across interphase  and metaphase chromosomes, while H3K27me3 is clearly different. The results suggest that there is extensive remodeling of the epigenome as cells enter mitosis, but that a high-level memory of some components of the interphase epigenome is retained into metaphase.
Standard protocols for preparation and staining of metaphase chromosomes require fixation in acidified organic solvents, a step that extracts the great majority of histones and other proteins . We have adopted an approach using unfixed chromosomes [9,15,16], a procedure that has the major advantage that histones remain in their native (that is, unfixed, undenatured) form and are therefore structurally compatible with the synthetic peptides used to raise anti-histone antisera [17,18]. We found that both the relative sizes and centromeric indices (arm ratios) of unfixed chromosomes were very similar to their counterparts fixed in methanol/acetic acid (Additional files 1 and 2), allowing us to use these properties as a first step in chromosome identification. Unfixed chromosomes are not amenable to conventional G-banding procedures. To distinguish morphologically similar chromosomes, we used the chromosome-specific banding patterns generated by the DNA counterstain DAPI (4,6-diamino-2-phenyl-indole). DAPI selectively stains regions that are AT-rich and GC-poor, and gives a banding pattern that resembles G-banding and is unique for each chromosome .
Unfixed metaphase chromosome spreads from human lymphoblastoid cells were immunostained with antibodies to histone H3 tri-methylated at lysine 4 (H3K4me3), a modification that has been associated with transcriptionally active, or potentially active, chromatin [18-21]. Centromeric heterochromatin was consistently unstained, while the arms of most chromosomes showed a characteristic pattern of brightly stained and weakly stained regions (Figure 1a, b). Using a combination of size, centromeric index and reverse DAPI banding (Figure (Figure1c),1c), we were able to identify all chromosomes and construct karyotypes (Figure 1d, e). There was consistently strong staining of both arms of chromosome 19, weak staining of chromosome 13 and distinctive banding of most chromosomes, with particularly prominent bands on chromosomes 1, 6, 9, 11 and 12. The immunofluorescent banding pattern was consistent between sister chromatids and homologues and reproducible from one spread to another, despite the inevitable differences in overall chromosome size. Alignments of chromosomes from five immunostained spreads are shown in Additional file 3.
Very similar immunostaining patterns were given by antisera to two other modifications also loosely associated with transcriptionally active chromatin, namely H3 acetylated at lysine 27 (H3K27ac) and H3 acetylated at lysine 9 (H3K9ac) [22,23] (Figure (Figure2a;2a; Additional files 4 and 5). Conversely, staining with a variety of antisera to acetylated H4 was more uniform. The acetylated H4 bands corresponded to those seen with antisera to H3K4me3 but the differential labeling of bands and interband regions was less extreme. A typical example is shown in Figure Figure2c.2c. H4K8ac is absent from both constitutive (centric) and facultative heterochromatin and our findings are generally consistent with previous studies on acetylated H4 [10,13].
H3 tri-methylated at lysine 27 (H3K27me3) is put in place by the methyltransferase Ezh2, a component of the Polycomb silencing complex PRC2 and has been associated with formation of facultative heterochromatin and gene silencing [24-26]. In female cells, one of the two X chromosomes generally stained more strongly than its homologue, and more strongly than the single X in male cells (Figure (Figure2b;2b; Additional file 6). The more intensely stained X is likely to be the inactive homologue . H3K27me3 was undetectable on blocks of constitutive centric heterochromatin (Figure (Figure2b;2b; Additional file 6) or on the Y heterochromatin in male cells (Additional file 7). There are distinctive regional variations in H3K27me3 staining intensity along the chromosome arms, but without the sharply defined banded distribution typical of H3K4me3 (Figure (Figure1).1). We find only limited overlap between the two modifications. For example, the short arm of chromosome 6 is relatively enriched in both modifications, but on closer inspection H3K27me3 has a more telomeric location (6pter-22.3) than H3K4me3, which is centrally located in the short arm (centered at 6p21), leaving the telomeric region relatively weakly stained (compare the multiple examples of chromosome 6 in Additional files 3 and 6). Also, the prominent H3K4me3 band on chromosome 11q just below the centromere (11q12.1-13.3) is not enriched in H3K27me3 (Figure (Figure2b).2b). Overall, we find that H3K27me3 is consistently enriched at telomeric regions, at least on the larger chromosomes (chromosomes 1 to 15). This distinctive staining pattern was seen with two different antisera to H3K27me3 (listed in Additional file 8). H3K27ac is a modification that may act as an antagonist of Polycomb-mediated silencing through suppression of H3K27 tri-methylation [4,24]. While the distribution of H3K27ac (Figure (Figure2a)2a) is clearly different from that of H3K27me3 (Figure (Figure2b),2b), H3K27me3 is not consistently excluded from regions rich in H3K27ac, or vice versa.
Immunostaining with antibodies to H4 tri-methylated at lysine 20 (H4K20me3) strongly and selectively labeled the centric heterochromatin of metaphase chromosomes from human lymphoblastoid cells (Figure (Figure2d),2d), consistent with previous results in other cell types . Absence of staining of centric heterochromatin by antisera to the other histone modifications tested here is clearly not due to a general inaccessibility of histone epitopes in heterochromatin. Chromosome arms were essentially unstained by antibodies to H4K20me3, with the exception of the Y chromosome in male cells, on which heterochromatic regions on the distal long arm were consistently stained (Figure (Figure2d2d).
Over the course of the work presented here, complete immunostained karyotypes for H3K4me3, H3K9ac, H3K27ac and H3K27me3 have been constructed from lymphoblastoid cell lines (LCLs) derived from two different individuals, one male and one female. At the present level of resolution, we have found no evidence for individual differences in chromosome banding. The same banding patterns have also been seen in occasional preparations from two other LCLs (results not shown). Analyses of other cell types have been less extensive, but immunostaining of chromosomes from human primary fibroblasts with antibodies to H3K4me3 revealed a banding pattern essentially the same as that seen in LCLs (Additional file 9). The banding patterns described are not restricted to a particular cell lineage. However, differences may occur among more widely divergent, or aberrant, cell types. Improved resolution of immunofluorescent bands, perhaps through analysis of extended, prophase chromosomes, may also reveal differences not apparent with the present approach.
Recent analyses have confirmed that most genes are clustered in a relatively small number of genomic regions [28-30]. These gene-rich regions are also enriched in CGIs, relatively CpG-rich DNA sequences found at and around the promoter regions of many genes and characterized by low levels of DNA methylation [31,32]. We constructed gene density/CGI maps for each human chromosome by calculating the gene and CGI content of 10-Mb windows across the chromosome. In Figure Figure3,3, the resulting histograms are aligned with a representative example of each chromosome immunostained for H3K4me3. There is a close and consistent correspondence between high levels of H3K4me3 and regions of relatively high gene/CGI content. This is true not only for regions of intense staining (for example, landmark bands on chromosomes 1q, 6p and 11q) but also for less strongly staining bands that do not stand out in the original spreads (for example, the bands distributed across chromosomes 3 and 12) (Figure (Figure1;1; Additional file 3). As expected from our earlier results, chromosomes immunostained with antibodies to H3K9ac and H3K27ac showed essentially the same close relationship between staining intensity and gene/CGI density (results not shown). In contrast, on chromosomes immunostained for H3K27me3, there was only limited overlap between gene/CGI-rich regions and staining intensity (Additional file 7).
To allow a quantitative analysis of the relationship between the distribution of histone modifications at metaphase and other chromosome properties, we measured the level of H3K4me3 across chromosome 1 by scanning. Typical scans of sister chromatids are shown in Figure Figure4a.4a. Fluorescence intensity is expressed as a percentage of the most highly fluorescent element and distance along the chromosome is expressed in megabases (chromosome 1 is 247 Mb long and we have assumed a linear relationship between positions on the metaphase chromosome and genomic DNA). To allow us to combine data from multiple scans, the chromosome was divided into 25 equal segments (each having a nominal length of 10 Mb) and the total fluorescence within each segment calculated. The fluorescence distribution (banding pattern) obtained by averaging scans from 12 chromosomes (24 chromatids) is shown in Figure Figure4b.4b. Comparison of these quantitative data with gene and CGI frequencies across chromosome 1, also grouped within 10-Mb windows (Figure (Figure3),3), shows that they are closely correlated (r = 0.70 and 0.68 respectively, P < 0.0002).
As a first step in exploring the link between H3K4me3 levels at metaphase and transcription in interphase, we used single color, high-density oligonucleotide arrays to measure transcript levels for 3,071 RefSeq genes across chromosome 1 in the same lymphoblastoid cells used for immunolabeling. Total transcript levels within 10-Mb windows across chromosome 1 are shown in Figure Figure4b.4b. There is a close correlation between interphase transcription and levels of H3K4me3 at metaphase, measured by immunofluorescence labeling (Figure (Figure4b;4b; r = 0.73, P < 0.00002). Both transcription and H3K4me3 immunofluorescence are strongest in regions of the chromosome depleted in major G bands (for example, 1pter-p33, 1q21-23; Figure 4a, b).
The genome-wide distribution of various histone modifications in asynchronous (mostly interphase) human lymphoblastoid cells has recently been defined by ChIP-seq  (see Materials and methods). The results can be aligned with immunostained metaphase chromosomes to provide an initial comparison of the interphase and metaphase epigenomes. Results for three modifications (H3K27ac, H3K4me3 and H3K27me3) on three chromosomes (chromosome 1, 6 and 11) are shown in Figure Figure5.5. At 10-Mb resolution, there is a close correspondence between the interphase and metaphase distributions of H3K27ac and H3K4me3, with clearly defined interphase peaks aligning with the major metaphase bands. The correspondence for H3K4me3 is particularly precise, with even the weakly stained double band on distal chromosome 1q evident in interphase (Figure (Figure5).5). Quantitative analysis using chromosome scanning data (Figure (Figure4)4) confirms the visual alignment of H3K4me3 levels across chromosome 1 at metaphase and interphase, with a strong correlation between them (r = 0.74, P < 0.00002; all pairwise correlations are presented in Additional file 10). In contrast, we find little correspondence between the distributions of H3K27me3 in interphase and metaphase. The chromosome-wide distribution of H3K27me3 in interphase at 10-Mb resolution is relatively homogeneous, the most prominent feature being its depletion across the block of centric heterochromatin on chromosome 1 (Figure (Figure5).5). There are no interphase peaks corresponding to the highly stained H3K27me3 bands present at metaphase.
Previous studies have shown that progression into mitosis is accompanied by an overall decrease in global histone acetylation levels, reduced acetate turnover and changes in the relative levels of acetylation at specific lysines [34,35]. In view of this, it is perhaps surprising that the high level distribution of histone acetylation across the interphase genome, as revealed by ChIP-seq, is retained in metaphase chromosomes (Figure (Figure5).5). A possible explanation comes from the finding that for both H3K27ac and H3K4me3, the differences between enriched and depleted regions are more extreme in metaphase chromosomes than in interphase chromatin. For example, the regions on chromosome 1p and 1q that lie between the brightly stained bands (distal 1p, proximal 1q) are virtually unstained and comparable to centric heterochromatin, a finding confirmed for H3K4me3 by quantitative scanning (Figure (Figure4a).4a). The equivalent regions at interphase show levels of modification well above that of centric heterochromatin (Figure (Figure5).5). While the different technologies used to derive the two sets of data may contribute to these differences, the comparison suggests that at least some histone modifications are preferentially removed from gene-poor chromosomal regions as cells enter mitosis.
H3 di-methylated at lysine 4 (H3K4me2) has been shown to be strongly enriched at promoters with the highest CpG content (CGI promoters), even when they are transcriptionally silent . It has been suggested that H3K4 methylation protects these promoters from silencing by CpG methylation, a proposition supported by in vitro experiments . In light of these findings, one could propose that H3K4me3 levels at metaphase are a simple reflection of CGI density. However, closer inspection of the chromosome labeling patterns suggests that banding is unlikely to be solely attributable to simple genomic features such as gene or CGI density. For example, the gene-rich, CGI-rich region 11q12.1-13.3 is consistently one of the most strongly stained regions in the genome with antisera to the three activating modifications tested. The region at 11pter-15.3 is similarly gene-rich and only slightly less CGI-rich (Figure (Figure3),3), yet stains much less strongly with antisera to H3K27ac and H3K4me3 (Additional files 3, 4 and 5). Another example is provided by the gene/CGI-rich regions across chromosome 12. The region on the q arm adjacent to the centromere labels with antisera to all three active modifications tested, but the labeling intensity is consistently less than, or at best equal to, labeling of the less gene/CGI-rich region on the distal q arm (Figure (Figure3;3; Additional files 3, 4 and 5). It is interesting to note that this strongly staining distal region has a higher density of short interspersed nuclear element (SINE) repeats (UCSC hg18 ) than the more gene-rich, centromere proximal region (Figure (Figure6).6). This is unusual because gene/CGI density and SINE density are very closely correlated across the genome (figures for chromosome 1 are shown in Additional file 10). On the basis of these examples, it could be argued that SINE density is more closely associated with levels of active histone modifications than gene/CGI density. This possibility is supported by the correlations derived from the chromosome 1 scanning data (Additional file 10).
Immunostaining of polytene chromosomes from the salivary glands of chironimid insects established the principle that levels of histone acetylation across the interphase genomes of higher eukaryotes show extreme regional variation, giving distinctive and reproducible immunofluorescent banding patterns [38,39]. Islands of acetylated histone H4 occurred within transcriptionally active and silent regions and within condensed (phase dense) and more open (phase light) chromatin, and were therefore not solely dependent on either transcription or chromatin compaction . In the absence of polytene chromosomes, it is only comparatively recently that the same principle has been shown to apply to the interphase genomes of mammals. By combining ChIP with a cloning strategy based on the serial analysis of gene expression (SAGE) technique, Roh et al.  identified over 46,000 regions enriched in H3 acetylated at lysine 9 and/or 14 in human T-cells. These regions, designated 'acetylation islands', were often associated with promoters, putative control elements and CGIs. At least some of the acetylated islands were dynamic features; activation of T cells with accompanying gene activation and chromatin remodeling resulted in the appearance of over 4,000 new islands and the disappearance of some pre-existing ones . There was a close correlation between the frequency of acetylated islands and gene density  and in a chromosome-by-chromosome presentation of the data (supplementary Figure Figure33 in ), regions of high acetylation (for example, on chromosomes 1q, 6p, 11q and 19) correspond to the brightly staining H3K9ac and H3K27ac metaphase bands presented here.
H3K27me3 also shows evidence of regional variation across the genome. An analysis of H3K27me3 across mouse chromosome 17 by ChIP-chip and application of a new algorithm for detecting broad regions of histone modification  showed that the modification tends to occur in large regions, designated BLOCs, of average size 43 kb. There are examples of H3K27me3 spreading across large domains in humans, where consistently high levels of H3K27me3 cover the 100- to 200-kb regions encompassing the human HOX gene clusters. At a higher level, H3K27me3 BLOCs were found to be more frequent in gene-rich, SINE-rich regions, along with high levels of H3 and H4 acetylation. The authors propose that these regions alternate across the chromosome with gene-poor, SINE-poor, long interspersed nuclear element (LINE)-rich regions with relatively high levels of H3K9me3 and H4K20me3, two histone modifications associated with constitutive heterochromatin [42,43]. As discussed by the authors, this model is not supported by mouse ChIP-seq data  analyzed in the same way, or with ENCODE data from human cell lines that showed no evidence for consistent co-localization of H3K27me3 and active histone modifications such as acetylated H3 and H4 and H3K4me3 . The data presented here show that in human metaphase chromosomes, H3K27me3 is preferentially located across defined regions of 10 Mb and above. These regions are not gene-rich, nor does H3K27me3 consistently co-localize with acetylated histones or H3K4me3. However, there is overlap between H3K27me3-rich and H3K4me3/H3K9ac/H3K27ac-rich regions (examples can be seen in Figure Figure5),5), showing that, at the highest level, the two chromatin types are not mutually exclusive. As yet we have not been able to align the H3K27me3 banding pattern with any genomic features. H3K27me3 bands do not correspond to the frequency of LINE repeats plotted as 10-Mb windows (results not shown), or to SINE and ALU repeats, which closely correlate, as expected, with gene/CGI density (Additional file 10).
The bands we describe are large, approximately 10 to 50 Mb, and presumably encompass many (perhaps several hundred) smaller chromatin domains, some associated with specific genes and gene clusters and their control elements. A crucial question is whether the bands have any functional significance in their own right, or whether they passively reflect the net level of histone modification among the subdomains that they contain. In assessing this, it is relevant that genes and their control elements make up only a small proportion of the chromatin within a band, with even the most gene-rich band having only approximately 30 genes/Mb. The histone modifications studied here are relatively common and therefore must be mostly located in intergenic chromatin. The difference in gene/CGI density between the most gene-rich and gene-poor domains at 10-Mb resolution is only about 6-fold and differences in repeat density are even less. It is questionable whether differences of this order are sufficient to account for the differences in staining intensity between bands and interbands, with the latter often essentially unstained (that is, comparable to centric heterochromatin). It is also interesting that the banding patterns given by three very different modifications (H3K4me3, H3K9ac and H3K27ac) are so similar. It may be that the banding given by H4K8ac, and other acetylated H4 isoforms, for which the difference in staining intensity between bands and interbands is less extreme, may be a closer reflection of gene/CGI density. It should also be borne in mind that, for some modifications at least, high-level chromosome banding may not be directly determined by DNA sequence elements but by other aspects of chromosome behavior. For example, if interphase chromosome territories are configured so that some regions are accessible to, or share a nuclear location with, subsets of histone modifying enzymes, then one would expect to see large chromosome domains displaying high levels of selected histone modifications, just as we observe.
Our results show that for active modifications (H3K4me3, H3K9ac and H3K27ac) immunofluorescent bands on metaphase chromosomes correspond to enriched regions in the interphase epigenome revealed by ChIP-based approaches [18,40]. The localized persistence of active modifications through mitosis may play a role in determining gene function in the following G1 phase, thereby contributing to the heritability of epigenetic states . This could be done by maintaining a general chromatin property, such as the open chromatin structure found across gene-rich regions of the chromosome . A close analysis of the human transcriptome map (HTM) by Versteeg and colleagues  showed that many highly expressed genes are clustered in about 40 gene-rich regions of 10 to 15 Mb in size, designated ridges. Weakly expressed genes tend to cluster in similarly sized, gene poor regions designated antiridges. A quantitative analysis of the properties of ridges and antiridges on chromosomes 1, 6 and 11 in six different human cell lines, all in G1 phase, showed that ridges were consistently less condensed, less spherical and further from the nuclear periphery than antiridges. These properties were not changed by the major differences between cell types in karyotype and gene expression pattern . Ridges often correspond in both position and extent to the metaphase chromosome bands rich in active histone modifications described here (the band at 60 to 70 Mb on chromosome 11p is a good example) [28,29]. It may be that the distribution of histone modifications at metaphase represents part of a mechanism by which the structural properties of gene-rich regions and subregions are maintained through mitosis. Large chromatin regions carrying specific combinations of modified histones could also help establish chromosome territories in the reforming G1 nucleus, perhaps to ensure optimum positioning of gene-rich chromatin .
The H3K27me3-rich bands on metaphase chromosomes appear to have no equivalents in interphase, indicating that they are generated by regional adjustment of this modification as cells progress through mitosis. In this context, it is interesting that the intensity of H3K27me3 staining, particularly of the smaller chromosomes such as 19 and 20, is more variable between cells than is that of the other modifications studied here (Additional file 6). Perhaps this is attributable to ongoing demethylation of H3K27 in metaphase cells, with the paler-staining chromosomes being derived from cells that had been blocked in metaphase for longer (up to 4 hours). As noted earlier, close comparison of interphase and metaphase distributions of active modifications suggests that here too there are targeted changes in modification levels, with a selective reduction in interband regions serving to enhance the banding pattern in metaphase. Taken together, these findings suggest that there is widespread remodeling of the epigenome during mitosis. The enzymatic mechanisms responsible for such remodeling can be investigated by selective disruption of candidate enzymes and it will be of particular interest to explore how these mechanisms might be subverted in disease states.
There is a characteristic distribution of histone modifications across the metaphase genome, giving each chromosome a distinctive immunofluorescent banding pattern. Bands of the active modifications H3K4me3, H3K9ac and H3K27ac are virtually indistinguishable, but differ from bands of H3K27me3. Bands are consistent between cell lines of the same type and, for H3K4me3 at least, between two different cell lineages. There is a close correlation between bands of active modifications and gene/CGI density, SINE density and transcript levels in interphase, though none of these parameters alone provides a complete explanation for the location and relative staining intensities of the different bands. The functional, or DNA sequence properties that correlate with the H3K27me3 bands remain mysterious. For H3K4me3, H3K9ac and H3K27ac, the metaphase banding resembles their distribution at interphase, whereas for H3K27me3 metaphase and interphase distributions are different. The results provide evidence of extensive remodeling of the epigenome as cells enter mitosis, even for modifications where the resemblance between interphase and metaphase distributions is clear.
Immortalized, Epstein Barr virus-transformed LCLs from healthy, karyotypically normal individuals were established in-house and maintained at 37°C, 5% CO2 in RPMI 1640 medium, 10% fetal bovine serum, supplemented with L-glutamine (2 mM) and penicillin/streptomycin (all additives from Gibco (Gibco, Grand Island, NY, USA)). These lines retain a normal diploid karyotype over many years in culture and provide a consistent source of experimental material. All results presented here are from two individuals, one male (line AH) and one female (line VM).
Rabbit polyclonal antisera used for immunolabeling are listed in Additional file 8. In-house antisera were prepared by immunization with synthetic peptides as previously described [35,48]. The specificity of all antisera was confirmed by inhibition ELISA using native histones immobilized on microtiter plates and a selection of synthetic peptides, as described previously .
Cells in exponential growth were treated for 2 hours with colcemid (KaryoMax, Gibco) at 0.1 μg/ml. Cells were pelleted by centrifugation at 1,200 rpm (Chillspin, MSE, London, UK) for 5 minutes at 4°C, washed twice with cold phosphate buffered saline, resuspended in 75 mM KCl at 2 × 105 cells/ml and left at room temperature for 10 minutes. Aliquots (200 μl) of the swollen cell suspension were spun onto glass slides at 1,800 rpm for 10 minutes in a Shandon Cytospin 4. Slides were then immersed for 10 minutes at room temperature in KCM buffer (120 mM KCl, 20 mM NaCl, 10 mM Tris/HCl pH 8.0, 0.5 mM EDTA, 0.1% Triton X-100). Immunolabeling was carried out for 1 hour at 4°C, as described previously , with antisera diluted 200- to 400-fold in KCM supplemented with 1 to 1.5% BSA (Sigma-Aldrich, Dorset, UK). For all labelings, the secondary antibody was fluorescein isothiocyanate (FITC)-conjugated goat anti-rabbit immunoglobulin (Sigma F1262) diluted 150-fold in KCM, 1% BSA. Slides were washed twice in KCM (5 minutes at room temperature), fixed in 4% (v/v) formaldehyde (10 minutes, room temperature), rinsed in deionized water and mounted in Vector Shield (Vector Lab, Peterborough, UK) supplemented with DAPI (Sigma) at 2 μg/ml.
Labeled slides were visualized on a Zeiss Axioplan 2 epifluorescence microscope and potentially suitable metaphases spreads captured using Smart Capture software (Digital Scientific, Cambridge, UK). The filmstrip obtained was screened for the best spreads and poor quality spreads discarded. The coordinates of each spread were recorded using an England Finder(tm) graticule. Each identified spread was recaptured at the West Midlands Regional Genetics Laboratory on a Zeiss microscope and images stored in the ISIS software (Metasystems, Altlussheim, Germany). Karyotyping was based on reverse DAPI staining, which gives a conventional G-banding pattern, using the standard chromosome identification criteria used for clinical diagnostic work. A printout of each karyotype was produced and each chromosome located in the spread and numbered. Immunofluorescent karyotypes were then reconstituted in Smart Type software (Digital Scientific) to reveal the distribution of histone modifications across identified chromosomes.
Twelve representative examples of chromosome 1 (that is, 24 chromatid arms) were scanned for fluorescence intensity using the GraphPolygon Extension of IPLab in SmartCapture. The centromere of each chromosome was manually marked by adding a line of blue pixels, and each chromatid arm was manually tracked from pter, to centromere, to qter. GraphPolygon data comprised color intensities averaged over the pixels on each segment perpendicular to the manually tracked medial axis of the chromatid (for a line of specified width, which is slightly less than the chromatid arm). Numerical pixel values were obtained for green (FITC), red (DAPI) and blue (manually annotated centromere). The relative longitudinal position of each segment of each chromatid was then normalized to the actual length of chromosome 1 (247 Mb). The segment data from each chromatid were first anchored at the centromere (120 Mb), whose position was determined by the maximal level of pixel staining in the blue channel. Color intensity of FITC and DAPI staining was normalized for each chromatid arm as a percentage of the maximum staining for the relevant color. The normalized longitudinal segment positions were then grouped into 25 10-Mb windows and the mean and standard deviation of FITC intensity and DAPI intensity were calculated for each window.
ChIP-seq data from the Bernstein laboratory at the Broad Institute of MIT and Harvard were downloaded from the ENCODE data coordination center at UCSC . Signal data for the lymphoblastoid line GM12878 was averaged for each 10-Mb window for comparison with our data. The CGI distributions presented in this paper are all calculated from the data recently generated and compiled by Illingworth et al.  using CXXC affinity purification and deep sequencing.
RNA was extracted and purified from log-phase lymphoblastoid cells using the RNeasy kit with DNase digestion (Qiagen, Crawley, West Sussex, UK), according to the manufacturer's instructions. cDNA was synthesized using the Superscript double-stranded DNA synthesis kit (Invitrogen, Paisley, UK), cleaned up by RNase A treatment (Invitrogen; 4 μg per 150 μl reaction) followed by phenol:chloroform extraction and ethanol precipitation. Samples were labeled with cy3, hybridized to a 12 × 135 k HD2 expression array (Roche Nimblegen, Madison, Wisconsin, USA); containing 3 probes per gene for 44,049 human genes) and scanned (GenePix 4000B, Molecular Devices, Sunnyvale, CA, USA). Data were extracted, normalized by quantile normalization and gene calls generated by the robust multichip average algorithm using Nimblescan (Roche Nimblegen). Normalized expression data from three biological replicates is available at the Gene Expression Omnibus , accession number [GEO:GSE24459]; the data cover the entire human genome and are limited to RefSeq gene sequences. For comparison with H3K4me 3 levels and gene/CGI density, transcription from RefSeq genes was summed within 10-Mb windows across chromosome 1.
BSA: bovine serum albumin; CGI: CpG island; ChIP: chromatin immunoprecipitation; DAPI: 4,6-diamino-2-phenyl-indole; ELISA: enzyme-linked immunosorbant assay; FITC: fluorescein isothiocyanate; LCL: lymphoblastoid cell line; LINE: long interspersed nuclear element; SINE: short interspersed nuclear element.
The authors declare that they have no competing interests.
ET, FM, JAH and BMT conceived and designed experiments; ET, FM, JAH and PP performed experiments; ET, FM, JAH, PP, LPO and BMT analyzed data; RSI, AMRT, VD and LPO contributed reagents, materials and analytical/technical expertise; ET, FM, JAH and BMT wrote the paper. All authors read and approved the final manuscript.
Figure showing relative lengths of human chromosomes with and without fixation in methanol acetic acid.
Figure showing centromeric indices of human chromosomes with and without fixation in methanol acetic acid.
Figure showing alignment of individual metaphase chromosomes immunostained for H3K4me3 from five different chromosome spreads.
Figure showing alignment of individual metaphase chromosomes immunostained for H3K27ac from six different chromosome spreads.
Figure showing alignment of individual metaphase chromosomes immunostained for H3K9ac from two different chromosome spreads.
Figure showing alignment of individual metaphase chromosomes immunostained for H3K27me3 from six different chromosome spreads.
Figure showing the correspondence between gene density, CpG island density and H3K27me3 levels across human metaphase chromosomes.
Table showing antibodies used for labeling and their origins.
Immunostained karyotype showing the distribution of H3K4me3 across human fibroblast chromosomes.
Table showing correlations between H3K4me3 levels across metaphase chromosome 1, gene-, CGI- and repeat-frequencies and various properties of the interphase epigenome.
We thank Rebecca Sie and Milan Fernando for skilful technical assistance, Thomas Jenuwein for generously supplying antisera, Rogier Versteeg for sharing data with us early in the project, Steve Kissane for advice on microarray hybridization protocols and Adrian Bird and Karl Nightingale for valuable comments on the manuscript. This work was supported by Cancer Research UK (programme grant C1015) and by the European Union via FP6 and the Epigenome Network of Excellence.