- Open Access
Aberrant CBFA2T3B gene promoter methylation in breast tumors
Molecular Cancervolume 3, Article number: 22 (2004)
The CBFA2T3 locus located on the human chromosome region 16q24.3 is frequently deleted in breast tumors. CBFA2T3 gene expression levels are aberrant in breast tumor cell lines and the CBFA2T3B isoform is a potential tumor suppressor gene. In the absence of identified mutations to further support a role for this gene in tumorigenesis, we explored whether the CBFA2T3B promoter region is aberrantly methylated and whether this correlates with expression.
Aberrant hypo and hypermethylation of the CBFA2T3B promoter was detected in breast tumor cell lines and primary breast tumor samples relative to methylation index interquartile ranges in normal breast counterpart and normal whole blood samples. A statistically significant inverse correlation between aberrant CBFA2T3B promoter methylation and gene expression was established.
CBFA2T3B is a potential breast tumor suppressor gene affected by aberrant promoter methylation and gene expression. The methylation levels were quantitated using a second-round real-time methylation-specific PCR assay. The detection of both hypo and hypermethylation is a technicality regarding the methylation methodology.
Allelic loss of heterozygosity (LOH) of the human chromosome 16q in several sporadic cancer types, including breast, prostate and ovary cancers, suggests this chromosome arm harbors tumor suppressive loci [1–3]. The most frequent region of allelic loss occurs within 3-megabases (Mb) at 16q24.3 between the marker D16S498 and the telomere [1, 4]. Sequencing of the 3-Mb region has identified approximately 100 genes . Eight of these have been excluded as potential tumor suppressors for breast cancer based on mutation analysis in tumor DNA . Recently, a messenger RNA (mRNA) expression survey was completed within 2.4-Mb of this region examining the expression profiles of over 75 genes in a panel of breast tumor cell lines . It was found that only three genes exhibited significantly aberrant expression profiles. These genes were highly expressed in some cell lines and lowly expressed in others. This led to the hypothesis that this aberrant expression may reflect a role for these genes as tumor suppressors in determining cancer phenotype and behavior. One of these genes was the core-binding factor, alpha subunit 2, translocation to 3; termed CBFA2T3.
CBFA2T3 encodes for a protein that belongs to the eight-twenty-one (ETO) family, which also includes the genes CBFA2T1, CBFA2T2 in mammalian cells and nervy in Drosophila . The mammalian members of this protein family are involved in therapy-related chromosomal translocations causing acute myeloid leukemia . The CBFA2T3 gene encodes two alternative transcripts, CBFA2T3A (NM_005187) and CBFA2T3B (NM_175931) (Figure 1A). CBFA2T3A reportedly functions as a nuclear transcriptional co-repressor via its interaction with histone deacetylase (HDAC) complexes . CBFA2T3B functions as a kinase anchorage protein in T lymphocytes and may play a role in inflammatory response . Recently, it was demonstrated that CBFA2T3B functions as a transcriptional repressor and exhibits in vitro characteristics consistent with tumor suppressor activity . This gene was found to be lowly expressed in a number of breast tumor cell lines and upon re-introduction it reduced their growth on plastic and in soft agar.
Mutation analysis in breast tumor cell lines and primary breast tumor samples failed to identify any CBFA2T3 sequence aberrations . It was recognized that aberrant promoter methylation might be the mechanism responsible for the altered expression of CBFA2T3 in breast tumors. The expression of several tumor suppressor genes has been found inactivated or reduced in tumors in association with promoter hypermethylation . Promoter hypermethylation can occur in conjunction with allelic loss and or mutation and is regarded as an alternative form of 'knockout' in biallelic inactivation. Accumulating evidence now suggests that promoter hypermethylation may affect genes that reside within regions of frequent allelic loss more often than mutation . Alternatively, several oncogenes have been found to be over-expressed in tumors in association with promoter hypomethylation . As such, aberrant promoter methylation is considered a fundamental process in developing cancers and has recently received considerable interest as a rapid non-invasive molecular screening tool for the early detection of tumor cells in a range of bodily fluids and biopsy specimens .
In this study, the methylation status of a high-density CpG island promoter region located five prime to the exon 1B sequence of the CBFA2T3B transcript is described (Figure 1A). We explored whether this region is aberrantly methylated in breast tumor cell lines and primary breast tumor samples and whether a correlation exists between methylation and gene expression. Both aberrant hypo and hypermethylation levels were detected in breast tumors in correlation with elevated and reduced expression. The phenomenon of hypo and hypermethylation relates to the amount of DNA used in the methylation methodology as discussed.
Isoform-specific analysis of CBFA2T3 gene expression levels
CBFA2T3 encodes two alternative transcripts, CBFA2T3A and CBFA2T3B (Figure 1A). It was recently demonstrated that CBFA2T3 expression levels are aberrant in breast tumor cell lines. In this study, the total expression levels were assayed using real-time RT-PCR and primers that span exons 4–5 of the CBFA2T3 gene. The tumor suppressor activity previously shown for CBFA2T3B, led us to examine the expression profile of this transcript in breast tumor cell lines using a TaqMan probe specific for exon 1B. CBFA2T3B displayed aberrant expression similar to the total (Figure 2A). This aberrant expression was found to be endogenously low in all samples examined. The expression of this gene was so low that it could not be reliably detected using Northern Blots or RNase protection (data not shown). In addition, expression levels of the CBFA2T3A transcript were assayed in breast tumor cell lines using a TaqMan probe specific for exon 1A. CBFA2T3A expressed at lower levels than CBFA2T3B, but due to the identification of several complex splice variants between exon 1A and exons 1B, 2, 3 and 4, further analysis is required (data not shown). An example of the raw data real-time RT-PCR expression analysis of the CBFA2T3 transcripts is shown (Additional file 1).
Qualitative analysis of CBFA2T3B promoter methylation levels
Based on the idea that this altered expression might develop via aberrant methylation, the promoter activity and methylation status of a high-density CpG island located five prime to the exon 1B sequence of CBFA2T3B was examined (Figures 1A and 1B). The promoter activity was assayed using chloramphenicol acetyl-transferase (CAT) ELISA and confirmed that a 1-kb region spanning the island was capable of promoting a 30-fold increase in CAT expression (Figure 2C). To determine if this promoter region is aberrantly methylated in breast tumors, 24 breast tumor cell lines, 20 primary breast tumors, 20 normal breast counterparts and 24 normal whole blood samples were screened using methylation-specific PCR (MSP). Four separate 100–200 bp regions spanning the promoter were amplified using primers specific for either unmethylated or methylated cytosines at the CpG sites shown (Figure 1B). The full results of this analysis are summarized in Additional file 2. In general, a low 'basal' level of methylation was detected in all samples at the various regions examined. An example of this basal methylation is shown for the normal blood samples at region two (Figure 3A). Unlike the bloods, several breast tumor cell lines, primary breast tumors and normal breast counterpart samples displayed complex high to low methylation levels. An example of this complex methylation is shown for the primary breast tumors and their normal counterparts at region four (Figure 3B). This analysis revealed that only few cell lines displayed clear hypermethylation (e.g. MDA-MB-231) or hypomethylation (e.g. BT-483) in association with reduced and elevated expression. To support this association, it was found that treatment of MDA-MB-231 cells with the demethylating agent 5-aza-2'-deoxycytidine (5-Aza-dC) was capable of increasing both CBFA2T3 and CBFA2T3B expression levels by > 100-fold relative to controls (Figure 2B). Overall, however, it was difficult to comprehend this complex promoter methylation and an association with expression thus suggesting that further analysis was required.
Bisulfite sequence analysis of CBFA2T3B promoter methylation patterns
To understand this complex promoter methylation, sodium bisulfite sequencing of the hypermethylated breast tumor samples was used to examine the pattern and frequency of methylation with this region. This analysis revealed that specific cytosines appeared more susceptible to methylation compared to others (Figure 4). A characteristic sinusoidal pattern of approximately six high to low frequency methylation levels every 40–150 bp was detected in cell lines and tumor samples (Figures 1C and 4). Based on this pattern, primers were designed specific for the real-time MSP quantitation of high-frequency cytosine methylation levels residing within a consensus Specificity protein (Sp1) binding site located approximately minus 450 bp from the transcriptional start of exon 1B (Figures 1B and 1C). As Sp1 proteins are commonly known to regulate gene transcription, it was considered that variable methylation at this site may be reflective of elevated or reduced expression and thus suitable for resolving and correlating the complex promoter methylation levels.
Quantitative analysis of CBFA2T3B promoter methylation levels
To assay CBFA2T3B promoter methylation levels at this Sp1 site, a bisulfite sequencing amplicon spanning this region was initially PCR amplified and column purified from 24 breast tumor cell lines, 55 primary breast tumors, 22 normal breast counterparts and 46 normal whole blood samples. Second-round real-time MSP was performed on the amplicons using internal forward primers to detect for either unmethylated or methylated cytosines at the Sp1 site. Standard curve dilutions of previously prepared internally primed clones representative of either unmethylated or methylated sequence were used to extrapolate methylation levels and normalize for differences in amplification efficiencies. The methylated cytosines were expressed as a fractional ratio of unmethylated cytosines to determine the methylation indices [mi = m/(m + u)]. On plotting the indices, a clear difference between the tumor, normal groups and complex promoter methylation levels was revealed (Figure 5). The normal blood samples maintained a specific basal methylation level and were similar to normal breast counterparts with methylation indices ranging from 0.006–0.09 and 0.002–0.08, respectively. Median methylation index levels were 0.02 in normal bloods and 0.01 in normal breast counterparts. In contrast, several breast tumor cell lines and primary breast tumors were highly variable relative to the normal samples with methylation indices ranging from 0.0002–0.8 and 0.0001–0.9, respectively. A Levene's test revealed a statistically significance difference in variances of methylation indices with the tumor groups being more varied than the normals (P = .001). It was predicted that 83–75% of breast tumor cell lines and 78–69% of primary breast tumors displayed aberrant methylation levels outside the methylation index interquartile ranges of the normal blood (0.02–0.04) and normal breast counterpart samples (0.008–0.03), respectively. Half of the aberrations detected in breast tumor cell lines were either hypo or hypermethylated relative to the normal breast counterpart interquartile ranges. Up to 22% of the primary breast tumors were hypomethylated and 47% hypermethylated relative to the normal breast counterpart interquartile ranges. An example of the raw data real-time MSP and melt curve analysis showing the aberrant methylation levels in breast tumor cell lines compared to normal whole blood samples is shown (Additional files 3 and 4).
Correlation of CBFA2T3B promoter methylation and gene expression levels
To correlate CBFA2T3B promoter methylation levels with gene expression, the methylation indices from 24 breast tumor cell lines and 20 primary breast tumor samples were plotted against their expression. A statistically significant inverse correlation between CBFA2T3B promoter methylation and exon 1B specific expression was established (r2 = 0.63; r = -0.8, P = .0002). Based on the possibility that five prime RNA degradation and secondary structures may have affected the exon 1B complementary DNA (cDNA) synthesis, a correlation between methylation and the total expression is alternatively shown (Figure 6). CBFA2T3B promoter hypermethylation and reduced expression inversely correlated with hypomethylation and elevated expression (r2 = 0.72; r = -0.9, P < .0001). At a hypermethylated index of around 0.9, approximately 4 mRNA molecules (mlcls) per 104 cells were detected compared to a hypomethylated index of 0.0001 and 120,000 mRNA mlcls per 104 cells. The number of CBFA2T3B promoter mlcls methylated per cell for each breast tumor cell line was also calculated by multiplying the methylation index values by the number of 16q24.3 DNA mlcls per cell as previously determined by FISH . These values were plotted against expression in aim to improve the original correlation (r2 = 0.77; r = -0.9, P < .0001) (Additional file 5). By plotting these data sets a power regression was derived which could be used to solve unknown x or y values.
In this study, it has been demonstrated that expression of the CBFA2T3B isoform is altered in breast tumors and that this correlates strongly with aberrant CpG island promoter methylation. Moreover, a comprehensive method for the detection, quantitation and correlation of promoter methylation and gene expression levels has been developed. MSP was used in combination with sodium bisulfite sequencing to identify sites within the CBFA2T3B promoter region displaying high-frequency methylation in breast tumors. Second-round real-time MSP was used to quantitate methylation levels at these sites in breast tumor cell lines, primary breast tumors, normal breast counterparts and normal whole blood samples. The CBFA2T3B promoter methylation levels were calculated as methylation indices and the indices from breast tumors were plotted against their expression.
Validated methylation detection using the MSP and real-time MSP methods
Throughout the development of this method, it was recognized that a pre-requisite for valid MSP amplification required that sufficient amounts of bisulfite modified DNA were used in order to detect methylation and avoid stochastic effects when quantitating low target mlcl numbers and parameters relating to particle distribution statistics. The amount of DNA modified and subsequently amplified is an important parameter in terms of actually detecting potential CpG methylation in a given specimen as discussed.
Most MSP studies use 50 ng of bisulfite modified DNA for amplification. As a normal diploid cell contains an average 6.6 ρg of DNA, this equates that approximately 15,000 alleles or mlcls are made available as starting templates for MSP. If for example, 1 mlcl in 1,000 of these is methylated (mi = 0.001) for a given gene at a specific CpG site(s), then approximately 15 mlcls will act as potential starting templates. Routinely, in unmodified DNA specimens a 'purified' 100–300 bp amplicon diluted to this mlcl number should be detected at around 28–30 cycles of real-time PCR under standard primer efficiencies. Initially, this is a late cycle threshold (CT) for the detection of amplification and is prone to stochastic effects. Under MSP conditions, the DNA has been bisulfite modified which introduces numerous variables that further reduce the overall probability of detection. Upon modification multiple sequence permutations at the CpG site(s) of interest may arise (e.g. from 3 CpG sites a total 8 possible C to TpG permutations exist as 23 = 8). If primers are designed to detect a permutation present at say 4 mlcls in the above example of 15 mlcls (i.e. at attogram starting amounts), then under MSP, a 100–300 bp amplicon at this mlcl number will not amplify until after 35 cycles or may not amplify at all. Moreover, because the initial cycles of MSP are asymmetric the starting template mlcl numbers are further reduced. In addition, it is well known that the DNA is substantially degraded following bisulfite modification with studies demonstrating up to 84–96% degradation . If 1 μg of DNA is modified, as routinely reported, and say 80% is degraded, then based on an unlikely recovery rate of 100%, only 20% of the original pool of templates used in the 50 ng will be available for amplification. Thus, in the above example if only 20% of the 15 mlcls and its possible permutations are available then they may not amplify at all. The pool of primer specific methylated templates could in fact be instantly diminished. When considering the levels of bisulfite-mediated degradation in 1 μg of DNA, even alleles methylated at levels as high as 1 in 10 might be difficult to detect or become spurious. For instance, if assuming that a 100 times more template is available (i.e. 1,500 mlcls) yet permutations exist (i.e. 400 mlcls) and the recover rate is only 4%, then as little as 16 mlcls will be available as starting templates for amplification.
Thus, several technical thresholds exist for the detection of methylation using MSP. The use of 1 μg of DNA for modification, 50 ng for amplification, sequence permutations, PCR efficiency and bisulfite-mediated DNA degradation all reduce the overall probability and validity of detection. To warrant potential detection of methylation in a given specimen at specific CpG sites above a threshold of say 1 in 1,000, it was found necessary to use at least 10X coverage of the 1 μg amount or 10 μg of DNA for modification and 400–500 ng for MSP amplification. Taking into account the above example, this should provide that at least 30–40 mlcls are made available as starting templates, which should amplify within 28 cycles. Specimens with methylation levels lower than this, such as 1 in 50,000, will not be reliably detected and require the use of second-round nested MSP amplification. In this case, the first-round will still require sufficient amounts as any initial stochastic effects may result in poor reproducibility for the second-round. Moreover, as demonstrated in this study, bisulfite sequence analysis of the CBFA2T3B promoter methylation patterns in the hypermethylated MSP samples indicated that the potential for variable methylation frequencies do exist. In this case, it was found that even with sufficient amounts of DNA the qualitative MSP was unreliable and presented complex methylation data. The use of MSP primers for sites that are not methylated or methylated at low to high frequencies, in combination with these other technical thresholds, created complex stochastic methylation data merely decipherable by examining several regions. As a result, this concealed the true methylation status in a majority of the samples under investigation. Only those samples with clear hypermethylation (e.g. MDA-MB-231) or hypomethylation (e.g. BT-483) were greatly reproducible and associated with reduced and elevated expression (Figure 2A and Additional file 2).
To overcome these technical thresholds, a second-round real-time MSP assay was developed. A sequencing amplicon displaying high-frequency methylation sites was initially PCR amplified under 10X coverage conditions and column purified from all samples. Second-round real-time MSP was performed on approximately 100 ρg of the amplicon using internal forward primers to detect for either unmethylated or methylated cytosines at the high-frequency sites. Using this approach, methylation was detectable in every sample examined. All samples amplified with 28 cycles and a methylation index was calculated for these. In comparison, when the standard MSP conditions were assayed by real-time using SYBR Green I, it was found that the CT of methylation detection was below 28 cycles for most samples (data not shown). Extremely low methylation levels, such as those in BT-483, were undetectable using this method but were detected using the second-round. In the latter case, up to 1 methylated mlcl in 10,000 (mi = 0.0001) could be reliably detected within 25 cycles of amplification (Figure 5 and Additional file 3). Notably, no wild-type, unmethylated or methylated cross-amplification was detected when using more DNA. If this occurs then there may be a problem with primer design or modification as the conversion should be complete. In fact, with more DNA a bias in amplification for the target sequence should be created.
In addition, several other technicalities were overcome. The use of only a single second-round primer to quantitate the methylation levels of the high-frequency sites limited the number of possible 5-methyl-CpG sequence permutations that could perturb the accuracy of quantitation. This coupled with the quantitation of methylation levels from absolute standard curves enabled the normalization of reaction efficiencies and calculation of absolute methylation indices [mi = m/(m + u)] and absolute methylation ratios (u/m). The methylation ratios are simply an alternative way of expressing the methylation data (Additional file 6). The calculation of these single methylation values based on the quantitative normalization discriminates against biased amplifications and comparisons of unmethylated to methylated 'band' intensities when using MSP. Moreover, the absolute quantitation of indices and ratios offers improvement over the relative methods, such as comparative cycle thresholds (ΔΔCT), as they have biological significance, are less consuming, more accurate, and do not require the dynamic range in amplification efficiencies of target and references to be similar to enable valid quantitation.
The concept of gene promoter hypo and hypermethylation
Thus, by using this method it was possible to detect CBFA2T3B promoter methylation in all samples. This phenomenon may actually be widespread for genes under the control of methylation as in accordance with the replication model of maintenance methylation [18–20]. This model might suggest that the CpG island is not just randomly methylated; the CpG island is 'always' methylated by memory yet propagated at variable levels for cell type-specific expression or at aberrant levels in association with cancers. Remarkably, in this study it was found that the methylated CBFA2T3B CpG island is propagated as a sinusoidal pattern at aberrant levels in both permanent cell lines and recently resected tumor specimens suggesting that a methylation memory does exist. In fact, it was recognized that the bisulfite sequencing detection levels of this sinusoidal pattern complemented the second real-time MSP quantitations. For example, a methylation ratio in SK-BR-3 of 10:1 (mi = 0.1) was concordant with a 1 in 10 or 10% 5-methylcytosine per 10 clone frequency (Figure 4 and Additional file 6). If by comparing the other bisulfite sequencing levels with their methylation ratios this might suggest that up to 30, 60, 90 or 5,000 clones of 4B, 5B, 3T and MCF-7, respectively, would require sequencing in order to detect 1 methylated sinusoidal mlcl.
Based on the idea that promoter methylation may occur ubiquitously, the concept arises; what is hypo and hypermethylation. In this study, it was demonstrated that hypo and hypermethylation are merely a prediction of levels outside a majority or 'interquartile range' of methylation levels found in the normal samples. For instance, of all the methylation levels detected in the cell lines, only 37.5% of these were predicted as either hypo or hypermethylated relative to the interquartile ranges of normal breast methylation. Moreover, the number of predicted hypo and hypermethylated samples are much lower when compared to the 'full range' of normal methylation ratios as shown in Additional file 6. This is considerably lower for the primary breast tumor samples (i.e. only 24% hypermethylated relative to normal breast) and is likely due to the heterogeneity of breast tumors in concealing the true tumor-related methylation levels. Regardless, by using a more sensitive detection method, the phenomenon of hypo and hypermethylation has appeared.
The phenomenon hypomethylation is a reflection of the low methylation levels that can be detected in the 'hypomethylated' samples relative to the methylated and hypermethylated samples. Alternatively, the phenomenon of hypermethylation is readily observed yet indicates that the detection of methylation alone does not simply represent hypermethylation. The CBFA2T3B gene is methylated in all samples and according to MSP in a high-percentage of samples. In effect, several studies which have examined the methylation status of potential tumor suppressors and cancer-related genes often demonstrate methylation in a high-percentage of samples and state hypermethylation based on the detection of methylation. Examples include, ras-effector nore1A (RASSF1A) [21, 22], stratifin (14-3-3σ) , p15 and p16 , O6-methylguanine-DNA methyltransferase (MGMT) , mismatch repair gene (hMLH1)  and hyperplastic colon polyps gene (HPP1) . In these cases, as with CBFA2T3B, it is probable that only a small percentage of this methylation has cancerous significance in terms of 'methylation-induced silencing', which is routinely confirmed by correlation with 'absence' of expression. Notably, several studies use insufficient amounts of RNA and cDNA for expression analysis and state the absence of amplification to represent inactivation rather than a level of reduction, particularly in cases of low endogenous expression. The need for quantitation to classify methylation levels has been recognized for genes such as glutathione S-transferase P1 (GSTP1)  and adenomatous polyposis coli (APC) , although in these cases the use of low DNA may nevertheless affect intra and interassay reproducibility.
In biological terms, several scenarios exist as to how these aberrant methylation levels might develop. In the case of hypomethylation, a correlation could be made with severely duplicated chromosome copies of 16q24.3 (e.g. BT-483, MCF-7). A possible scenario here is that the hypomethylation is apparent because of a duplicated copy number and thus, for example, not a direct cause of reduced DNA methyltransferase activity or over-expression per se. In the case of hypermethylation, a correlation could be made with 16q24.3 LOH (e.g. MDA-MB-231, MDA-MB-468). An emerging scenario here is that the hypomethylation induces 16q LOH to promote aberrant DNA methyltransferase activity and hypermethylation . Accumulating evidence suggests that the hypermethylation itself, or 'aberrant' methylation, may be targeted to constantly methylated CpG islands (i.e. methylation induces methylation), and in addition targeted to transcriptionally inert CpG islands . In this study, it was found that the CBFA2T3B promoter region is in fact constantly methylated at a median methylation index of 0.02 (i.e. 2 mlcls in 100 are methylated). When comparing this median methylation index with a median gene expression index from the breast tumor cell lines, it was calculated that only approximately 20 mRNA mlcls per 98 'active' alleles are transcribed. These calculations are shown in Additional file 5. This suggests that the CBFA2T3B gene is not only constantly methylated but also largely transcriptionally inert with the remaining active alleles possibly trans-factor dependent for expression.
Although the nature of such targeted aberrations are unknown, several studies demonstrating that altered maintenance and or de novo DNA methyltransferase activities can induce tumorigenesis, clearly demonstrates that the controlled methylation of CpG islands is crucial for normal cell development [30–34]. Accordingly, an in silico prediction of several Sp1, homeotic, epidermal and insulin growth factor recognition sites within the CBFA2T3B promoter region may implicate a role for this element in epithelial development. How the CBFA2T3B CpG island is maintained and dispersed at specified levels within a population of cells is unknown but likely relates to the methylation machinery in controlling distribution within a cell type-specific population (i.e. methylation memory and or phenotype). In this study, it was shown that aberrant deviations outside these specified levels occur profoundly in a majority of breast tumors, particularly the pure tumorigenic cell lines, and as such might comprise an element in tumor formation.
Overall, this study lends further support to the idea that CBFA2T3B is aberrantly regulated in breast cancer. Additional clues supporting a role for this gene in tumor suppression may reside within its protein structure. CBFA2T3B contains a characteristic zinc finger myeloid-nervy-DEAF-1 (zf-MYND) domain and nervy homology regions . Several studies demonstrate that these regions function in transcriptional co-repression via interaction with HDAC and or nuclear co-repressor complexes [9, 35–38]. Accumulating evidence suggests the zf-MYND domain, which is also common to developmental proteins RP-8, DEAF-1, suppressin, Blu, BS69, PDCD2 and Bop, may interact with co-repression complexes to regulate cell-cycle transcription during cell type-specific differentiation [39–46]. Abnormal regulation may be central to tumor formation as supported by reports that several zf-MYND-like proteins display tumor suppressive activity.
Recently, much emphasis has been placed on the development of methylation-based tumor biomarkers for early breast cancer detection to predict disease outcome and strategies for therapy . Interesting data indicates that the use of nipple aspirate fluids may provide for a rapid non-invasive source of screening material for methylation biomarker analysis . Similar studies are underway to evaluate if the methylation status of CBFA2T3B presents biomarker utility. Although in this case, because the normal breast fluids and breast tumors themselves are histologically complex tissues containing a variety of cell types, it is recognized that further studies are required to refine the methylation index interquartile ranges. This may involve large-scale comparisons between normal and tumor cells captured using laser microdissection. Ideally, such normal controls would be resected at autopsy or reduction mammoplasty from non-risk category individuals.
Sample collection and nucleic acid isolation
24 breast tumor cell lines were obtained from the American Type Culture Collection and cultured under recommended conditions. 46 normal whole blood samples, 55 primary breast tumor samples with pathologically classified grade III lesions and 22 adjacent normal breast counterparts samples were obtained with clinical research approval from the Flinders Medical Centre, Department of Surgery. Breast tumor cell line, primary breast tumor and normal breast counterpart genomic DNA was isolated using GenElute for Mammalian Tissues (Sigma). Whole blood genomic DNA was isolated using the QIAamp DNA Blood Kit (Qiagen). Breast tumor cell line and primary breast tumor total RNA was isolated using RNAqueous-4PCR (Ambion). Nucleic acid concentrations were determined using RiboGreen (Molecular Probes). Breast tumor cell line, primary breast tumor and normal tissue (Clontech) total RNA extracts were DNase I treated (Ambion).
Real-time reverse transcription-PCR (RT-PCR)
10–20 μg of total RNA was oligo(dT)16 reverse transcribed at 55°C for 2 h using MMLV (Promega) with addition of RNAguard (Promega) and 5% DMSO. CBFA2T3, CYPA, ATP5A and SYK expression levels were assayed by real-time RT-PCR using SYBR Green I (BMA). CBFA2T3 isoform expression levels were assayed using TaqMan probes specific for exons 1A and 1B (GeneWorks). Real-time RT-PCR was performed on a Rotor-Gene 2000 (Corbett Research) using standard 25 μl HotStar Taq conditions (Qiagen) on cDNA equivalent to 100–1,500 ng total RNA. 0.35X final SYBR Green I or 200 nM probe was used for detection. Amplifications were at 95°C for 10 min, 45 cycles at 94°C for 20 s, annealing temperature for 30 s and 72°C for 30 s. Primer sequences and annealing temperatures are shown (Additional file 7). Unknown expression levels were extrapolated from standard curve dilutions of column purified cDNA amplicons. Replicate standard curve assays (n ≥ 2) were used with CT coefficient variations averaging < 15% over six orders within replicates and between dilutions. mRNA mlcl numbers were quantitated from samples at the parameter CT from standard curves with known mlcls/μl calculated from the dilution mass in μg/μl ÷ M.W. of ssRNA transcript (× 6.02 × 1017 mlcls/μmole). mRNA mlcls per cell were calculated at the CT concentration ÷ amount of total RNA (ρg) per reaction multiplied by 4 based on a 4 ρg total RNA per cell estimation. mRNA was shown as raw expression data (n > 4) and or normalized against CYPA or ATP5A. Several other housekeeping and cancer-related gene expression levels were quantitated to ensure the mRNA mlcl per cell estimations were compatible with other methods (data not shown).
Sodium bisulfite modification
10–20 μg of genomic DNA was digested overnight at 37°C with restriction enzymes XbaI, XhoI, HindIII and EcoRI and cleaned using nucleotide purification columns (QIAvac 24, Qiagen). The digested DNA was pooled and bisulfite modified . 10 μg of DNA was diluted in 500 μl of water and denatured with 55 μl of 2 M NaOH for 20 min at 37°C. DNA was mixed with 300 μl of 10 mM hydroquinone (Sigma), 5.2 ml of 3.6 M NaHSO3 (pH 5.0) (Sigma), overlaid with paraffin oil and deaminated in the dark for 16 h at 55°C. DNA was desalted using Qiagen purification columns, eluted in 500 μl water and desulfonated with 55 μl of 3 M NaOH for 20 min at 37°C. DNA was neutralized and precipitated with 800 μl of 10 M ammonium acetate (pH 7.0), 20 μl linear acrylamide and 5 ml cold 100% EtOH. Modified DNA was pelleted, resuspended in 40 μl 1 mM Tris-Cl (pH 8.0) and the concentrate stored at -80°C.
Methylation-specific PCR (MSP) and sodium bisulfite sequencing
Four 100–200 bp regions spanning approximately 1-kb of the CBFA2T3B promoter region were amplified from ≥ 300 ng of bisulfite modified DNA under standard HotStar Taq conditions using primers specific for either unmethylated or methylated cytosines. Hypermethylated breast tumor samples were used for bisulfite sequencing to generate a CBFA2T3B promoter region methylation map. Four 300–400 bp regions spanning the promoter were amplified under standard conditions as described above using primers simultaneous for both unmethylated and methylated cytosines. The amplicons were sub-cloned into pGEM (Promega) and 5 to 10 clones sequenced using BigDye (Applied Biosystems). Primer sequences and annealing temperatures are shown (Additional file 7).
MDA-MB-231 was treated with the demethylating agent 5-Aza-dC (Sigma). Approximately 2.0 × 105 cells per T75 flask were seeded in 10 ml RPMI-1640 supplemented with 10% FCS, 15 mM HEPES, 10 mg/liter PGS and cultured for 48 h at 37°C with 5% CO2. Cells were treated with 50 μm 5-Aza-dC for 120 h and replenished with fresh medium and 5-Aza-dC every 24 h. These concentrations are not inhibitory to cell growth . Concentrations ranging from 1–5 μM 5-Aza-dC had no affect on expression levels (data not shown). Control untreated cells were cultured in parallel and supplemented with DMSO. Replicate T75 flasks for both treatments and controls were performed (n = 4). Total RNA was isolated and analyzed for CBFA2T3 isoform and SYK expression levels using real-time RT-PCR. MDA-MB-231 cells were also treated with the HDAC inhibiting agents trichostatin A (Sigma), sodium butyrate (Sigma) and apicidin (Calbiochem). Only trichostatin A elicited re-expression levels similar to 5-Aza-dC (data not shown).
Real-time methylation-specific PCR (MSP)
A sequencing amplicon was initially tested for bisulfite PCR amplification and cloning bias by scoring percent methylation frequencies of overlapping amplicons. Bisulfite PCR bias was also tested by amplification on proportional mixtures of hypo and hypermethylated bisulfite treated DNA. This amplicon was PCR amplified and column purified from all samples. Second-round real-time MSP was performed on the amplicon using internal forward primers to detect for either unmethylated or methylated cytosines at the Sp1 site displaying high-frequency cytosine methylation. The methylated and unmethylated primer specificities were tested by real-time amplification on serial dilutions of unmethylated and methylated clones, respectively. Second-round amplicons were also sequenced to ensure specificity. Second real-time MSP was performed on a Rotor-Gene 2000 using internal forward primers under standard 25 μl HotStar Taq conditions with approximately 100 ρg of first-round cleaned amplicon and 0.35X SYBR Green I. Amplifications were at 95°C for 10 min, 45 cycles at 94°C for 20 s, annealing temperature for 30 s and 72°C for 30 s. Unknown unmethylated and methylated cytosine levels were extrapolated at the Sp1 site from standard curve dilutions of internally PCR amplified and column purified cloned amplicons originally sequenced and found to be representative of either unmethylated (e.g. BT-483) or methylated (e.g. MDA-MB-231) sequence. Unmethylated and methylated mlcl numbers were quantitated from all samples at the parameter CT from standard curves with known mlcls/μl calculated from the dilution mass in μg/μl ÷ M.W. of dsDNA clone sequence (× 6.02 × 1017 mlcls/μmole). CBFA2T3B promoter methylation levels were expressed as methylation indices [mi = m/(m + u)] and ratios (u/m) [49, 50].
CBFA2T3B promoter methylation index medians and interquartile ranges were determined for each group of tissue samples. Statistical comparisons between groups were performed using the two-way Levene's test for Equality of Variances, i.e. H0: median methylation index variances are the same in each group; HA: median methylation index variances are not the same in each group. Variances in the median methylation indices between each group were considered statistically significant when P ≤ .05. Correlations between CBFA2T3B promoter methylation and gene expression levels were determined by calculating a Spearman's rank coefficient. All statistical analyses were performed using GraphPad Prism Version 4.0.
Cleton-Jansen AM, Callen DF, Seshadri R, Goldup S, Mccallum B, Crawford J, Powell JA, Settasatian C, van Beerendonk H, Moerland EW, Smit VT, Harris WH, Millis R, Morgan NV, Barnes D, Mathew CG, Cornelisse CJ: Loss of heterozygosity mapping at chromosome arm 16q in 712 breast tumors reveals factors that influence delineation of candidate regions. Cancer Res. 2001, 61: 1171-1177.
Suzuki H, Komiya A, Emi M, Kuramochi H, Shiraishi T, Yatani R, Shimazaki J: Three distinct commonly deleted regions of chromosome arm 16q in human primary and metastatic prostate cancers. Genes Chromosomes Cancer. 1996, 17: 225-233. 10.1002/(SICI)1098-2264(199612)17:4<225::AID-GCC4>3.3.CO;2-S
Hansen LL, Jensen LL, Dimitrakakis C, Michalas S, Gilbert F, Barber HR, Overgaard J, Arzimanoglou II: Allelic imbalance in selected chromosomal regions in ovarian cancer. Cancer Genet Cytogenet. 2002, 139: 1-8. 10.1016/S0165-4608(02)00620-9
Whitmore SA, Crawford J, Apostolou S, Eyre H, Baker E, Lower KM, Settasatian C, Goldup S, Seshadri R, Gibson RA, Mathew CG, Cleton-Jansen AM, Savoia A, Pronk JC, Auerbach AD, Doggett NA, Sutherland GR, Callen DF: Construction of a high-resolution physical and transcription map of chromosome 16q24.3: a region of frequent loss of heterozygosity in sporadic breast cancer. Genomics. 1998, 50: 1-8. 10.1006/geno.1998.5316
Crawford J, Ianzano L, Savino M, Whitmore S, Cleton-Jansen AM, Settasatian C, d'apolito M, Seshadri R, Pronk JC, Auerbach AD, Verlander PC, Mathew CG, Tipping AJ, Doggett NA, Zelante L, Callen DF, Savoia A: The PISSLRE gene: structure, exon skipping, and exclusion as tumour suppressor in breast cancer. Genomics. 1999, 56: 90-97. 10.1006/geno.1998.5676
Powell JA, Gardner AE, Bais AJ, Hinze SJ, Baker E, Whitmore S, Crawford J, Kochetkova M, Spendlove HE, Doggett NA, Sutherland GR, Callen DF, Kremmidiotis G: Sequencing, transcript identification and quantitative gene expression profiling in the breast cancer loss of heterozygosity region 16q24.3 reveal three potential tumor-suppressor genes. Genomics. 2002, 80: 303-310.
Davis JN, McGhee L, Meyers S: The ETO (MTG8) gene family. Gene. 2003, 303: 1-10. 10.1016/S0378-1119(02)01172-1
Gamou T, Kitamura E, Hosoda F, Shimizu K, Shinohara K, Hayashi Y, Nagase T, Yokoyama Y, Ohki M: The partner gene of AML1 in t(16;21) myeloid malignancies is a novel member of the MTG8 (ETO) family. Blood. 1998, 91: 4028-4037.
Hoogeveen AT, Rossetti S, Stoyanova V, Schonkeren J, Fenaroli A, Schiaffonati L, van Unen L, Sacchi N: The transcriptional corepressor MTG16a contains a novel nucleolar targeting sequence deranged in t(16;21)-positive myeloid malignancies. Oncogene. 2002, 21: 6703-6712. 10.1038/sj.onc.1205882
Schillace RV, Andrews SF, Liberty GA, Davey MP, Carr DW: Identification and characterization of myeloid translocation gene 16b as a novel a kinase anchoring protein in T lymphocytes. J Immunol. 2002, 168: 1590-1599.
Kochetkova M, McKenzie OL, Bais AJ, Martin JM, Secker GA, Seshadri R, Powell JA, Hinze SJ, Gardner AE, Spendlove HE, O'Callaghan NJ, Cleton-Jansen AM, Cornelisse C, Whitmore SA, Crawford J, Kremmidiotis G, Sutherland GR, Callen DF: CBFA2T3 (MTG16) is a putative breast tumor suppressor gene from the breast cancer loss of heterozygosity region at 16q24.3. Cancer Res. 2002, 62: 4599-4604.
Baylin SB, Herman JG: DNA hypermethylation in tumorigenesis: epigenetics joins genetics. Trends Genet. 2000, 16: 168-174. 10.1016/S0168-9525(99)01971-X
Jones PA, Baylin SB: The fundamental role of epigenetic events in cancer. Nat Rev Genet. 2002, 3: 415-428. 10.1038/nrg962
Counts JL, Goodman JI: Hypomethylation of DNA: a possible epigenetic mechanism involved in tumor promotion. Prog Clin Biol Res. 1995, 391: 81-101.
Laird PW: The power and the promise of DNA methylation markers. Nat Rev Cancer. 2003, 3: 253-266. 10.1038/nrc1045
Callen DF, Crawford J, Derwas C, Cleton-Jansen AM, Cornelisse CJ, Baker E: Defining regions of loss of heterozygosity of 16q in breast cancer cell lines. Cancer Genet Cytogenet. 2002, 133: 76-82. 10.1016/S0165-4608(01)00565-9
Grunau C, Clark SJ, Rosenthal A: Bisulfite genomic sequencing: systematic investigation of critical experimental parameters. Nucleic Acids Res. 2001, 29: E65- 10.1093/nar/29.13.e65
Holliday R, Pugh JE: DNA modification mechanisms and gene activity during development. Science. 1975, 186: 226-232.
Bird A: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16: 6-21. 10.1101/gad.947102
Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet. 2003, 33: 245-254. 10.1038/ng1089
Dammann R, Yang G, Pfeifer GP: Hypermethylation of the CpG island of Ras association domain family 1A (RASSF1A), a putative tumor suppressor gene from the 3p21.3 locus, occurs in a large percentage of human breast cancers. Cancer Res. 2001, 61: 3105-3109.
Lo KW, Kwong J, Hui AB, Chan SY, To KF, Chan AS, Chow LS, Teo PM, Johnson PJ, Huang DP: High frequency of promoter hypermethylation of RASSF1A in nasopharyngeal carcinoma. Cancer Res. 2001, 61: 3877-3881.
Ferguson AT, Evron E, Umbricht CB, Pandita TK, Chan TA, Hermeking H, Marks JR, Lambers AR, Futreal PA, Stampfer MR, Sukumar S: High frequency of hypermethylation at the 14-3-3 sigma locus leads to gene silencing in breast cancer. Proc Natl Acad Sci U S A. 2000, 97: 6049-6054. 10.1073/pnas.100566997
Wong IH, Lo YM, Yeo W, Lau WY, Johnson PJ: Frequent p15 promoter methylation in tumor and peripheral blood from hepatocellular carcinoma patients. Clin Cancer Res. 2000, 6: 3516-3521.
Palmisano WA, Divine KK, Saccomanno G, Gilliland FD, Baylin SB, Herman JG, Belinsky SA: Predicting lung cancer by detecting aberrant promoter methylation in sputum. Cancer Res. 2000, 60: 5954-5958.
Simpkins SB, Bocker T, Swisher EM, Mutch DG, Gersell DJ, Kovatich AJ, Palazzo JP, Fishel R, Goodfellow PJ: MLH1 promoter methylation and gene silencing is the primary cause of microsatellite instability in sporadic endometrial cancers. Hum Mol Genet. 1999, 8: 661-666. 10.1093/hmg/8.4.661
Young J, Biden KG, Simms LA, Huggard P, Karamatic R, Eyre HJ, Sutherland GR, Herath N, Barker M, Anderson GJ, Fitzpatrick DR, Ramm GA, Jass JR, Leggett BA: HPP1: a transmembrane protein-encoding gene commonly methylated in colorectal polyps and cancers. Proc Natl Acad Sci U S A. 2001, 98: 265-270. 10.1073/pnas.011415298
Goessl C, Krause H, Muller M, Heicappell R, Schrader M, Sachsinger J, Miller K: Fluorescent methylation-specific polymerase chain reaction for DNA-based detection of prostate cancer in bodily fluids. Cancer Res. 2000, 60: 5941-5945.
Usadel H, Brabender J, Danenberg KD, Jeronimo C, Harden S, Engles J, Danenberg PV, Yang S, Sidransky D: Quantitative adenomatous polyposis coli promoter methylation analysis in tumor tissue, serum, and plasma DNA of patients with lung cancer. Cancer Res. 2002, 62: 371-375.
Chen RZ, Pettersson U, Beard C, Jackson-Grusby L, Jaenisch R: DNA hypomethylation leads to elevated mutation rates. Nature. 1998, 395: 89-93. 10.1038/25779
Rhee I, Jair KW, Yen RW, Lengauer C, Herman JG, Kinzler KW, Vogelstein B, Baylin SB, Schuebel KE: CpG methylation is maintained in human cancer cells lacking DNMT1. Nature. 2000, 404: 1003-1007. 10.1038/35010000
Rhee I, Bachman KE, Park BH, Jair KW, Yen RW, Schuebel KE, Cui H, Feinberg AP, Lengauer C, Kinzler KW, Baylin SB, Vogelstein B: DNMT1 and DNMT3b cooperate to silence genes in human cancer cells. Nature. 2002, 416: 552-556. 10.1038/416552a
Robert MF, Morin S, Beaulieu N, Gauthier F, Chute IC, Barsalou A, MacLeod AR: DNMT1 is required to maintain CpG methylation and aberrant gene silencing in human cancer cells. Nat Genet. 2003, 33: 61-65. 10.1038/ng1068
Eden A, Gaudet F, Waghmare A, Jaenisch R: Chromosomal instability and tumors promoted by DNA hypomethylation. Science. 2003, 300: 455- 10.1126/science.1083557
Wang J, Hoshino T, Redner RL, Kajigaya S, Liu JM: ETO, fusion partner in t(8;21) acute myeloid leukemia, represses transcription by interaction with the human N-CoR/mSin3/HDAC1 complex. Proc Natl Acad Sci U S A. 1998, 95: 10860-10865. 10.1073/pnas.95.18.10860
Lutterbach B, Westendorf JJ, Linggi B, Patten A, Moniwa M, Davie JR, Huynh KD, Bardwell VJ, Lavinsky RM, Rosenfeld MG, Glass C, Seto E, Hiebert SW: ETO, a target of t(8;21) in acute leukemia, interacts with the N-CoR and mSin3 corepressors. Mol Cell Biol. 1998, 18: 7176-7184.
Gelmetti V, Zhang J, Fanelli M, Minucci S, Pelicci PG, Lazar MA: Aberrant recruitment of the nuclear receptor corepressor-histone deacetylase complex by the acute myeloid leukemia fusion partner ETO. Mol Cell Biol. 1998, 18: 7185-7191.
Amann JM, Nip J, Strom DK, Lutterbach B, Harada H, Lenny N, Downing JR, Meyers S, Hiebert SW: ETO, a target of t(8;21) in acute leukemia, makes distinct contacts with multiple histone deacetylases and binds mSin3A through its oligomerization domain. Mol Cell Biol. 2001, 21: 6470-6483. 10.1128/MCB.21.19.6470-6483.2001
Owens GP, Hahn WE, Cohen JJ: Identification of mRNAs associated with programmed cell death in immature thymocytes. Mol Cell Biol. 1991, 11: 4177-4188.
Pederson JA, LaFollette JW, Gross C, Veraksa A, McGinnis W, Mahaffey JW: Regulation by homeoproteins: a comparison of deformed-responsive elements. Genetics. 2000, 156: 677-686.
LeBoeuf RD, Burns JN, Bost KL, Blalock JE: Isolation, purification, and partial characterization of suppressin, a novel inhibitor of cell proliferation. J Biol Chem. 1990, 265: 158-165.
Agathanggelou A, Dallol A, Zochbauer-Muller S, Morrissey C, Honorio S, Hesson L, Martinsson T, Fong KM, Kuo MJ, Yuen PW, Maher ER, Minna JD, Latif F: Epigenetic inactivation of the candidate 3p21.3 suppressor gene BLU in human cancers. Oncogene. 2003, 22: 1580-1588. 10.1038/sj.onc.1206243
Ansieau S, Leutz A: The conserved Mynd domain of BS69 binds cellular and oncoviral proteins through a common PXLXP motif. J Biol Chem. 2002, 277: 4906-4910.
Masselink H, Bernards R: The adenovirus E1A binding protein BS69 is a corepressor of transcription through recruitment of N-CoR. Oncogene. 2000, 19: 1538-1546. 10.1038/sj.onc.1203421
Scarr RB, Sharp PA: PDCD2 is a negative regulator of HCF-1 (C1). Oncogene. 2002, 21: 5245-5254. 10.1038/sj.onc.1205647
Gottlieb PD, Pierce SA, Sims RJ, Yamagishi H, Weihe EK, Harriss JV, Maika SD, Kuziel WA, King HL, Olson EN, Nakagawa O, Srivastava D: Bop encodes a muscle-restricted protein containing MYND and SET domains and is essential for cardiac differentiation and morphogenesis. Nat Genet. 2002, 31: 25-32.
Krassenstein R, Sauter E, Dulaimi E, Battagli C, Ehya H, Klein-Szanto A, Cairns P: Detection of breast cancer in nipple aspirate fluid by CpG island hypermethylation. Clin Cancer Res. 2004, 10: 28-32.
Li B, Goyal J, Dhar S, Dimri G, Evron E, Sukumar S, Wazer DE, Band V: CpG methylation as a basis for breast tumor-specific loss of NES1/kallikrein 10 expression. Cancer Res. 2001, 61: 8014-8021.
Lo YM, Wong IH, Zhang J, Tein MS, Ng MH, Hjelm NM: Quantitative analysis of aberrant p16 methylation using real-time quantitative methylation-specific polymerase chain reaction. Cancer Res. 1999, 59: 3899-3903.
Pfeifer GP, Steigerwald SD, Hansen RS, Gartler SM, Riggs AD: Polymerase chain reaction-aided genomic sequencing of an X chromosome-linked CpG island: methylation patterns suggest clonal inheritance, CpG site autonomy, and an explanation of activity state stability. Proc Natl Acad Sci U S A. 1990, 87: 8252-8256.
AJB completed the work and manuscript. AEG aided in bisulfite sequencing. OLDM performed promoter analysis. DFC, GRS and GK supervised the work. All authors read and approved the final manuscript.