Aberrant CBFA2T3B gene promoter methylation in breast tumors
© Bais et al; licensee BioMed Central Ltd. 2004
Received: 28 June 2004
Accepted: 10 August 2004
Published: 10 August 2004
Skip to main content
© Bais et al; licensee BioMed Central Ltd. 2004
Received: 28 June 2004
Accepted: 10 August 2004
Published: 10 August 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.
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.
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.
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.
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.
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).
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).
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.
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).
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.
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.