Role of aldo-keto reductases and other doxorubicin pharmacokinetic genes in doxorubicin resistance, DNA binding, and subcellular localization
© Heibein et al.; licensee BioMed Central Ltd. 2012
Received: 19 January 2012
Accepted: 17 August 2012
Published: 31 August 2012
Since proteins involved in chemotherapy drug pharmacokinetics and pharmacodynamics have a strong impact on the uptake, metabolism, and efflux of such drugs, they likely play critical roles in resistance to chemotherapy drugs in cancer patients.
To investigate this hypothesis, we conducted a whole genome microarray study to identify difference in the expression of genes between isogenic doxorubicin-sensitive and doxorubicin-resistant MCF-7 breast tumour cells. We then assessed the degree of over-representation of doxorubicin pharmacokinetic and pharmacodynamic genes in the dataset of doxorubicin resistance genes.
Of 27,958 Entrez genes on the array, 7.4 per cent or 2,063 genes were differentially expressed by ≥ 2-fold between wildtype and doxorubicin-resistant cells. The false discovery rate was set at 0.01 and the minimum p value for significance for any gene within the “hit list” was 0.01. Seventeen and 43 per cent of doxorubicin pharmacokinetic genes were over-represented in the hit list, depending upon whether the gene name was identical or within the same gene family, respectively. The most over-represented genes were within the 1C and 1B families of aldo-keto reductases (AKRs), which convert doxorubicin to doxorubicinol. Other genes convert doxorubicin to other metabolites or affect the influx, efflux, or cytotoxicity of the drug. In further support of the role of AKRs in doxorubicin resistance, we observed that, in comparison to doxorubicin, doxorubincol exhibited dramatically reduced cytotoxicity, reduced DNA-binding activity, and strong localization to extra nuclear lysosomes. Pharmacologic inhibition of the above AKRs in doxorubicin-resistant cells increased cellular doxorubicin levels, restored doxorubicin cytotoxicity and re-established doxorubicin localization to the nucleus. The properties of doxorubicinol were unaffected.
These findings demonstrate the utility of using curated pharmacokinetic and pharmacodynamic knowledge bases to identify highly relevant genes associated with doxorubicin resistance. The induction of one or more of these genes was found to be correlated with changes in the drug’s properties, while inhibiting one specific class of these genes (the AKRs) increased cellular doxorubicin content and restored drug DNA binding, cytotoxicity, and subcellular localization.
KeywordsDoxorubicin Resistance Gene profiling PharmGKB Pharmacokinetics Aldo-keto reductases Cytotoxicity Lysosome DNA-binding Drug localization
Doxorubicin is a DNA-binding, topoisomerase II inhibitor [1, 2], which is among the most effective chemotherapy drugs in cancer treatment [1, 3]. However, intrinsic or acquired resistance to doxorubicin in patient tumours is common, resulting in treatment failure and disease progression. Multiple mechanisms for doxorubicin resistance have been identified in vitro, such as the increased expression of drug transporters [4–7], alterations in doxorubicin metabolism  or localization [9, 10], and defects in the drug’s ability to induce apoptosis . Unfortunately, progress in restoring drug sensitivity for drug-resistant tumours, particularly by inhibiting drug efflux transporters, has been incremental at best [12, 13]. This limited progress demands that a more nuanced approach be taken, including the identification of all proteins that likely affect the pharmacokinetics and pharmacodynamics of doxorubicin.
Genome profiling is a method that can provide data on gene expression and/or allelic variations across biological samples, often using whole genome approaches. This promises to be a great aid to oncologists in identifying and treating drug-resistant tumours. Unfortunately, this task is a difficult one, given the variability associated with patient data sets and the large number of “false positives” inherent in such approaches from by-stander effects. One method to improve the identification of genes relevant to a specific phenomenon such as doxorubicin resistance is to pair knowledge of metabolic or signal transduction pathways to gene expression data . In this study, we use full genome microarray analysis to compare gene expression between MCF-7 cells selected for maximal resistance to doxorubicin (MCF-7DOX2-12 cells) and equivalent cells selected for the same number of passages in the absence of drug (MCF-7CC12 cells). After identifying genes having altered expression in doxorubicin-resistant cells, we then used a well-known, curated pharmacogenomics knowledgebase (PharmGKB) to identify which of these genes play a role in doxorubicin pharmacokinetics or pharmacodynamics, as these were more likely to have a direct effect on doxorubicin efficacy. This combination of full genome microarray analysis identifying genes differentially expressed upon acquisition of doxorubicin resistance with an assessment of over-representation of doxorubicin pharmacokinetic or pharmacokinetic genes in the dataset provided significant insight into new pathways associated with doxorubicin resistance. Moreover, extensive comparisons between the biochemical properties of doxorubicin and one of its metabolites (doxorubicinol) provided us with significant insight into how a simple hydroxylation reaction can strongly affect the biochemical and cellular properties of doxorubicin, including dramatically reduced cytotoxicity, diminished DNA- binding activity, altered cellular accumulation of the drug and altered subcellular localization.
Differentially expressed genes upon acquisition of doxorubicin resistance
Using full genome Agilent microarrays and Partek Genomics Suite, 2063 genes from a total of 27958 Entrez genes on the array (7.4%) were found to be differentially expressed by ≥2-fold between MCF-7CC12 cells MCF-7DOX2-12 cells. The false discovery rate was set at 0.01 and the minimum p value for significance for any gene within the “hit list” was 0.01. The microarray data was deposited in the NCBI Gene Expression Omnibus (GEO) database, accession number GSE27254) in accordance with MIAME standards . Access to the microarray data can be obtained via the following url: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dbezngycywquuhm&acc=GSE27254.
Several AKRs are over-expressed in MCF-7DOX2-12cells
As previously demonstrated using a much smaller microarray platform (1720 genes) , the “1C” family of AKRs was observed to be over-expressed upon acquisition of doxorubicin resistance. Moreover, as shown in Additional file 1: Table S1, a variety of AKR family members were among the most differentially expressed genes upon acquisition of doxorubicin resistance in MCF-7 cells. In these microarray studies, AKR1B1, AKR1B10, AKR1C1, and AKR1C3 all had strongly elevated expression (10.0-, 13.4-, 4.45-, and 4.71-fold, respectively). As stated previously, the product of the AKR family of genes facilitates the conversion of doxorubicin to doxorubicinol . Such a strong overexpression of multiple AKR transcripts in MCF-7DOX2-12 cells suggests that the AKRs may play a major role in doxorubicin resistance.
Doxorubicinol is 1 million-fold less cytotoxic than doxorubicin in MCF-7 cells
5β-cholanic acid restores sensitivity of MCF-7DOX2-12cells to doxorubicin
As illustrated in Figure 4, treatment of MCF-7DOX2-12 cells with both doxorubicin and 5β-cholanic acid, a potent inhibitor of AKR1B10 , AKR1C2, and AKR1C3 activity , almost fully restored doxorubicin sensitivity to that of MCF-7CC12 cells (IC50 of 22.9 ± 9 nM). In contrast, treatment of MCF-7CC12 cells with 5β-cholanic acid and doxorubicin had little effect on doxorubicin sensitivity (IC50 of 14.1 ± 8.0 nM), suggesting insufficient AKR activity in these cells to affect doxorubicin sensitivity. Addition of 5β-cholanic acid had no effect on sensitivity of MCF-7CC12 cells to doxorubicinol (IC50 of 7.5 ± 10.6 mM). However, addition of 5β-cholanic acid to MCF-7DOX2-12 cells did appear to increase their sensitivity to doxorubicinol to a barely detectable range (IC50 of 3.2 ± 1.7 mM), suggesting a possible ability of the inhibitor to affect further metabolism of doxorubicinol in doxorubicin-resistant cells.
Restoration of doxorubicin sensitivity is accompanied by restored nuclear localization in MCF-7DOX2-12cells
We then assessed whether co-treatment of cells with 5β-cholanic acid altered doxorubicin or doxorubicinol localization (Figure 5). Interestingly, 200 μM 5β-cholanic acid was able to completely restore doxorubicin localization to the nucleus of MCF-7DOX2-12 cells, suggesting that the conversion of doxorubicin to doxorubicinol does alter the drug’s ability to reach or bind its target. The same concentration of 5β-cholanic acid, however, had no effect on doxorubicinol localization in MCF-7CC12 and MCF-7DOX2-12 cells.
Doxorubicinol fails to accumulate in MCF-7CC12 and MCF-7DOX2-12cells
Intracellular levels of doxorubicin are significantly altered upon treatment of MCF-7DOX2-12cells with 5β-cholanic acid and/or cyclosporine A
Treatment of MCF-7CC12 cells with 5β-cholanic acid and the pan ABC transporter inhibitor cyclosporine A increased cellular doxorubicin content by 51% and 80%, respectively (Figure 6B). Addition of both agents increased doxorubicin content to almost twice that of untreated cells, but none of the above differences in doxorubicin content were considered statistically significant. In contrast, 5β-cholanic acid or cyclosporine A significantly increased doxorubicin content in MCF-7DOX2-12 cells by 2.8-fold (Figure 6B). Treatment of MCF-7DOX2-12 cells with both 5β-cholanic acid and cyclosporine A increased cellular doxorubicin content to levels 4.4-fold higher than untreated cells (Figure 6B). These differences relative to untreated cells were found to be highly significant, and are likely due to the increased expression of AKRs  and ABC drug transporters known to be over-expressed in MCF-7DOX2-12 cells, including Abcc1 .
Doxorubicinol binds to DNA with lower affinity than doxorubicin
Use of the binomial statistic to interpret the significance of pathways in gene expression data
DNA microarray, high throughput quantitative PCR, and other gene profiling approaches have been highly useful in identifying differences in gene expression between cells or tumours responding to chemotherapy agents and those that do not. Unfortunately, the false discovery rate for such approaches is quite high, largely due to the identification of a large number of “passenger genes” unrelated to drug response. A wide variety of pathway analysis tools exist today, some manually curated, and some created primarily through machine learning. PharmGKB , Ariadne Pathway Studio , Reactome , Ingenuity Pathway Analysis (http://www.ingenuity.com/), GenMAPP , and DAVID , are examples of available tools which can be used to map changes in gene expression to alterations in biochemical pathways. The difficulty with this approach is the sheer size of the data sets, the large number of documented pathways, and the complex statistics required to determine the significance of findings. In this study we elected to use a simple model to examine the biology of doxorubicin resistance, namely looking for “over-representation” of doxorubicin pharmacokinetic and pharmacodynamic genes in datasets of genes having altered expression in doxorubicin resistance.
In order to assess the feasibility of this approach and to survey the broadest number of genes, we used Agilent full genome microarrays containing 27,958 Entrez gene probes, unlike our previous study of only 1720 gene probes . This approach helped to uncover a number of AKRs induced during selection for doxorubicin resistance, including AKR1C1, AKR1B1, AKR1B10, and AKR1C3. Their expression was elevated between 4.5- and 13.4-fold (Additional file 1: Table S1). Given that the probes for these AKRs on the Agilent microarrays were not isoform-specific, we used RTqPCR with isoform-specific primers (Table 4) to determine that, upon selection for doxorubicin resistance, transcripts for AKR1C2, AKR1C3, and AKR1B10 were overexpressed 3.6-, 9.1-, and 10.4-fold, respectively (Figure 3). In addition to the AKRs, other over-represented genes (Figure 2) provide further insight into other proteins that likely contribute to doxorubicin resistance. For example, NQO1 codes for NAD(P)H dehydrogenase quinone 1, which plays a role in converting doxorubicin to doxorubicin deoxyaglycone or to doxorubicin semiquinone (Figure 2). Its 3-fold increase in expression might therefore increase the conversion of doxorubicin to these metabolites as well. Transcripts for the drug efflux pump Abcc1 were also upregulated 8.3-fold, as well as transcripts for other ATP-binding cassette (ABC) transporters such as Abcd3, Abcg2, and Abca1. In addition, a gene (SLC22A15) homologous to the solute carrier protein Slc22a16 (which promotes doxorubicin uptake into cells ) was found to be down regulated by 2.8-fold. The combined changes in the expression of ABC transporters and solute carrier proteins would be expected to reduce doxorubicin accumulation into cells. The gene for catalase (CAT) was found to be upregulated 3.6-fold in MCF-7DOX2-12 cells. Since its gene product helps protect cells from oxidative damage by reactive oxygen species , its elevated expression would protect cells from reactive oxygen species known to be generated by doxorubicin. Genes associated with the cardiotoxicity of doxorubicin (through negative effects on mitochondrial function when converted to doxorubicinol) also have altered expression in breast tumour cells upon selection for doxorubicin resistance, including ACO1, ATPS, CYCS, and ATP2B4 (Figure 2 and Additional file 1: Table S1).
Of the above-described changes in gene expression, the greatest were for the AKRs. Evidence provided in this study supports their substantial role in doxorubicin resistance in tumour cells in vitro, and possibly in the tumours of cancer patients. While many of the changes in gene expression identified in our microarray study likely play a bona fide role in doxorubicin resistance (given their roles in cells), some of the identified genes may not be the “drivers” of drug resistance, but change expression through the altered expression of the driver genes.
Role of the AKRs in resistance to doxorubicin
A role for AKRs in xenobiotic and anthracycline metabolism has already been well established in the literature [21, 35–37]. We also published previously that aldo-keto reductases (AKRs) are overexpressed upon acquisition of anthracycline resistance, that doxorubicin localization to the nucleus is altered in doxorubicin-resistant cells, and that inhibition of AKRs restores doxorubicin sensitivity in doxorubicin-resistant cells . However, the current study significantly extends these observations in many respects. For example, it reveals that the expression of other members of the AKR family is elevated as breast tumour cells acquire resistance to doxorubicin. This would further increase the production of doxorubicinol and its possible conversion to other downstream metabolites. Moreover, our study provides a detailed comparison between doxorubicin and doxorubicinol in terms of their cytotoxicity, subcellular localization, and DNA binding activity.
Interestingly, despite having identical fluorescence capacities, cellular levels of doxorubicinol in both MCF-7CC12 and MCF-7DOX2-12 cells was considerably lower than that of doxorubicin (as measured by cellular fluorescence intensity after drug administration and washing away free drug not taken up by cells). This decreased doxorubicinol uptake may be because hydroxylated doxorubicin is more polar and less able to traverse the hydrophobic plasma membrane. Moreover, even if the confocal microscope settings are modified to allow greater sensitivity to detect cellular doxorubicinol, doxorubicinol was found not to be localized to the nucleus in both MCF-7CC12 and MCF-7DOX2-12 cells. This indicates that the differential localization of doxorubicin between MCF-7CC12 and MCF-7DOX2-12 cells may be due to the strongly elevated conversion of doxorubicin to doxorubicinol (or other fluorescent metabolites) in MCF-7DOX2-12 cells. This may be why “doxorubicin” had an altered location in anthracycline-resistant cells in our previous study. The fluorescence observed in lysosomes may be that of doxorubicin, but also of doxorubicinol and other fluorescent doxorubicin metabolites. Consistent with this view, and not reported in our previous study, the administration of the AKR inhibitor 5β-cholanic acid significantly restored “doxorubicin” localization to the nucleus. More likely the inhibitor prevented doxorubicin conversion to doxorubicinol, permitting more doxorubicin to be retained within the nucleus.
What could account for the decreased localization of doxorubicin to the nucleus? We report in the current study that doxorubicinol has significantly lower ability to bind to DNA than doxorubicin (altered Bmax and Kapp). The conversion of doxorubicin to doxorubicinol by AKRs would result in reduced binding to DNA and hence less ability of the drug to remain associated with the nucleus. In our previous study, we did not differentiate between the cellular localization of doxorubicin and doxorubicinol.
One surprising finding in our study was the lack of detection of significant doxorubicinol in MCF-7DOX2-12 cells (Figure 6). This was despite the elevated expression of a number of AKRs in the cell line (Figure 3), which would be expected to covert doxorubicin to doxorubicinol. And yet, the addition of 5β-cholanic acid with doxorubicin increased the cellular content of doxorubicin (Figure 6), supporting the observation that 5β-cholanic acid is able to block the conversion of doxorubicin to doxorubicinol. What may account for the discrepancy in these points of view? One possibility is that 5β-cholanic acid blocks the efflux of doxorubicin by drug transporters (possibly Abcc1), thereby increasing the retention of doxorubicin in cells. One argument against this hypothesis is that both 5β-cholanic acid and cyclosporine A increased cellular doxorubicin content (Figure 6), the latter being a known inhibitor of Abcc1 function . The combination of both agents increased cellular doxorubicin content further, suggesting that they were acting by distinct mechanisms. Moreover, unlike 5β-cholanic acid (Figure 4), addition of cyclosporine A had no effect on the cytotoxicity of doxorubicin in MCF-7DOX2-12 cells, as measured in a clonogenic assay . Finally, another inhibitor of AKR catalytic activity with a structure very distinct from cyclosporine A (flufenamic acid) also restored doxorubicin cytotoxicity and nuclear localization in MCF-7DOX2-12 cells (data not shown). This suggests that it is the ability of these agents to inhibit AKR activity that is responsible for the restoration of drug cytotoxicity. An alternative argument is that the doxorubicinol, once formed, is further metabolized, such that the metabolite is not retained in the method used to extract cellular doxorubicin and doxorubicinol for HPLC-based measurements. Thus, doxorubicinol would not be seen to accumulate in MCF-7DOX2-12 cells.
Despite the ability of both cyclosporin A and 5β-cholanic acid to increase cellular doxorubicin content in MCF-7DOX2-12 cells, why was only the latter agent able to appreciably restore doxorubicin cytotoxicity? Increasing the cellular content of doxorubicin by the cyclosporine-mediated reduction of drug efflux may not sufficiently increase its cytotoxicity if the additional cellular doxorubicin is rapidly converted to doxorubicinol by the elevated expression of AKRs and/or if the additional doxorubicin is sequestered into lysosomes. In contrast, AKR inhibition may block all conversion of doxorubicin to doxorubicinol, such that any drug entering the cell remains as doxorubicin and is able to rapidly reach the nucleus, before being sequestered.
Using a full genome approach, this study provides important new insight into pharmacokinetic and pharmacodynamic pathways that are altered upon selection of cells for resistance to doxorubicin. In addition to our previously reported finding of increased expression of the AKR 1C isoforms , the current study reveals other changes in gene expression that would be expected to affect the cytotoxicity of doxorubicin. This includes genes that may: decrease uptake of doxorubicin (SLC22A15), enhance efflux of doxorubicin (ABCC1, ABCG2, ABCD3, ABCA1), enhance conversion of doxorubicin to doxorubicinol (AKR1B10, AKR1B1), doxorubicin deoxyaglycone or doxorubicin semiquinone (NQO1), and inhibit the ability of doxorubicin to damage tumour cells through the generation of reactive oxygen species (CAT). Moreover, this study provides an in-depth comparison of the biochemical properties of doxorubicin versus doxorubicinol. While the former is highly cytotoxic, has high DNA binding affinity, and localizes to the nucleus in wildtype breast tumour cells, doxorubicinol is over a million times less cytotoxoic, has significantly reduced DNA binding activity, and is retained in the cytoplasm or lysosomes of cells. We also show that the administration of AKR inhibitors with doxorubicin in MCF-7DOX2 cells substantially restores both drug localization to the nucleus and drug cytotoxicity. Interestingly, doxorubicinol is highly cardiotoxic, and it is believed that doxorubicinol is responsible for the cardiotoxicity associated with doxorubicin chemotherapy [39, 40]. Since the AKR inhibitor 5β-cholanic acid is a well-tolerated naturally occurring bile acid in humans, and since flufenamic acid has been used in clinical trials with manageable toxicities , there may be significant value in conducting clinical trials in which either 5β-cholanic acid or flufenamic acid are co-administered with doxorubicin during chemotherapy. Results in this study would suggest that these AKR inhibitors may increase tumour levels of doxorubicin and block cardiotoxicity induced by doxorubicin conversion to doxorubicinol. This may dramatically improve the therapeutic index of doxorubicin when administered to cancer patients and improve the duration of clinical response for this otherwise highly effective chemotherapy drug.
Supplies and reagents
Supplies and reagents used in this study came from a variety of sources. Unless otherwise noted, Sigma (St. Louis, MO) was the supplier.
MCF-7 breast adenocarcinoma cells were obtained from the American Tissue Culture Collection (ATCC) (lot HTB-22) and selected for resistance to doxorubicin (Pfizer Pharmaceuticals, St. Laurent, QC) as previously described . Briefly, doxorubicin-sensitive, wildtype MCF-7 cells were grown in progressively increasing concentrations of doxorubicin from 1000x below the IC50 for the drug in parental MCF-7 cells (dose 1) to its maximally tolerated dose (dose 12) in 1.5- or 3-fold increments, with retention of cells surviving the greater of the two doses. Cells selected for survival in the varying doses of doxorubicin were termed MCF-7DOX2 cells. A “co-cultured control” cell line was selected under identical conditions in the absence of drug (MCF-7CC cells). These cells served as a control to help identify changes in gene expression due to long-term cell culture. The highest dose level to which cells were selected are indicated in the subscript of the cell line name. For example, MCF-7DOX2-12 cells refers to cells selected to the 12th dose level of doxorubicin. The 2 in the subscript is to prevent confusion with a previously isolated doxorubicin-resistant cell line in our laboratory (MCF-7DOX). All cells used in this study were selected to dose level 12 (0.1 μM doxorubicin). Cells were grown in high-glucose DMEM medium (Fisher Scientific, Nepean, ON) supplemented with penicillin-streptomycin (Fisher) and 10% fetal bovine serum (Fisher) in 75 cm2 tissue culture flasks (Sarstedt Canada, Montreal, QC), unless otherwise noted. Cells were maintained at 37°C in air supplemented with 5% CO2 in a humidified environment. Cells were passaged weekly, with a medium change once between passages. Drug-resistant cells were maintained in medium containing doxorubicin at their selection dose.
Changes in gene expression between MCF-7CC12 and MCF-7DOX2-12 cells were identified by microarray analysis using Agilent 4x44k whole human genome arrays (product number G4112F; Agilent Technologies, Mississauga, ON). These arrays enabled us to determine the level of expression of 27,958 human Entrez genes (close to the entire genome). Five hundred ng of total RNA, isolated with a Qiagen RNeasy kit (Mississauga, ON), was used for each sample. The RNA was then labeled with Cy3 or Cy5 using an Agilent Quick Amp labeling kit (Product #5190-0444). Hybridization was performed as per the manufacturer’s protocol. Experiments were repeated using multiple batches of labeled RNA, with both forward and reverse-labeling to account for dye bias, for a total of 16 two-colour arrays. The microarrays were scanned, and feature extraction and background intensity corrections were performed with Agilent software (v.10.7.3.1). Using Partek Genomics suite (St. Louis, MO) to perform a 4 way ANOVA using the Method of Moments , a list of genes significantly over- or under-expressed in MCF-7DOX2-12 cells relative to MCF-7CC12 cells. The false discovery rate was set at 0.01, with only genes changing expression by ≥ 2-fold being noted. The four variables assessed in the 4 way ANOVA were the cell line (MCF-7CC12 cells versus MCF-7DOX2-12 cells), the dye used (Cy3 versus Cy5), the experimental batch of arrays (to address batch effects) and the arrays themselves to address random effects. The input file was the data from all 16 two-colour arrays comparing gene expression between MCF-7DOX2-12 and MCF-7CC12 cells. The model used was: Yijklm = μ + Cell line (i) + Dye (j) + Exp batch (k) + arrays (random effect) (kl) + εijklm, where Yijklm represents the mth observation on the ith Cell line jth Dye kth Exp batch lth arrays, μ is the common effect for the whole experiment, εijklm represents the random error present in the mth observation, on the ith Cell line, jth Dye, kth Exp batch, lth arrays. The errors εijklm were assumed to be normally and independently distributed , with mean 0 and standard deviation δ for all measurements. Arrays and Exp batch were considered random effects. Normalized expression was transformed to the base 2.0, with p values reported for significance of differences in the expression of each gene. The output of the analysis was a p value for the significance of the observed fold change in expression of a particular gene and for the mean ratio of expression of a particular gene between the two cell lines. The microarray data was deposited in the NCBI Gene Expression Omnibus (GEO) database, accession number GSE27254) in accordance with MIAME standards . The url to access this data is: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=dbezngycywquuhm&acc=GSE27254.
Agilent raw data was converted to .MEV format using TIGR Express Converter v2.1.
The MEV files were then normalized using TIGR MIDAS v2.22 and a LOWESS normalization filter with flip-dye consistency checking where appropriate.
A one-class Significance Analysis of Microarrays (SAM) analysis was performed using TIGR MEV v4.6.1.
For pathway analysis, the list of significant differentially-expressed genes generated by Partek Genomics Suite was compared to a curated list of genes, transcripts or proteins shown to be involved in doxorubicin pharmacokinetics and pharmacodynamics in tumour cells or cardiomyocytes, available from the Pharmacogenomics Knowledge Base (PharmGKB) [16, 28]. These lists were then compared in S-Plus (v8.0) using the binomial statistical test as previously described .
RNA extraction, reverse transcription, and quantitative polymerase chain reaction
Given that the 1C and 1B AKR isoforms are highly conserved and that the probes for the 1C and 1B AKR transcripts on the Agilent 4X44 microarrays were not isoform-specific, we designed isoform-specific 1C and 1B primers to accurately quantify the levels of expression of the various 1C and 1B transcripts. These primers (Table 4) and isoform-specific primers for the carbonyl reductases (which also convert doxorubicin to doxorubicinol) were used in RTqPCR experiments. Total RNA was extracted from the MCF-7CC12 and MCF-7DOX2-12 cell lines using a Qiagen RNeasy kit (Mississauga, ON), reverse transcribed, and the cDNAs amplified using an ABI 7900HT quantitative PCR machine and SYBR Green 1 detection chemistry as described previously . After RNA extraction and before reverse transcription, RNA was quantified and quality ensured using an Agilent Bioanalyzer 2100 RNA nano kit. 2 μg of RNA was then DNase I (Invitrogen) treated. For reverse transcription, either MMLV reverse transcriptase (Invitrogen, Burlington, ON) and an oligo-dT20 primer, or Superscript III reverse transcriptase (Invitrogen) and an AKR1C-specific reverse transcription primer were used, all according to the manufacturer’s protocol. The cDNA was then stored at −20°C until analysis via the ΔCT method. All experiments were performed according to MIQE standards . Primers used for reverse transcription or to amplify specific cDNAs are described in Additional file 4: Table S4.
Protein extraction and quantification
Ten cm plates of cells (cultured for 2 passages without drug to 80% confluence) were rinsed twice in Dulbecco’s PBS, removed of excess liquid using a pipette tip, and placed on ice. Three hundred to 1000 μl of RIPA buffer (10 mM Tris–HCl, 1% sodium deoxycholate, 0.1% SDS, 1% Triton X-100, 150 mM NaCl, pH 7.5) with added Complete™ protease inhibitors (Roche Diagnostics, Laval, QC) were added to each plate, and the plates scraped with a cell scraper. The resulting crude lysate was passed through a 21 gauge needle 5 times. The lysate was incubated on ice for 30 minutes and then centrifuged for 20 minutes at 13.2 × 103 ×g at 4°C. The supernatant was retained and stored at −80°C. Protein concentrations in extracts were measured using a BCA protein quantification kit (Pierce Bioscience) using standard solutions of Bovine Serum Albumin (BSA).
Forty μg of protein were diluted in 6x Laemelli loading buffer and loaded onto a 10% SDS-PAGE gel with a 4% stacking gel. The gel was also loaded with 5μL of Bio-Rad (Mississagua, ON) dual colour protein marker. The gels were subjected to electrophoresis in a Bio-Rad mini tetra system at 80 V for 30 minutes and then at 120 V for an additional 1 h. The gels were removed and the gel proteins transferred to nitrocellulose membranes (GE Healthcare, Baie d’Urfe QC) using a Bio-Rad semi-dry electroblotting apparatus for 1 h at 12v. Membranes were then stained with 0.5% Ponceau S in 1% acetic acid to confirm transfer efficiency and even protein loading. The membranes were blocked in 5% skim milk powder (Carnation) in 0.1% TBST (20 mM Tris–HCl, 150 mM NaCl, 0.05% Tween 20, pH 7.5) for 1 h at room temperature, and then incubated overnight at 4°C with an isoform specific mouse monoclonal anti-AKR1C3 antibody  (Sigma) at a 1:2,000 dilution in 5% skim milk powder in TBST. The membranes were then washed with TBST for 15 minutes, and incubated in HRP-conjugated goat anti-mouse secondary antibody (Santa Cruz Biotechnology Inc., Santa Cruz, CA) at a 1:10,000 dilution in 5% skim milk powder in TBST for 1 h at room temperature. Membranes were again washed in TBST for 15 minutes and subjected to 3 × 5 minute final washes in TBS before being covered in ECL solution (Santa Cruz Biotechnology Inc.) and imaged using a gel documentation system (Alpha Innotech) for 10 minutes. An identical procedure using a mouse monoclonal anti-β-tubulin antibody (Santa Cruz Biotechnology Inc.) at a 1:10,000 dilution was used to monitor β-tubulin levels in the extracts as a loading control.
Cells were plated on #1 coverslips placed in 6-well plates and allowed to adhere overnight before being treated. Treatment consisted of addition of 0.5 μM doxorubicin or doxorubicinol (Toronto Research Chemicals Inc., North York, ON), and either 200 μM 5β-cholanic acid (Steraloids, Newport, RI) or DMSO as a vehicle control. After the 24 h treatment, DRAQ5 (Biostatus, Leicestershire, UK) was added to the culture media for 15 minutes as a nuclear counterstain. The coverslips were rinsed gently in 3 sequential PBS washes and sealed onto standard microscope slides using clear nail polish. After the nail polish dried, cells were observed using a Zeiss LSM 510 META confocal laser scanning microscope using an argon-ion laser at a 488 nm wavelength band for excitation of doxorubicin and doxorubicinol and using a 560 nm long-pass filter to detect intrinsic fluorescence of doxorubicin and its metabolites. A 633 nm laser with a 650 nm long-pass filter was used to detect DRAQ5 fluorescence.
High-performance liquid chromatography
Cells (8.0 × 106 cells per 10 cm plate) were allowed to adhere overnight, after which they were treated with 0.5 μM doxorubicin or 0.5 μM doxorubicinol (with or without 5 μM cyclosporine A and/or 200 μM 5β-cholanic acid) for 24 h. After this time period, the media was decanted (with 0.5 mL reserved for HPLC analysis), and the plates were rinsed twice in PBS. One mL of a 0.2 M Na2HPO4 solution, pH 8.5, was added to the plates and the cells were scraped off of the plate. A 0.5 ml volume of the same solution was added to the 0.5 mL of reserved media. Each sample was then added to 4 mL of a 9:1 v/v chloroform:n-heptanol mixture in a polypropylene 15 mL centrifuge tube and shaken on a mixer for 20 minutes, after which the samples were centrifuged for 10 minutes at 2000× g at 20°C. The bottom organic layer was then aspirated from the tube using a glass 5 mL pipette and dispensed into a new 15 mL centrifuge tube containing 250uL of 0.1 M ortho-phosphoric acid. Each tube was then mixed on a vortex mixer for 30 seconds before being centrifuged for 2 minutes at 2000× g. The top 200μL of the upper aqueous layer was then removed and stored at −80 degrees Celsius for later analysis.
Separations were performed using a revised gradient elution based on a previously described isocratic method  on a Waters Alliance e2695 system with a Waters 2475 fluorescence detector set at 480 nm excitation and 560 nm emission. Chromatographic conditions were the following: column: YMC CN 25 × 5 mm column; Eluent A: 10 mM NaH2PO4 pH 4.0, Eluent B: HPLC grade CH3CN; flow rate: 1.0 mL/min. The gradient program was as follows: 0 min = 20% B 80% A, 10 min = 50% B 50% A, 11 to 24 min = 20% B 80% A. The slope for each gradient change was linear.
DNA binding affinity assay
The relative DNA binding affinity of doxorubicin and doxorubicinol was determined by using a fluorescent intercalator displacement assay . Briefly, a quartz cuvette was filled with 3 mL of Tris buffer (0.1 M Tris, 0.1 M NaCl, pH 8.0; approximating conditions inside the nucleus) to which 4.4 μM ethidium bromide was added. A fluorescence reading (excitation: 545 nm, emission: 590 nm) was taken using a Perkin Elmer LS-50 fluorimeter; this constituted the baseline reading. Pre-sheared salmon sperm DNA (8.8 μM in base pairs) was then added to the cuvette, incubated for 5 minutes, and again the fluorescence was determined; this constituted the maximal or 100% reading. Aliquots of doxorubicin or doxorubicinol (0.067 μM) were added to the cuvette, incubated for 5 minutes, and the corresponding reading recorded. The background reading was subtracted for each reading and then divided by the maximal reading to determine per cent of maximal binding. These data were fit to curves to determine Kapp and Bmax values.
Measurement of drug sensitivity
Drug sensitivity was assessed using a variation  of the standard clonogenic assay . Briefly, for each condition, 12 × 25 cm2 flasks were plated with 2.5 × 105 cells and left to adhere overnight. The next day, each flask was treated with a different concentration of doxorubicin, decreasing in 3-fold increments, from 3.0 × 10-6 M to 5.13 × 10-11 M, with a final flask receiving no doxorubicin. After 24 h, cells were trypsinized, pelleted, and resuspended in 300uL of medium which was then combined with 2.7 mL of methyl cellulose growth medium [2.6% methyl cellulose, (Shin-Etsu) and 30% FBS in IMDM (Princess Margaret Hospital)]. After being mixed thoroughly, the suspension was allowed to settle for 30 minutes before 1.2 ml of cells were introduced into 6-well tissue culture plates. Plates were incubated for 2 weeks and then 10 randomly selected fields in each well were counted at 40x magnification.
Graphpad Prism (v5.0) was used for all statistical tests unless otherwise noted. Differences between treatment means were assessed using either a Student’s unpaired t-test or an unpaired 1-way Analysis of Variance (ANOVA) with Tukey’s Honestly Significant Difference (HSD) post-hoc test where appropriate. A p value ≤ 0.05 was considered significant.
Source of funding
Canadian Institutes of Health Research (Grant MOP-8993 to A.M.P.)
ATP-binding cassette protein C1
Analysis of variance
Measure of the number of receptors for a particular drug
High performance liquid chromatography
Concentration required to inhibit growth by 50%
Apparent dissociation constant
Minimum information associated with a microarray experiment
NAD(P)H dehydrogenase, quinine 1
Polymerase chain reaction
This work was supported by a grant from the Canadian Institutes of Health Research (MOP-89993). Core funding to support this work is also acknowledged from the Northern Cancer Foundation.
- Schneider YJ, Baurain R, Zenebergh A, Trouet A: DNA-binding parameters of daunorubicin and doxorubicin in the conditions used for studying the interaction of anthracycline-DNA complexes with cells in vitro. Cancer Chemother Pharmacol. 1979, 2: 7-10.View ArticlePubMed
- Foglesong PD, Reckord C, Swink S: Doxorubicin inhibits human DNA topoisomerase I. Cancer Chemother Pharmacol. 1992, 30: 123-125. 10.1007/BF00686403.View ArticlePubMed
- Lopez M: Anthracyclines in the adjuvant treatment of breast carcinoma: thirty years later. Clin Ter. 2006, 157: 165-177.PubMed
- Kartner N, Riordan JR, Ling V: Cell surface P-glycoprotein associated with multidrug resistance in mammalian cell lines. Science. 1983, 221: 1285-1288. 10.1126/science.6137059.View ArticlePubMed
- Doyle LA, Ross DD: Multidrug resistance mediated by the breast cancer resistance protein BCRP (ABCG2). Oncogene. 2003, 22: 7340-7358. 10.1038/sj.onc.1206938.View ArticlePubMed
- Frank NY, Margaryan A, Huang Y, Schatton T, Waaga-Gasser AM, Gasser M, Sayegh MH, Sadee W, Frank MH: ABCB5-mediated doxorubicin transport and chemoresistance in human malignant melanoma. Cancer Res. 2005, 65: 4320-4333. 10.1158/0008-5472.CAN-04-3327.View ArticlePubMed
- Borst P, Evers R, Kool M, Wijnholds J: A family of drug transporters: the multidrug resistance-associated proteins. J Natl Cancer Inst. 2000, 92: 1295-1302. 10.1093/jnci/92.16.1295.View ArticlePubMed
- Gavelova M, Hladikova J, Vildova L, Novotna R, Vondracek J, Krcmar P, Machala M, Skalova L: Reduction of doxorubicin and oracin and induction of carbonyl reductase in human breast carcinoma MCF-7 cells. Chem Biol Interact. 2008, 176: 9-18. 10.1016/j.cbi.2008.07.011.View ArticlePubMed
- Hurwitz SJ, Terashima M, Mizunuma N, Slapak CA: Vesicular anthracycline accumulation in doxorubicin-selected U-937 cells: participation of lysosomes. Blood. 1997, 89: 3745-3754.PubMed
- Rajagopal A, Simon SM: Subcellular localization and activity of multidrug resistance proteins. Mol Biol Cell. 2003, 14: 3389-3399. 10.1091/mbc.E02-11-0704.PubMed CentralView ArticlePubMed
- Aas T, Borresen AL, Geisler S, Smith-Sorensen B, Johnsen H, Varhaug JE, Akslen LA, Lonning PE: Specific P53 mutations are associated with de novo resistance to doxorubicin in breast cancer patients. Nat Med. 1996, 2: 811-814. 10.1038/nm0796-811.View ArticlePubMed
- Fischer V, Einolf HJ, Cohen D: Efflux transporters and their clinical relevance. Mini Rev Med Chem. 2005, 5: 183-195.View ArticlePubMed
- Leonard GD, Fojo T, Bates SE: The role of ABC transporters in clinical practice. Oncologist. 2003, 8: 411-424. 10.1634/theoncologist.8-5-411.View ArticlePubMed
- Cho RJ, Campbell MJ: Transcription, genomes, function. Trends Genet. 2000, 16: 409-415. 10.1016/S0168-9525(00)02065-5.View ArticlePubMed
- Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29: 365-371. 10.1038/ng1201-365.View ArticlePubMed
- Thorn CF, Oshiro C, Marsh S, Hernandez-Boussard T, McLeod H, Klein TE, Altman RB: Doxorubicin pathways: pharmacodynamics and adverse effects. Pharmacogenet Genomics. 2011, 21: 440-446. 10.1097/FPC.0b013e32833ffb56.PubMed CentralView ArticlePubMed
- Veitch ZW, Guo B, Hembruff SL, Bewick AJ, Heibein AD, Eng J, Cull S, Maclean DA, Parissenti AM: Induction of 1C aldoketoreductases and other drug dose-dependent genes upon acquisition of anthracycline resistance. Pharmacogenet Genomics. 2009, 19: 477-488. 10.1097/FPC.0b013e32832c484b.View ArticlePubMed
- Minotti G, Menna P, Salvatorelli E, Cairo G, Gianni L: Anthracyclines: molecular advances and pharmacologic developments in antitumor activity and cardiotoxicity. Pharmacol Rev. 2004, 56: 185-229. 10.1124/pr.56.2.6.View ArticlePubMed
- Penning TM, Drury JE: Human aldo-keto reductases: Function, gene regulation, and single nucleotide polymorphisms. Arch Biochem Biophys. 2007, 464: 241-250. 10.1016/j.abb.2007.04.024.PubMed CentralView ArticlePubMed
- Kassner N, Huse K, Martin HJ, Godtel-Armbrust U, Metzger A, Meineke I, Brockmoller J, Klein K, Zanger UM, Maser E, Wojnowski L: Carbonyl reductase 1 is a predominant doxorubicin reductase in the human liver. Drug Metab Dispos. 2008, 36: 2113-2120. 10.1124/dmd.108.022251.View ArticlePubMed
- Novotna R, Wsol V, Xiong G, Maser E: Inactivation of the anticancer drugs doxorubicin and oracin by aldo-keto reductase (AKR) 1C3. Toxicol Lett. 2008, 181: 1-6. 10.1016/j.toxlet.2008.06.858.View ArticlePubMed
- Chang BK, Brenner DE, Gutman R: Cellular pharmacology of doxorubicinol alone and combined with verapamil in pancreatic cancer cell lines. Anticancer Res. 1989, 9: 341-345.PubMed
- Bernardini N, Giannessi F, Bianchi F, Dolfi A, Lupetti M, Zaccaro L, Malvaldi G, Del TM: Comparative activity of doxorubicin and its major metabolite, doxorubicinol, on V79/AP4 fibroblasts: a morphofunctional study. Exp Mol Pathol. 1991, 55: 238-250. 10.1016/0014-4800(91)90004-H.View ArticlePubMed
- Endo S, Matsunaga T, Mamiya H, Ohta C, Soda M, Kitade Y, Tajima K, Zhao HT, El-Kabbani O, Hara A: Kinetic studies of AKR1B10, human aldose reductase-like protein: endogenous substrates and inhibition by steroids. Arch Biochem Biophys. 2009, 487: 1-9. 10.1016/j.abb.2009.05.009.View ArticlePubMed
- Steckelbroeck S, Oyesanmi B, Jin Y, Lee SH, Kloosterboer HJ, Penning TM: Tibolone metabolism in human liver is catalyzed by 3alpha/3beta-hydroxysteroid dehydrogenase activities of the four isoforms of the aldo-keto reductase (AKR)1C subfamily. J Pharmacol Exp Ther. 2006, 316: 1300-1309.View ArticlePubMed
- Coley HM, Amos WB, Twentyman PR, Workman P: Examination by laser scanning confocal fluorescence imaging microscopy of the subcellular localisation of anthracyclines in parent and multidrug resistant cell lines. Br J Cancer. 1993, 67: 1316-1323. 10.1038/bjc.1993.244.PubMed CentralView ArticlePubMed
- Hembruff SL, Laberge ML, Villeneuve DJ, Guo B, Veitch Z, Cecchetto M, Parissenti AM: Role of drug transporters and drug accumulation in the temporal acquisition of drug resistance. BMC Cancer. 2008, 8: 318-10.1186/1471-2407-8-318.PubMed CentralView ArticlePubMed
- Klein TE, Chang JT, Cho MK, Easton KL, Fergerson R, Hewett M, Lin Z, Liu Y, Liu S, Oliver DE, Rubin DL, Shafa F, Stuart JM, Altman RB: Integrating genotype and phenotype information: an overview of the PharmGKB project. Pharmacogenetics Research Network and Knowledge Base. Pharmacogenomics J. 2001, 1: 167-170. 10.1038/sj.tpj.6500035.View ArticlePubMed
- Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio–the analysis and navigation of molecular networks. Bioinformatics. 2003, 19: 2155-2157. 10.1093/bioinformatics/btg290.View ArticlePubMed
- Joshi-Tope G, Gillespie M, Vastrik I, D'Eustachio P, Schmidt E, de BB, Jassal B, Gopinath GR, Wu GR, Matthews L, Lewis S, Birney E, Stein L: Reactome: a knowledgebase of biological pathways. Nucleic Acids Res. 2005, 33: D428-D432.PubMed CentralView ArticlePubMed
- Dahlquist KD, Salomonis N, Vranizan K, Lawlor SC, Conklin BR: GenMAPP, a new tool for viewing and analyzing microarray data on biological pathways. Nat Genet. 2002, 31: 19-20. 10.1038/ng0502-19.View ArticlePubMed
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: 3-10.1186/gb-2003-4-5-p3.View Article
- Okabe M, Unno M, Harigae H, Kaku M, Okitsu Y, Sasaki T, Mizoi T, Shiiba K, Takanaga H, Terasaki T, Matsuno S, Sasaki I, Ito S, Abe T: Characterization of the organic cation transporter SLC22A16: a doxorubicin importer. Biochem Biophys Res Commun. 2005, 333: 754-762. 10.1016/j.bbrc.2005.05.174.View ArticlePubMed
- Schriner SE, Linford NJ, Martin GM, Treuting P, Ogburn CE, Emond M, Coskun PE, Ladiges W, Wolf N, Van RH, Wallace DC, Rabinovitch PS: Extension of murine life span by overexpression of catalase targeted to mitochondria. Science. 2005, 308: 1909-1911. 10.1126/science.1106653.View ArticlePubMed
- Jin Y, Penning TM: Aldo-keto reductases and bioactivation/detoxication. Annu Rev Pharmacol Toxicol. 2007, 47: 263-292. 10.1146/annurev.pharmtox.47.120505.105337.View ArticlePubMed
- Ohara H, Miyabe Y, Deyashiki Y, Matsuura K, Hara A: Reduction of drug ketones by dihydrodiol dehydrogenases, carbonyl reductase and aldehyde reductase of human liver. Biochem Pharmacol. 1995, 50: 221-227. 10.1016/0006-2952(95)00124-I.View ArticlePubMed
- Martin HJ, Maser E: Role of human aldo-keto-reductase AKR1B10 in the protection against toxic aldehydes. Chem Biol Interact. 2009, 178: 145-150. 10.1016/j.cbi.2008.10.021.View ArticlePubMed
- Qadir M, O'Loughlin KL, Fricke SM, Williamson NA, Greco WR, Minderman H, Baer MR: Cyclosporin A is a broad-spectrum multidrug resistance modulator. Clin Cancer Res. 2005, 11: 2320-2326. 10.1158/1078-0432.CCR-04-1725.View ArticlePubMed
- Olson RD, Mushlin PS, Brenner DE, Fleischer S, Cusack BJ, Chang BK, Boucek RJ: Doxorubicin cardiotoxicity may be caused by its metabolite, doxorubicinol. Proc Natl Acad Sci USA. 1988, 85: 3585-3589. 10.1073/pnas.85.10.3585.PubMed CentralView ArticlePubMed
- Mordente A, Minotti G, Martorana GE, Silvestrini A, Giardina B, Meucci E: Anthracycline secondary alcohol metabolite formation in human or rabbit heart: biochemical aspects and pharmacologic implications. Biochem Pharmacol. 2003, 66: 989-998. 10.1016/S0006-2952(03)00442-8.View ArticlePubMed
- Kapadia L, Elder MG: Flufenamic acid in treatment of primary spasmodic dysmenorrhoea. A double-blind crossover study. Lancet. 1978, 1: 348-350.View ArticlePubMed
- Eisenhart C: The assumptions underlying the analysis of variance. Biometrics. 1947, 3: 1-21. 10.2307/3001534.View ArticlePubMed
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, Sturn A, Snuffin M, Rezantsev A, Popov D, Ryltsov A, Kostukovich E, Borisovsky I, Liu Z, Vinsavich A, Trush V, Quackenbush J: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34: 374-378.PubMed
- Hembruff SL, Villeneuve DJ, Parissenti AM: The optimization of quantitative reverse transcription PCR for verification of cDNA microarray data. Anal Biochem. 2005, 345: 237-249. 10.1016/j.ab.2005.07.014.View ArticlePubMed
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.View ArticlePubMed
- Lin HK, Steckelbroeck S, Fung KM, Jones AN, Penning TM: Characterization of a monoclonal antibody for human aldo-keto reductase AKR1C3 (type 2 3alpha-hydroxysteroid dehydrogenase/type 5 17beta-hydroxysteroid dehydrogenase); immunohistochemical detection in breast and prostate. Steroids. 2004, 69: 795-801. 10.1016/j.steroids.2004.09.014.View ArticlePubMed
- Boger DL, Fink BE, Brunette SR, Tse WC, Hedrick MP: A simple, high-resolution method for establishing DNA binding affinity and sequence selectivity. J Am Chem Soc. 2001, 123: 5878-5891. 10.1021/ja010041a.View ArticlePubMed
- Chadderton A, Villeneuve DJ, Gluck S, Kirwan-Rhude AF, Gannon BR, Blais DE, Parissenti AM: Role of specific apoptotic pathways in the restoration of paclitaxel-induced apoptosis by valspodar in doxorubicin-resistant MCF-7 breast cancer cells. Breast Cancer Res Treat. 2000, 59: 231-244. 10.1023/A:1006344200094.View ArticlePubMed
- Pollard JW, Walker JM: Animal Cell Culture. 1990, Clifton NJ: Humana Press, 5View Article
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/381/prepub
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.