|
|
||||||||
Divisions of
1 Functional Genomics
2 Cardiology, Jichi Medical School, 3311-1 Yakushiji, Kawachigun, Tochigi 329-0498, Japan
3 Department of Molecular Biology, Cancer Research Institute, Sapporo Medical University, Sapporo 060-8556, Japan
4 CREST, Japan Science and Technology Agency, Saitama 332-0012, Japan
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Acetylation of histones is mediated by histone acetyltransferases (HATs) and takes place on the
-amino group of conserved lysine residues located in the NH2-terminal tail of core histones (Carrozza et al. 2003). Such histone modification is tightly linked to transcriptional regulation by either chromatin remodeling or the provision of binding sites for other proteins. HAT activity in cells is rapidly counteracted by the activity of histone deacetylases (HDACs) (Verdin et al. 2003), with the result that the turnover time of histone acetylation is as short as a few minutes (Waterborg 2002).
The importance of histone acetylation in the regulation of gene expression has been demonstrated for a variety of cellular processes including cell differentiation, cell cycle progression, DNA repair, and carcinogenesis (Kouzarides 1999; Yasui et al. 2003). Furthermore, a histone acetylation profile (HAP) of a certain type of acute myeloid leukaemia was shown to be directly related to the sensitivity of the malignancy to chemotherapeutic agents (Grignani et al. 1998; Lin et al. 1998). Elucidation of the HAP of human malignancies may therefore provide a basis for novel approaches to cancer treatment. This notion is supported by the observation that inhibitors of HDAC activity are effective in the treatment of haematological malignancies such as multiple myeloma and myelodysplastic syndrome (Gore et al. 2001; Catley et al. 2003). The means to establish a genome-wide HAP have been lacking, however (Gabrielli et al. 2002).
To characterize the HAP of any given cell or tissue type, we have now developed a new screening method, termed differential chromatin scanning (DCS). Application of the DCS method to a human gastric cancer cell line has readily identified hundreds of target fragments of histone deacetylase (HDAC), including a number of cancer-related genes such as tumour protein p73 (TP73), cell division cycle 34 (CDC34), and bone morphogenetic protein 7 (BMP7).
With its high fidelity, DCS enables us a genome-wide screening of DNA fragments embedded in histones with differential acetylation level between a given pair of samples.
| Results |
|---|
|
|
|---|
The DCS procedure is schematically shown in Fig. 1. We verified the screening ability of this method by identifying HDAC targets in the genome of the human gastric cancer cell line MKN28. The cells were treated with trichostatin A (TSA), a specific inhibitor of HDAC activity, for use as a tester sample, whereas MKN28 cells not subjected to TSA treatment were used as a driver. Inhibition of HDAC activity by TSA would be expected to increase the acetylation level of target histones compared with that apparent in the driver sample. We thus attempted to rapidly isolate genomic fragments bound to histones that were acetylated only in the tester sample.
|
The nonspecific binding of residual RNA was minimized by treating the DNA solution with RNase A, and the DNA fragments were then rendered blunt-ended. The DNA was digested maximally with RsaI to obtain fragments with a relatively uniform size of several hundred base pairs. A TAG adaptor was ligated to both ends of the DNA fragments, and subsequent PCR amplification with a TAG primer yielded amplicons with an XmaI/SmaI site at each end.
The tester DNA was then digested with XmaI (thereby generating cohesive ends), whereas the driver DNA was digested with SmaI (generating blunt ends). The tester DNA was ligated with the first subtraction adaptor (Toyota et al. 1999) through its cohesive ends and was then annealed to an excess amount of the driver DNA. Under this condition, DNA fragments present only in the tester sample undergo self-annealing and thereby generate a binding site for the first subtraction primer at both ends. Subsequent PCR amplification with this primer thus selectively amplified the tester-specific DNA fragments.
To exclude DNA fragments that possess endogenous (probably nonspecific) binding sites for the first subtraction primer, we digested the first subtraction products with XmaI and ligated the resulting molecules with the second subtraction adaptor. A second round of subtraction amplification was then performed with the second subtraction primer, yielding DNA fragments that were associated with acetylated histones specifically in the tester cells.
HDAC targets in a cancer cell line
From the products of the second round of subtraction PCR, we randomly selected 288 DNA fragments and determined their nucleotide sequence. The size of the genome fragments was 50 bp in 265 clones (mean size,
270 bp), and these DNA sequences were subjected to subsequent analyses. The sequences were first screened with the BLAT search program (Kent 2002) against the nucleotide sequence database assembled as of July 2003 by the Genome Bioinformatics Group of the University of California at Santa Cruz (http://genome.ucsc.edu/). Among the 265 DNA fragments examined, 200 sequences showed > 95% identity to the human genome sequence and 198 of these were located either within a protein-coding gene (demonstrated or predicted) or in the vicinity (within 5 kbp) of such a gene (Table 1 and Table S1).
|
TP73 as a target of HDAC
Interestingly, certain of these 32 clones corresponded to sequences either within or close to human cancer-related genes, including those for TP73, CDC34, and BMP7. One of the clones (10sBKR_D10) was, for instance, mapped to a region located
1.2 kbp downstream of the last exon of TP73 (GENBANK accession number, NM_005427
[GenBank]
) (Fig. 2A). TP73 is a close homologue of a well-known tumour suppressor, TP53. The TP73 gene is localized at 1p36 locus which is subject to recurrent loss of heterozygosity in various human cancers (Kaghad et al. 1997). An array of isoforms for TP73 are generated in vivo through alternative splicing mechanism of its mRNA, and those forms exert a variety of (sometimes, opposite) actions on tumour growth (Benard et al. 2003).
|
Other targets of HDAC
One of our clones (10sBKR_E11) was mapped to a region which contains the first intron and the second exon of the granzyme M (GZMM) gene (GENBANK accession number, NM_005317
[GenBank]
), and is also located
5 kbp downstream of the last exon of CDC34 (GENBANK accession number, NM_004359
[GenBank]
) (Fig. 3A). Therefore, it was possible that histone acetylation of this region may affect the transcription of both genes. We first quantified the tester/driver ratio for the relative amount of the 10sBKR_E11 fragment in the acetylated histone H3 immunoprecipitates. As shown in Fig. 3B, histone at this region was more heavily acetylated in the tester cells compared to the driver cells.
|
CDC34 encodes a ubiquitin E2 ligase, and its activity is essential for G1- to S-phase transition in cell cycle (Pagano et al. 1995). CDC34 catalyses ubiquitination of the inhibitors for cyclin-dependent kinases (CDKI) which are thereby subject to proteolysis. Therefore, it was interesting to find that CDC34 gene was a target of HDAC in a cancer cell line.
We selected seven additional genes for quantification of the relative amounts of their mRNAs. Of the total of nine genes examined, seven genes were preferentially expressed in the TSA-treated cells compared with the non-treated cells (tester/driver ratio 1.5) (Table 1), consistent with the notion that histone acetylation results in the recruitment of transcription factors and consequent transcriptional activation.
Among the 136 clones that mapped within or close to known genes, only eight pairs were assigned to the same genes. Given the high fidelity of our DCS screening, these data suggest that HDAC complexes simultaneously act at hundreds of independent target histones within the same cell. To visualize directly the genome-wide HAP, we mapped to human chromosome figures our genomic clones whose chromosomal positions are known (Fig. 4). The HDAC targets were distributed widely throughout the human genome, although some hot spots for deacetylation were apparent. Seven of our DCS clones were, for instance, mapped to the same 8q24.3 locus. Detailed mapping indicated that all these clones were located within a region spanning 2.5 Mbp. Therefore, it may be possible that regional alterations of chromatin structure lead to a coordinated transcriptional regulation of genes within the affected region.
|
| Discussion |
|---|
|
|
|---|
Recently, similar genome-wide comparison methods have been developed by coupling chromatin immunoprecipitation (ChIP) with either comparative genomic hybridization (CGH) or analysis with a DNA microarray containing thousands of CpG islands or promoter fragments (Weinmann et al. 2002; Ballestar et al. 2003; Odom et al. 2004). However, the CGH analysis is not able to provide information on the precise genome sequences differentially immunoprecipitated, and the CpG array can only identify genome fragments that contain CpG islands. Interestingly, the majority of our clones identified by DCS do not satisfy the criteria of CpG islands (for instance, none of the clones corresponding to TP73, CDC34-GZMM, or BMP7 was rich in the CpG sequence). Additionally, screening of whole promoter fragments is still hampered by the incomplete human genome annotation and by its inability to investigate the enhancer regions. The emerging number of non-coding RNAs also makes it difficult to precisely map promoter sequences in human genome. Therefore, for the identification of HAT-HDAC targets, our DCS method is likely to provide a non-biased means of genome-wide screening.
Although the study in this manuscript has revealed only the HDAC targets in a cancer cell line, the DCS method should be useful to directly identify HAT targets as well. For instance, a DCS comparison between chemotherapy-sensitive cancer cells and chemotherapy-resistant ones would enable to isolate genes whose associated histones are acetylated in a chemoresistance-dependent manner. In addition, DCS method can be applied to other non-malignant cells/tissues. We could indeed identify HDAC targets in cardiomyocytes (R.K., personal communication).
Furthermore, modification of our method may open up the possibility of a general ChIP-based subtraction strategy. For instance, DCS performed with antibodies to a given histone- or DNA-binding protein, rather than with antibodies to acetylated histone, should readily allow the isolation of target genomic fragments differentially recognized by such proteins.
| Experimental procedures |
|---|
|
|
|---|
MKN28 cells were obtained from the Japanese Collection of Research Bioresources (http://cellbank.nihs.go.jp/) and were maintained in RPMI 1640 medium (Invitrogen) supplemented with 10% foetal bovine serum. For preparation of the tester sample, cells were incubated for 24 h with 300 nM TSA (Wako, Tokyo, Japan). Satoh et al. (2002) have already demonstrated that, in MKN28 cells, histones of several genes become acetylated in a TSA-dependent manner.
The DCS method
Both formaldehyde treatment of cells and immunoprecipitation with antibodies to acetylated histone H3 were performed with the ChIP assay kit (#06-599; Upstate Biotechnology). Sonication of DNA was performed in one 10-s pulse to minimize the extent of DNA shearing (the resulting DNA fragments should migrate as a smear corresponding to a size of 1020 kbp in gel electrophoresis). Immunoprecipitation was achieved with protein A-Sepharose beads (Sigma) suspended in TE (10 mM Tris-HCl, 7.4, and 0.1 mM EDTA) containing 0.2 mg/mL tRNA (Roche Biochemicals) and 0.5 mg/mL bovine serum albumin (New England Biolabs), instead of with the protein A beads supplied with the kit. DNA fragments were recovered from the immunoprecipitates and treated with RNase A to remove residual RNA. The DNA fragments were digested with RsaI (New England Biolabs), and ligated to the TAG adaptor (5'-CCACCGCCATCCGAGCCTTTCTGCCCGGG-3'/3'-GAAAGACGGGCCC-5'). After PCR amplification with the TAG primer (5'-CCACCGCCATCCGAGCCTTTCTGC-3'), the tester DNA and the driver DNA were digested with XmaI and SmaI, respectively. Only the tester DNA (0.5 µg) was ligated to the first subtraction adaptor (5'-GTGAGGGTCGGATCTGGCTGGCTC-3'/3'-CGACCGAGGGCC-5'). The adaptor-ligated tester DNA was then annealed with 40 µg of the driver DNA at 67 °C for 2024 h before PCR amplification with the first subtraction primer (5'-GTGAGGGTCGGATCTGGCTGGCTC-3'). After digestion of single-stranded DNA with mung bean nuclease (New England Biolabs), the amplified products were subjected to digestion with XmaI followed by a second round of subtraction PCR with the second subtraction adaptor (5'-GTTAGCGGACACAGGGCGGGTCAC-3'/3'-GCCCAGTGGGCC-5') and second subtraction primer (5'-GTTAGCGGACACAGGGCGGGTCAC-3'). The final products were digested with XmaI and ligated into pBlueScript (Stratagene) for nucleotide sequencing. The detailed protocol for the DCS method is available as a supplemental protocol online.
DNA quantification
Genome fragments immunoprecipitated by antibodies to acetylated histone H3 were subjected to PCR with a QuantiTect SYBR Green PCR Kit (Qiagen). The amplification protocol comprised incubations at 94 °C for 15 s, 62 °C for 30 s, and 72 °C for 60 s. Incorporation of the SYBR Green dye into PCR products was monitored in real time with an ABI PRISM 7700 sequence detection system (PE Applied Biosystems), thereby allowing determination of the threshold cycle (CT) at which exponential amplification of PCR products begins. The CT values for DNA molecules corresponding to the GAPDH gene and genome fragments of interest were used to calculate the abundance of the latter fragments relative to that of the former. The oligonucleotide primers for PCR were 5'-TCGGTGCGTGCCCAGTTGAACC-3' and 5'-ATGCGGCTGACTGTCGAACAGGAG-3' for GAPDH, 5'-ACCAGTCGCTGCTTTTAAATAAGG-3' and 5'-GCTAGGAGCTTCCCCTACTAACAT-3' for 10sBKR_D10, and 5'-CTGATGTGCGTTTTGAAGGACT-3' and 5'-GAAACTTTTAGGCTGATACGTGTG-3' for 10sBKR_E11.
mRNA quantification
Total RNA was prepared from the tester or driver cells with an RNeasy Mini column (Qiagen), treated with RNase-free DNase (Qiagen), and subjected to reverse transcription with PowerScript reverse transcriptase (BD Biosciences Clontech) with an oligo(dT) primer. Portions of the resulting cDNAs were subjected to PCR with a QuantiTect SYBR Green PCR Kit. The amplification protocol comprised incubations at 94 °C for 15 s, 65 °C for 30 s, and 72 °C for 60 s. The annealing temperature of PCR was changed to 60 °C for the cDNA of CDC34 and 66 °C for that of GZMM. The oligonucleotide primers for PCR were 5'-AGAACATCATCCCTGCCTCTACT-3' and 5'-ATATTTGGCAGGTTTTTCTAGACG-3' for GAPDH, 5'-GACGAGGACACGTACTACCTTCA-3' and 5'-GTAGGTGACTCGGCCTCTGTAG-3' for TP73, 5'-AAGATGTGGCACCCTAACATCTAC-3' and 5'-AGGGAGATCACACTCAGGAGAAT-3' for CDC34, and 5'-ATGTGTAACAACAGCCGCTTCT-3' and 5'-CTTGAAGATGTCAGTGCAGACC-3' for GZMM.
| Supplementary material |
|---|
|
|
|---|
| Acknowledgements |
|---|
| Footnotes |
|---|
* Correspondence: E-mail: hmano{at}jichi.ac.jp
| References |
|---|
|
|
|---|
Baylin, S.B. & Herman, J.G. (2000) DNA hypermethylation in tumorigenesis: epigenetics joins genetics. Trends Genet. 16, 168174.[CrossRef][Medline]
Benard, J., Douc-Rasy, S. & Ahomadegbe, J.C. (2003) TP53 family members and human cancers. Hum. Mutat. 21, 182191.[CrossRef][Medline]
Burgess-Beusse, B., Farrell, C., Gaszner, M., et al. (2002) The insulation of genes from external enhancers and silencing chromatin. Proc. Natl. Acad. Sci. USA
99 (Suppl. 4), 1643316437.
Carrozza, M.J., Utley, R.T., Workman, J.L. & Cote, J. (2003) The diverse functions of histone acetyltransferase complexes. Trends Genet. 19, 321329.[CrossRef][Medline]
Catley, L., Weisberg, E., Tai, Y.T., et al. (2003) NVP-LAQ824 is a potent novel histone deacetylase inhibitor with significant activity against multiple myeloma. Blood
102, 26152622.
Gabrielli, B.G., Johnstone, R.W. & Saunders, N.A. (2002) Identifying molecular targets mediating the anticancer activity of histone deacetylase inhibitors: a work in progress. Curr. Cancer Drug Targets 2, 337353.[CrossRef][Medline]
Gore, S.D., Weng, L.J., Zhai, S., et al. (2001) Impact of the putative differentiating agent sodium phenylbutyrate on myelodysplastic syndromes and acute myeloid leukemia. Clin. Cancer Res. 7, 23302339.
Grignani, F., De Matteis, S., Nervi, C., et al. (1998) Fusion proteins of the retinoic acid receptor-alpha recruit histone deacetylase in promyelocytic leukaemia. Nature 391, 815818.[CrossRef][Medline]
Kaghad, M., Bonnet, H., Yang, A., et al. (1997) Monoallelically expressed gene related to p53 at 1p36, a region frequently deleted in neuroblastoma and other human cancers. Cell 90, 809819.[CrossRef][Medline]
Kent, W.J. (2002) BLATthe BLAST-like alignment tool. Genome Res.
12, 656664.
Kouzarides, T. (1999) Histone acetylases and deacetylases in cell proliferation. Curr. Opin. Genet. Dev. 9, 4048.
van Leeuwen, F. & Gottschling, D.E. (2002) Genome-wide histone modifications: gaining specificity by preventing promiscuity. Curr. Opin. Cell Biol. 14, 756762.[CrossRef][Medline]
Lin, R.J., Nagy, L., Inoue, S., Shao, W., Miller, W.H. Jr & Evans, R.M. (1998) Role of the histone deacetylase complex in acute promyelocytic leukaemia. Nature 391, 811814.[CrossRef][Medline]
Odom, D.T., Zizlsperger, N., Gordon, D.B., et al. (2004) Control of pancreas and liver gene expression by HNF transcription factors. Science
303, 13781381.
Pagano, M., Tam, S.W., Theodoras, A.M., et al. (1995) Role of the ubiquitin-proteasome pathway in regulating abundance of the cyclin-dependent kinase inhibitor p27. Science 269, 682685.
Satoh, A., Toyota, M., Itoh, F., et al. (2002) DNA methylation and histone deacetylation associated with silencing DAP kinase gene expression in colorectal and gastric cancers. Br. J. Cancer 86, 18171823.[CrossRef][Medline]
Toyota, M., Ho, C., Ahuja, N., et al. (1999) Identification of differentially methylated sequences in colorectal cancer by methylated CpG island amplification. Cancer Res.
59, 23072312.
Verdin, E., Dequiedt, F. & Kasler, H.G. (2003) Class II histone deacetylases: versatile regulators. Trends Genet. 19, 286293.[CrossRef][Medline]
Waterborg, J.H. (2002) Dynamics of histone acetylation in vivo. A function for acetylation turnover? Biochem. Cell Biol. 80, 363378.[CrossRef][Medline]
Weinmann, A.S., Yan, P.S., Oberley, M.J., Huang, T.H. & Farnham, P.J. (2002) Isolating human transcription factor targets by coupling chromatin immunoprecipitation and CpG island microarray analysis. Genes Dev.
16, 235244.
Yasui, W., Oue, N., Ono, S., Mitani, Y., Ito, R. & Nakayama, H. (2003) Histone acetylation and gastrointestinal carcinogenesis. Ann. NY. Acad. Sci. 983, 220231.[CrossRef][Medline]
Received: 21 August 2004
Accepted: 16 September 2004
This article has been cited by other articles:
![]() |
R. Kaneda, S. Ueno, Y. Yamashita, Y. L. Choi, K. Koinuma, S. Takada, T. Wada, K. Shimada, and H. Mano Genome-Wide Screening for Target Regions of Histone Deacetylases in Cardiomyocytes Circ. Res., August 5, 2005; 97(3): 210 - 218. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | ADVANCED SEARCH | TABLE OF CONTENTS |