A carrier-assisted ChIP-seq method for estrogen receptor-chromatin interactions from breast cancer core needle biopsy samples
© Zwart et al.; licensee BioMed Central Ltd. 2013
Received: 5 September 2012
Accepted: 22 March 2013
Published: 8 April 2013
The Estrogen Receptor alpha (ERα) is the key transcriptional regulator in luminal breast cancer and is therefore the main target for adjuvant treatment of this subtype. Luminal gene signatures are dictated by the transcriptional capacities of ERα, which are a direct consequence of the receptors binding preference at specific sites on the chromatin. The identification of ERα binding signatures on a genome-wide level has greatly enhanced our understanding of Estrogen Receptor biology in cell lines and tumours, but the technique has its limitations with respect to its applicability in limited amounts of tumour tissue.
Here, we present a refinement of the ChIP-seq procedures to enable transcription factor mapping on limited amounts of tissue culture cells as well as from a limited amount of tumor tissue derived from core needle biopsies. Our approach uses a carrier that can be removed prior to DNA amplification and sequencing.
We illustrate the applicability of this refined technology by mapping the ERα genome-wide chromatin binding landscape in core needle biopsy material from primary breast tumours. With this, our refined technology permits for a high-resolution transcription factor mapping even from clinical samples.
Breast cancer is the most common diagnosed malignancy in women, with over 1.4 million new cases worldwide annually . 75% of all breast cancers are from the luminal subtype  for which tumour growth is thought to be dependent on the activity of the Estrogen Receptor alpha (ERα). This growth-dependency renders this pathway the main target for adjuvant endocrine treatment in luminal breast cancer.
The physiological behaviour of the ERα involves binding of the receptor to its natural ligand estradiol, after which the receptor associates to the chromatin, recruits its coregulators and alters the transcriptional activity of responsive genes, leading to increased cell proliferation and tumour growth. In endocrine treatment of breast cancer, the activation of the receptor can be inhibited through multiple ways, each of which resulting in an inhibition of ERα-driven cell proliferation. Even though these endocrine intervention therapies greatly increase the disease-free survival and overall survival of breast cancer patients , resistance to treatment is commonly observed. This resistance can be mediated through a multitude of different mechanisms, including differential expression of kinases [4, 5], coregulators [3, 6] and transmembrane receptors [7, 8].
Importantly, it is becoming apparent that mechanisms of resistance may result from intrinsic effects on ERα/chromatin interactions, as shown both in cell lines [9, 10] and in tumour samples . Since the ERα/chromatin interactome directly determines the estradiol-mediated effects on gene expression , distinct binding patterns may be a key-defining factor for deviating gene expression patterns as well . A major hurdle to assay such transcription factor/chromatin interactions, is the frequently very limited starting material due to the minimal amount of tissue that is obtained during diagnostic work-up.
Recently approaches have been described to map histone modifications in a limited amount of starting material , but no such protocol has been described for transcription factors. Similarly, new methods for amplification of limited material have been described , but these assume that sufficient DNA was recovered from the ChIP to enable amplification in the first place. One of the major technical limitations has been enriching sufficient DNA (via ChIP) from limited starting material. This problem has been addressed using ‘carrier chromatin’ , a strategy that is not amenable to global sequencing approaches. Here, we describe a procedure to enhance transcription factor ChIP-seq by incorporating a modified carrier method to enable ChIP-seq from limited starting material. We apply this method to map the ERα/chromatin interactome in as few as 10,000 tissue culture cells and apply this technology to map the ERα chromatin binding landscape in core needle biopsy samples from breast cancer patients.
Results and discussion
Carrier selection for small cell number ChIPs
Even though the increased enrichment at the pS2 promoter was indicative of increased ChIP efficiency, this cannot be generalized towards all ERα/chromatin interactions. Therefore, the ERα ChIP-samples were further processed and prepared for sequencing. Only the two extreme conditions (no carrier versus mRNA/Histones carrier) were sequenced and compared to a publicly available dataset of proliferating MCF7 cells, performed on ~20 million cells . The carrier mRNA and histone were degraded prior to amplification of enriched ChIP DNA. An example of difference in binding capacity in the presence and absence of carrier mRNA/histones is shown in Figure 1B and all binding events are illustrated as a heatmap in Figure 1C. Importantly, we ranked all ERα binding peaks observed in the sample with ~20 million cells and found that ~60% of these were recapitulated when using the mRNA and histones as carrier from 10,000 cells. This is a substantial improvement compared to ERα ChIP-seq without carrier, where only signal was detected at ~20% of peaks. Heatmap visualisations and quantifications indicated that peak calling on 10,000 cells with RNA/histones carrier provided comparable signal intensity distributions as the saturated ChIP. This in contrast to the 10,000 cells analysis without carriers (Additional file 2: Figure S2) where mainly background was found. Grouping the ChIP distributions on the basis of signal intensities showed that peaks with high-intensity signal under the RNA/histones carrier ChIP conditions, consistently provided high signal distributions for the saturated ChIP, and vice versa (Additional file 3: Figure S3). No subsets of false-positive and false-negative regions were observed (Additional file 3: Figure S3). Interestingly, the dynamic range of signal intensity as detected in the saturated ChIP was dampened under the carrier conditions (Additional file 3: Figure S3A). Cumulatively, these data show that the addition of carrier chromatin facilitates transcription factor ChIP-seq from 10,000 cells, without negatively effecting the amplification or global fidelity of the process.
Carrier-ChIP from core needle biopsies of human breast tumours
ERα binding events are rarely found at promoters (5%) and the majority of binding sites is normally enriched at introns and distal intergenic regions . We analysed the binding site localizations and found that the mRNA/histones assisted ChIP-seq showed comparable genomic distributions compared to cell line data (Figure 3B). Motif analyses were performed for the mRNA/histone assisted ChIP condition (Figure 3C). As expected, ER DNA binding motifs were the most strongly enriched motifs.
Since ERα is the driving transcription factor in luminal breast cancer, we interrogated the genes proximal to the ERα binding events, and determined whether these were previously identified as basal or luminal gene signatures in independent cohorts of breast cancer patients  (Figure 3D). For the glycogen carrier condition, none of these hallmark genes were found to have a proximal ERα binding event for any of the tumour samples tested (data not shown). For the mRNA/histone-assisted condition, an enrichment of luminal genes was found over basal-signature genes (Figure 3D). Expression of these signature genes with proximal ERα binding events was found to correlate with a favourable outcome after endocrine treatment in a meta-analysis of breast cancer patients with ERα-positive disease  (Figure 3E), indicating that our carrier ChIP-seq protocol enabled us to identify direct ERα target genes with clinical implications from core needle biopsy samples.
The genome-wide assessment of ERα/chromatin interactions has been performed in tissue culture cell lines  and primary breast cancer tissue . Even though ChIP-seq is a very powerful technology to delineate transcription factor-chromatin interactions and to identify specific locations of histone marks on a genome-wide level, the technology does have its intrinsic limitations. The availability of a specific antibody and applicable for immunoprecipitation from formaldehyde-fixed material are major limitations. In addition, there are potential issues related to epitope masking. Another major limitation of ChIP-seq analyses is the need for a certain amount of starting material. Reports have described ChIP-seq analyses on histone marks on limited starting material , but these represent relatively stable covalent interactions, as compared to the dynamic and unstable interactions between transcription factors and chromatin. In directed ChIP analyses, the benefit of carrier during the immunoprecipitation procedure has proven to be substantial . In that method, non-relevant chromatin (e.g. drosophila DNA), is added as a carrier during the procedure, in order to enhance the enrichment from the limited but relevant cellular material. While this approach can be applied for qPCR-ChIP, it is not readily applicable for ChIP-seq, which will sequence all DNA, regardless of the species, resulting in substantial dilution of the relevant enriched DNA. How the addition of the carriers exactly enhances the ChIP procedure on a molecular level still needs to be resolved. Most likely, the bulky material of the carrier helps to retain the small amounts of relevant chromatin throughout the isolation procedure. In addition, it might prove an additional benefit of blocking nonspecific antibody/protein interactions, thereby increasing substrate specificity of the antibody in the reaction to reduce background signal.
Here, we introduce a novel mode of carrier-assisted Chromatin Immunoprecipitation that can be applied to transcription factor mapping by ChIP-seq. The enrichment from limited cellular material is enhanced by the addition of mRNA and recombinant histone protein, which function as a carrier. Importantly, both components of the carrier can subsequently be removed and provide an approach that is cost effective, easily applied and applicable in next-generation sequencing reactions.
We applied this approach to map ERα in 10,000 asynchronous MCF7 cells, which results in an improved DNA recovery, when compared to no carrier. Although this recapitulated ~60% of the signal observed when performing ERα ChIP-seq in saturating numbers of cells, these possessed all the expected characteristics of robust ERα binding events. The dynamic range of the ChIP signal intensity was dampened under carrier conditions, possibly implying that weaker ChIP enrichment regions are enhanced in our method. Furthermore, as a proof-of-principle, we applied our carrier-assisted ChIP protocol to limited material derived from core needle biopsy samples, isolated from breast cancer patients. Again, the presence of the mRNA/histone carrier resulted in increased ERα ChIP-seq signal, allowing for discovery of clinically and biologically relevant ERα regulatory domains from these samples. Several approaches have been recently described [15, 20] that enhance amplification from limited DNA. Although these are powerful methods, they assume that there is sufficient DNA enriched by ChIP, but currently, methods do not exist that enable the ChIP enrichment from limited material for integration with these novel amplification approaches. We believe that our modified carrier ChIP-method improves the ability for purifying transcription factor-associated DNA from limited material, such as tumour samples. The potential exists for coupling our modified ChIP approach with the novel amplification methods for greater ChIP-seq fidelity and signal from primary material, potentially permitting high-resolution transcription factor mapping even from clinical samples.
ChIP was performed as described before , but with a number of modifications. During the Immunoprecipitation, 100 μg/ml glycogen (Invitrogen) or 20 μg/ml recombinant histone 2B (M2505S; New England Biolabs) and 1 μg/ml human mRNA (Invitrogen) was added as carriers for the ChIP. During the entire procedure, non-sticky eppendorf tubes were used (13-698-794; Fisher Scientific). 5 μg antibody used was raised against ERα (SC-543; Santa Cruz), conjugated to 50 μl protein A Dynabeads (Invitrogen). After immunoprecipitation, ten RIPA washes were performed (samples transferred to a fresh eppendorf during the last wash), followed by one TBS wash and reverse crosslinking, as described .
ChIP DNA was amplified as described , using 18 PCR amplification cycles. Sequences generated by the Illumina GAIIx genome analyzer (using 36-bp reads) were aligned against NCBI Build 36.3 of the human genome using MAQ (http://maq.sourceforge.net/) with default parameters. Peaks were called using model-based analysis for ChIP-Seq (MACS) run using default parameters. Read counts and duplication rates for all cell line and tumor samples are shown in Additional file 5: Table S1. The number of clonal reads in the carrier sample was higher and may be a result of increased duplication reads of remaining non-ligated adapters.
Motif analysis, heatmaps and genomic distributions of binding events
ChIP-seq data snapshots were generated using the Integrative Genome Viewer IGV 2.1 (http://www.broadinstitute.org/igv/). Motif analyses were performed through the Cistrome (cistrome.org), applying the SeqPos motif tool . The genomic distributions of binding sites were analysed using the cis-regulatory element annotation system (CEAS) . The genes closest to the binding site on both strands were analysed. If the binding region is within a gene, CEAS software indicates whether it is in a 5′UTR, a 3′UTR, a coding exon, or an intron. Promoter is defined as 3 kb upstream from RefSeq 5′ start. If a binding site is >3 kb away from the RefSeq transcription start site, it is considered distal intergenic.
Tumour tissue handling
Tissue for ChIP-seq analyses was from four newly diagnosed breast cancer patients who were treated at the Netherlands Cancer Institute. Three specimens were collected under a prospective pre-operative endocrine therapy study approved by the institutional ethics committee. These patients provided written informed consent. Also one sample, taken from anonymous left over material, which would be discarded otherwise, was used. Since we are using anonymous, coded leftover material which is not traced back to the patients and does therefore not interfere with care and/ or prognosis, no ethical approval is required according to Dutch legislation [the Medical Research Involving Human Subjects Act; http://www.ccmo-online.nl/main.asp?pid=10&sid=30&ssid=51) and our institutional medical ethical review board. This conforms as described . Three biopsies (14-gauge, 15 μm thick) were taken from each tumour, two of which where frozen with liquid nitrogen and one was formalin-fixed for direct clinical assessment. Tumour cell percentage was determined using Hematoxylin and Eosin (H&E), and tissue was only used for ChIP-seq analysis when tumour cell percentage was higher then 50%. All samples stained positive for Estrogen Receptor. Further information on the tumor samples and patient characteristics can be found in Additional file 4: Table S2.
Gene expression analyses and survival data
Genes with an ERα binding site within 20 kb from the transcription start site were interrogated. Estradiol-responsive genes were used from a publically available database . Luminal and basal genes were previously identified . For survival analyses, a publically available meta-analysis was used. ERα positive-tumours were selected from breast cancer patients who received adjuvant endocrine treatment .
All genomic data is deposited at ArrayExpress under the accession number E-MTAB-1534.
We thank James Hadfield for Solexa sequencing, Gordon Brown, Rory Stark and Suraj Menon for bioinformatic help. We thank Sara Dietz for help on the initial development of the protocol.
We acknowledge the support of The University of Cambridge, Cancer Research UK, the Dutch Cancer Society, A Sisters Hope and The Netherlands Cancer Institute. J.S.C. is supported by an ERC Starting Grant (242664) and W.Z. is supported by a KWF Dutch Cancer Society fellowship.
- Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. International journal of cancer Journal international du cancer. 2010, 127 (12): 2893-2917. 10.1002/ijc.25516.View ArticlePubMedGoogle Scholar
- Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA: Molecular portraits of human breast tumours. Nature. 2000, 406 (6797): 747-752. 10.1038/35021093.View ArticlePubMedGoogle Scholar
- Osborne CK, Bardou V, Hopp TA, Chamness GC, Hilsenbeck SG, Fuqua SA, Wong J, Allred DC, Clark GM, Schiff R: Role of the estrogen receptor coactivator AIB1 (SRC-3) and HER-2/neu in tamoxifen resistance in breast cancer. J Natl Cancer Inst. 2003, 95 (5): 353-361. 10.1093/jnci/95.5.353.View ArticlePubMedGoogle Scholar
- Holm C, Rayala S, Jirstrom K, Stal O, Kumar R, Landberg G: Association between Pak1 expression and subcellular localization and tamoxifen resistance in breast cancer patients. J Natl Cancer Inst. 2006, 98 (10): 671-680. 10.1093/jnci/djj185.View ArticlePubMedGoogle Scholar
- Michalides R, Griekspoor A, Balkenende A, Verwoerd D, Janssen L, Jalink K, Floore A, Velds A, van’t Veer L, Neefjes J: Tamoxifen resistance by a conformational arrest of the estrogen receptor alpha after PKA activation in breast cancer. Cancer Cell. 2004, 5 (6): 597-605. 10.1016/j.ccr.2004.05.016.View ArticlePubMedGoogle Scholar
- Redmond AM, Bane FT, Stafford AT, McIlroy M, Dillon MF, Crotty TB, Hill AD, Young LS: Coassociation of estrogen receptor and p160 proteins predicts resistance to endocrine treatment; SRC-1 is an independent predictor of breast cancer recurrence. Clinical cancer research: an official journal of the American Association for Cancer Research. 2009, 15 (6): 2098-2106. 10.1158/1078-0432.CCR-08-1649.View ArticleGoogle Scholar
- Shou J, Massarweh S, Osborne CK, Wakeling AE, Ali S, Weiss H, Schiff R: Mechanisms of tamoxifen resistance: increased estrogen receptor-HER2/neu cross-talk in ER/HER2-positive breast cancer. J Natl Cancer Inst. 2004, 96 (12): 926-935. 10.1093/jnci/djh166.View ArticlePubMedGoogle Scholar
- Molina MA, Codony-Servat J, Albanell J, Rojo F, Arribas J, Baselga J: Trastuzumab (herceptin), a humanized anti-Her2 receptor monoclonal antibody, inhibits basal and activated Her2 ectodomain cleavage in breast cancer cells. Cancer Res. 2001, 61 (12): 4744-4749.PubMedGoogle Scholar
- Lupien M, Meyer CA, Bailey ST, Eeckhoute J, Cook J, Westerling T, Zhang X, Carroll JS, Rhodes DR, Liu XS: Growth factor stimulation induces a distinct ER(alpha) cistrome underlying breast cancer endocrine resistance. Genes Dev. 2010, 24 (19): 2219-2227. 10.1101/gad.1944810.PubMed CentralView ArticlePubMedGoogle Scholar
- Hurtado A, Holmes KA, Ross-Innes CS, Schmidt D, Carroll JS: FOXA1 is a key determinant of estrogen receptor function and endocrine response. Nat Genet. 2011, 43 (1): 27-33. 10.1038/ng.730.PubMed CentralView ArticlePubMedGoogle Scholar
- Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, Brown GD, Gojis O, Ellis IO, Green AR: Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012, 481 (7381): 389-393.PubMed CentralPubMedGoogle Scholar
- Fullwood MJ, Liu MH, Pan YF, Liu J, Xu H, Mohamed YB, Orlov YL, Velkov S, Ho A, Mei PH: An oestrogen-receptor-alpha-bound human chromatin interactome. Nature. 2009, 462 (7269): 58-64. 10.1038/nature08497.PubMed CentralView ArticlePubMedGoogle Scholar
- Zwart W, Theodorou V, Kok M, Canisius S, Linn S, Carroll JS: Oestrogen receptor-co-factor-chromatin specificity in the transcriptional regulation of breast cancer. EMBO J. 2011, 30 (23): 4764-4776. 10.1038/emboj.2011.368.PubMed CentralView ArticlePubMedGoogle Scholar
- Adli M, Zhu J, Bernstein BE: Genome-wide chromatin maps derived from limited numbers of hematopoietic progenitors. Nat Methods. 2010, 7 (8): 615-618. 10.1038/nmeth.1478.PubMed CentralView ArticlePubMedGoogle Scholar
- Shankaranarayanan P, Mendoza-Parra MA, Walia M, Wang L, Li N, Trindade LM, Gronemeyer H: Single-tube linear DNA amplification (LinDA) for robust ChIP-seq. Nat Methods. 2011, 8 (7): 565-567. 10.1038/nmeth.1626.View ArticlePubMedGoogle Scholar
- O’Neill LP, VerMilyea MD, Turner BM: Epigenetic characterization of the early embryo with a chromatin immunoprecipitation protocol applicable to small cell populations. Nat Genet. 2006, 38 (7): 835-841. 10.1038/ng1820.View ArticlePubMedGoogle Scholar
- Schmidt D, Wilson MD, Spyrou C, Brown GD, Hadfield J, Odom DT: ChIP-seq: using high-throughput sequencing to discover protein-DNA interactions. Methods. 2009, 48 (3): 240-248. 10.1016/j.ymeth.2009.03.001.PubMed CentralView ArticlePubMedGoogle Scholar
- Carroll JS, Meyer CA, Song J, Li W, Geistlinger TR, Eeckhoute J, Brodsky AS, Keeton EK, Fertuck KC, Hall GF: Genome-wide analysis of estrogen receptor binding sites. Nat Genet. 2006, 38 (11): 1289-1297. 10.1038/ng1901.View ArticlePubMedGoogle Scholar
- Gyorffy B, Lanczky A, Eklund AC, Denkert C, Budczies J, Li Q, Szallasi Z: An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010, 123 (3): 725-731. 10.1007/s10549-009-0674-9.View ArticlePubMedGoogle Scholar
- Adli M, Bernstein BE: Whole-genome chromatin profiling from limited numbers of cells using nano-ChIP-seq. Nat Protoc. 2011, 6 (10): 1656-1668. 10.1038/nprot.2011.402.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9 (9): R137-10.1186/gb-2008-9-9-r137.PubMed CentralView ArticlePubMedGoogle Scholar
- He HH, Meyer CA, Shin H, Bailey ST, Wei G, Wang Q, Zhang Y, Xu K, Ni M, Lupien M: Nucleosome dynamics define transcriptional enhancers. Nat Genet. 2010, 42 (4): 343-347. 10.1038/ng.545.PubMed CentralView ArticlePubMedGoogle Scholar
- Ji X, Li W, Song J, Wei L, Liu XS: CEAS: cis-regulatory element annotation system. Nucleic Acids Res. 2006, 34 (Web Server issue): W551-W554.PubMed CentralView ArticlePubMedGoogle Scholar
- van Diest PJ: No consent should be needed for using leftover body material for scientific purposes. For. BMJ. 2002, 325 (7365): 648-651. 10.1136/bmj.325.7365.648.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.