Skip to main content

ALDH isozymes downregulation affects cell growth, cell motility and gene expression in lung cancer cells

Abstract

Background

Aldehyde dehydrogenase isozymes ALDH1A1 and ALDH3A1 are highly expressed in non small cell lung cancer. Neither the mechanisms nor the biologic significance for such over expression have been studied.

Methods

We have employed oligonucleotide microarrays to analyze changes in gene profiles in A549 lung cancer cell line in which ALDH activity was reduced by up to 95% using lentiviral mediated expression of siRNA against both isozymes (Lenti 1+3). Stringent analysis methods were used to identify gene expression patterns that are specific to the knock down of ALDH activity and significantly different in comparison to wild type A549 cells (WT) or cells similarly transduced with green fluorescent protein (GFP) siRNA.

Results

We confirmed significant and specific down regulation of ALDH1A1 and ALDH3A1 in Lenti 1+3 cells and in comparison to 12 other ALDH genes detected. The results of the microarray analysis were validated by real time RT-PCR on RNA obtained from Lenti 1+3 or WT cells treated with ALDH activity inhibitors. Detailed functional analysis was performed on 101 genes that were significantly different (P < 0.001) and their expression changed by ≥ 2 folds in the Lenti 1+3 group versus the control groups. There were 75 down regulated and 26 up regulated genes. Protein binding, organ development, signal transduction, transcription, lipid metabolism, and cell migration and adhesion were among the most affected pathways.

Conclusion

These molecular effects of the ALDH knock-down are associated with in vitro functional changes in the proliferation and motility of these cells and demonstrate the significance of ALDH enzymes in cell homeostasis with a potentially significant impact on the treatment of lung cancer.

Background

Aldehyde dehydrogenases (ALDHs) are a group of NAD(P)+-dependent enzymes involved in the metabolism of a wide variety of aliphatic and aromatic aldehydes [1, 2]. Many disparate aldehydes are ubiquitous in nature and are toxic at low levels because of their chemical reactivity. Thus levels of metabolic-intermediate aldehydes must be carefully regulated which explains the existence of several distinct ALDH families in most studied organisms with wide constitutive tissue distribution [1, 2]. A systematic nomenclature scheme for the ALDH gene superfamily based on divergent evolution has been developed [3] and continues to be updated on paper [4] and on the internet by Dr. Vasilis Vasiliou and his group http://www.aldh.org. According to the latest database, the human genome contains 19 ALDH functional genes and three pseudogenes [4].

The role of some of these ALDHs in endobiotic and xenobiotic metabolism has been reviewed extensively before and the specific metabolic pathways affected have been detailed [2]. Many allelic variants within the ALDH gene family have been identified, resulting in pharmacogenetic heterogeneity between individuals which, in most cases, results in distinct phenotypes [2, 5] including intolerance to alcohol and increased risk of ethanol-induced cancers (ALDH2 and ALDH1A1), Sjogren-Larson Syndrome (ALDH3A1), type II hyperprolinemia (ALDH4A1), 4-hydroxybutyric aciduria (ALDH5A1), developmental delay (ALDH6A1), hyperammonemia (ALDH18A1), and late onset of Alzheimer's disease (ALDH2). Furthermore, knockouts of ALDH1A2 and ALDH1A3 in mouse are embryonic lethal and newborn lethal, respectively [6–8]. Changes in ALDH activity have been observed during experimental liver and urinary bladder carcinogenesis and in a number of human tumors [9].

One of the well studied pathways of ALDH activity is drug resistance to oxazaphosphorines. We have been interested in the role of ALDH 1A1 in drug resistance, first in hematopoietic progenitors and more recently in lung cancer. ALDH1A1, ALDH3A1, and ALDH5A1 have been shown to catalyze the oxidation of aldophosphamide [10–12]. We and others have shown that overexpression of ALDH1A1 and ALDH3A1 results in resistance to 4-hydroperoxycyclophosphamide (4-HC), an active derivative of cyclophosphamide (CP) [9–11, 13, 14]. More recently, ALDH3A1 was recognized as an oxidative stress response protein and thus can protect against the oxidative damage caused by other chemotherapy drugs such as etoposide [15]. We have also shown that down regulation of each enzyme by RNA antisense (AS) [16], all-trans retinoic acid (ATRA) [17] or siRNA [18] results in increased sensitivity to 4-HC.

Tetraethylthiuram disulfide (TT) (disulfiram, also known as Antabuse), an ALDH inhibitor, has been reported to affect the growth of multiple tumor cells, inhibit cancer cell invasiveness, and induce apoptosis using in vitro assays [19]. These effects were thought to be due to different mechanisms including inhibition of proteasome activity [20], increase Cu uptake with pro-oxidant effects [21, 22], inhibition of NF κB [23–25], inhibition of the relaxation activity of DNA topoisomerases I and II [26], and inhibition of caspases [27].

All of the above studies indicate the biologic and clinical significance of these enzymes and, therefore, the need to better define the regulatory mechanisms involved in determining their level of expression in normal and malignant tissues. Multiple studies, mainly in animal models, have been published on the regulation of the various ALDH isozymes [28–31]. Functional genomics aim at analyzing the regulation of genes in response to physiological changes. Microarray technology revolutionized the analysis of gene expression in biological processes to enable the assessment of gene activity on a genome-wide scale. In order to be able to perform such experiment in relation to ALDH1A1 and ALDH3A1, we have aimed at achieving "knock-down" of these enzymes using siRNA approach in vitro. Indeed, we achieved > 95% "knock-down" of ALDH activity in A549 lung cancer cell line using lentiviral vectors to permanently express siRNA sequences specific to each one of the enzymes. The data presented here reveal multiple genes affected by the "knock-down" of ALDH activity which will ultimately aid in identifying the molecular events related to the activity and regulation of ALDHs expression be it in oxidative stress, response to carcinogenic aldehydes or malignant transformation.

Methods

Cell Lines

A549 and H522 cell lines were obtained from ATCC and maintained in -80°C freezer or cultured in RPMI 1640 medium containing 10% FBS before experiments. Cells were maintained in 5% CO2 incubator at 37°C and used for experiments in the exponential phase of their growth. These cell lines were specifically chosen for the described experiments because they are known to have high levels of ALDH1A1 and ALDH3A1 [32]. These cells were transduced with lentiviral vectors (described below) containing specific siRNA sequences against ALDH1A1 (Lenti 1 cells), ALDH3A1 (Lenti 3 cells), both vectors (Lenti 1+3 cells), and against the green fluorescent protein (GFP) gene (GFP cells, used as a control). Overall total of 5 cell lines from A549 and H522 were used throughout the experiments described here including the parent wild type cell line (WT). Microarray gene profiling was performed only on A549 cells. In order to further validate that the effects seen in the ALDH knock down cells are directly related to ALDH activity, we used known ALDH inhibitors as described below in A549 and H522 cells as well as one other lung cancer cell line, H1299, that lacks any significant expression of ALDH1A1 or ALDH3A1 [32].

Effects of ALDH Knock-Down on Cell Growth and Motility

The assays described below were used to study the effects of ALDH knock down by siRNA on cell growth and migration as well as colony formation. Triplicates of 2 × 105 cells/ml/well were plated in 6-well tissue culture plates and grown for 72 hours. After this time, the cells were trypsinized, collected and manually counted using hematocytometer and microscope. Viability was determined by trypan blue exclusion. For these experiments, we used the five cell lines from H522 or A549. Means ± SD of growth rates were calculated and differences in growth rates were compared using the Student t-test.

A colorimetric assay kit (MTT based) for the non-radioactive quantification of cell proliferation and viability was purchased from Roche Diagnostics Corporation (Indianapolis, IN) and used for the cell proliferation experiments. In a 96-well plate, 1 × 103 cells/100 μL/well of each cell line were cultured for 96 hours. After the 96 hour incubation period, the labeling reagent was added to each well, followed by the solubilization solution four hours later. The cells were left overnight in the incubator, and absorbance was measured the following day using a Molecular Devices Kinetic Microplate reader at a wavelength of 570 nm as per the manufacturer's protocol. In this experiment, absorbance directly correlates with number of viable cells per well. Mean ± SD of the absorbance values were calculated for each experimental group and compared using the Student t-test.

For each of the cell lines, triplicates of 200 cells/ml/well in a six-well plate were cultured in RPMI-1640 + 10% FBS. The cells were allowed to grow for 5 days. The number of colonies adhered to the bottom of the plate was then counted using inverted microscope. The mean number of colonies was calculated and compared among all five cell lines originating from the A549 or H522 cell lines.

We measured the effect of down regulation of ALDH on cell migration and motility by using the in vitro scratch wound migration assay [33]. A549 monolayer confluent cells from the five different cell lines were scratched using sterile 200 μl pipette tip in 6-well plastic dishes, and after 20 hr of culture in RPMI-1640 +10% FBS, the migration ability of the cells was evaluated by the width of the wound under inverted microscope fitted with digital camera. Images of the scratch wound were taken immediately after the scratch and 20 hr later. The migration distances (in cm) of the cells were measured on printed images taken at the same magnification and the difference of the width of wounds at 0 and 20 hr was determined and the % wound healing was calculated.

SiRNA Plasmid Constructs

SiRNA constructs were made using pSilencer 2.0-U6 as vector backbone from Ambion. To generate the lentiviral siRNA constructs, the pSilencer 2.0-U6 siRNA plasmid was digested with Pvu II and the U6-siRNA cassette was cloned into the self-inactivating (SIN) pTYF-EFnlacZ lentiviral vector in the Not I site which was blunted with T4 DNA polymerase. For the constructions of ALDH LV siRNA vectors, the target site for ALDH1A1 (Genbank NM_000689) is: 5'-gtagccttcacaggatcaa-3' (nt 777–795) was chosen, and two primers were used to generate the siRNA construct; primer 1: 5'ATA CGC GGA TCC CGT AGC CTT CAC AGG ATC AAT TCA AGA G AT TGA TC-3', primer 2: 5'CGC TAG ACT AGT TAA AAA AGT AGC CTT CAC AGG ATC AAT CTC TTG AA TTG ATC -3'; the target site for ALDH3A1 (Genbank NM_000691) is: 5'-gaagatgattgcagagaca-3' (nt1167-1185), and two primers were used to generate the siRNA construct: primer 1: 5' ATA CGC GGA TCC CGA AGA TGA TTG CAG AGA CAT TCA AGA GAT GTC T-3', primer 2: 5' CGC TAG ACT AGT TAA AAA AGA AGA TGA TTG CAG AGA CAT CTC TTG AA TGT CT -3'. As a control, we similarly generated U6 promoter-driven LV-siRNA targeting GFP, as described before [34]. The lentiviral constructs with the U6-siRNA inserted in the reverse orientation were chosen for vector preparation (see Figure 1).

Figure 1
figure 1

Lentiviral constructs. Schematic presentation of the Lentiviral constructs used in our experiments to express the various siRNA sequences against ALDH1A1, ALDH3A1, and GFP.

Lentiviral Vector Production, Titration and Transduction

To produce lentiviral vectors, 293T cells were co-transfected with the three vector plasmids: pNHP, pHEFVSVG, and pTYF siRNA vector, and the virus supernatants were concentrated and titered as previously described [35]. For lentiviral transduction, the cells were plated at a concentration of 105 cells per well in 12-well plate to give 60–80% confluence after overnight culture and then infected at MOI of 50–100 in the presence of polybrene (8 μg/ml). Transduction efficiency was determined by the LacZ reporter gene assay as previously described [35].

After successful transduction with lentiviral constructs, we plated cells for liquid colony culture assay at 200 cells in 1 ml per 35 mm plate. On approximately day 5 of culture, we picked single colonies using inverted microscope and sterile pipette for replating and establishing clones for each of Lenti 1, Lenti 3, Lenti 1+3, and GFP cell lines. The 100% purity of the clones was validated by LacZ staining done on similarly cultured colonies from these established clones.

ALDH Activity and Protein Measurements

In order to measure the effects of the siRNA expression, we used the spectrophotometric methods for measurement of ALDH enzymatic activity as well Western blot analysis to detect changes in the protein levels of ALDH1A1 and ALDH3A1. Western blot analysis was performed as described before [17, 18]. In brief, lysates from each of the 5 cell lines mentioned above were size separated in parallel on two 12% denaturing SDS-polyacrylamide gels (BioRad, Hercules, CA), electrotransferred onto nitrocellulose membranes, blocked with 5% milk in TBS, and probed using chicken anti human ALDH1A1 and ALDH3A1 polyclonal antibodies provided generously by Dr. L. Sreerama (St Cloud University, Minneapolis, MN) and Dr. NE Sladek (University of Minnesota, Minneapolis, MN). The specificity of these antibodies has been documented by Dr. Sladek's group [36, 37]. Blots were incubated with chicken anti-human ALDH1A1 and ALDH3A1 primary antibodies at 1:200 or 1:300 dilution, respectively, for 1 hour at room temperature. After washing, the secondary antibody (horseradish peroxidase-labeled rabbit anti-chicken antibody; Sigma Chemical Co., St Louis, MO) was used at 1:4000 dilution for 1 hour. Chemiluminscence method (SuperSignal, Pierce, Rockford, IL) was used for the final visualization of the protein bands on X-ray film (Super Rx, Fuji Photo Film, Tokyo, Japan). After washing and blocking, the same blots were labeled again for visualization of actin as a loading control using anti-actin antibody (Oncogene Research Products, Cambridge, MA).

The cell lysates used for Western blot analysis, were also freshly used to measure ALDH enzyme activity using the spectrophotometric assay as described before [17, 18]. Briefly, the aliquots of 600 μl lysing buffer were incubated at 37°C in Beckman DLC 64 spectrophotometer cuvettes with the addition of cell lysate, 5 mM NAD+ and 5 mM propionaldehyde as a substrate. The rate of change in absorbance at 340 nm was measured in 3 replicates over 5 min. A control reaction in which the substrate was not added monitored the endogenous rate of NAD+ reduction. The ALDH activity was expressed in nmoles/107 cells.min.

RNA Isolation and cDNA Synthesis

A5 49 cell lines (WT, GFP, and Lenti 1+3) were plated in 25 cm2 Corning cell culture flasks at a density of 7.5 × 105 cells/ml with RPMI 1640 medium containing 10% FBS at 37°C. Four culture sets were prepared for each experimental group. Four different clones were used for GFP and Lenti 1+3. After 24 hr incubation, media was removed from the dish and RLT buffer was added directly to the cells. Cells were scraped into this buffer and subsequently homogenized using the QIAshredder kit (Qiagen, Valencia, CA). Total RNA was then harvested using the RNeasy RNA isolation kit (Qiagen). Integrity of total RNA was validated using a Bioanalyzer (Agilent Technologies) as well as OD 260/280 ratios.

Two micrograms of total RNA were utilized for the cDNA synthesis using the cDNA synthesis kit (Affymetrix 900431) and according to the Affymetrix technical protocols. First-strand synthesis was initiated using Superscript II Reverse Transcriptase and the T7-(dt) primer (Affymetrix 900431), which contains the T7 promoter sequence followed by an oligo (dt) tract. Second strand synthesis was carried out using Escherichia coli DNA polymerase I, E. coli DNA ligase, and RNase H (Affymetrix 900431). Five microliters of purified ds cDNA was used to generate biotin-labeled cRNA using the IVT labeling kit (Affymetrix 900449). Biotin labeled cRNA was then purified from unincorporated nucleotides using the RNA cleanup kit (Affymetrix 900371). Quantity and quality of cRNA was determined by OD260/280 ratios and by using a Bioanalyzer (Agilent Technologies).

Probe and Preparation/Staining of Array

Before using the array, biotin-labeled cRNA was fragmented for 35 minutes at 94°C. The following reagents were mixed to give the following final concentrations: fragmented cRNA (0.05 μg/ul), control oligonucleotide B2 (50 pM, Affymetrix), 20X eukaryotic hybridization controls (BioB at 1.5 pM, BioC at 5 pM, BioDn at 25 pM, and CreX at 100 pM, Affymetrix 900454), herring sperm DNA (0.1 mg/ml, Promega), acetylated BSA (50 mg/ml, GIBCO-BRL), and 2X hybridization buffer. Probe was heated to 99°C for 5 min, then transferred to 45°C for 5 min, and finally spun at maximum in a microcentrifuge for 5 min. Arrays used for this study were the Human Genome U133 Plus 2.0 Array (Affymetrix) representing over 47,000 transcripts and variants including 38,500 well-characterized human genes. Chips were interrogated with material harvested from A549 WT, Lenti 1+3 and GFP cells and each cell treatment was replicated four times using 4 chips for each experimental group, making the data from each array truly independent biological replicas. Each array was injected with 200 ul of target and rotated in a 45°C oven for 16 hours. Probe was removed from the array and the array was washed and stained using the Euk-WS2 fluidics protocol and the streptavidin-phycoerythrin (SAPE) reagent (Affymetrix).

Microarray Scanning,, Normalization, Expression Filters, Cluster Analysis and Statistical Analysis

Scanned images (*.dat files) were analyzed with GCOS software to generate a detection call for each probe set. Hybridization single intensities between GeneChips were normalized and an expression matrix was derived with dChip [38] using the PM only model model-based expression index. Probe sets whose hybridization single intensities were never measured above background levels were removed from the data set to reduce the noise level in the data set. Probe sets whose hybridization single intensities were measured at or below background levels were identified using the Affymetrix detection call algorithm [39].

Supervised learning, cross validation, and visualization

The modeled-based expression matrix of the probe sets passing the initial expression filter were analyzed in log space with BRB Array Tools 3.2.2 (developed by Richard Simon and Amy Peng Lam and is available online at the Biometric Research Branch of the National Cancer Institute website http://linus.nci.nih.gov/BRB-ArrayTools.html to identify genes that differentiated among the treatment classes at the P < 0.001 significance threshold. The ability of gene identified at the P < 0.001 significance threshold to be differentiated among treatment classes was assessed by using leave-one-out-cross-validation (LOOCV) using a nearest-neighbor prediction model. For visualization purposes the modeled-based expression matrix of probes sets identified as significant among the treatment groups were mean centered and variance normalized using algorithms implemented in dChip.

Functional Analysis of Transcripts

Probe sets that were significantly different at P < 0.001 and demonstrated ≥ 2 fold change in the level of gene expression were subject to comprehensive, detailed, six-level functional analysis, whereby identity was verified using GenBank accession number, aliases were compiled, and biochemical and biological functions were determined. After the identity of the probe set was confirmed, multiple internet bioinformatics and National Center for Biotechnology Information (NCBI) search engines, including PubMed searches, were used to assign the biological function. Probe sets without confirmed identities were placed in the "Unknown" category. Probe sets that could be identified but that did not have ascertainable function were classified as "Other"

Validation of Selected Gene Expression Levels by Quantitative Real-Time RT-PCR

In order to confirm microarray data, we performed real-time RT-PCR [40] for 7 selected genes measuring the change in their level of expression. Primer/probe sets for HMOX1, MAP3K8, MGST1, CXXC5, LHX8, RAB20, STRA6 and 18S rRNA (reference control) were designed by Assays-on-Demand (Applied Biosystems, Foster City, CA, USA) and are listed in Table 1. Total RNA was obtained as described above from all five A549 cells lines described under "Cell Lines".

Table 1 Primers used for quantitative real time RT-PCR.

Two micrograms of total RNA were reverse transcribed to cDNA in a 50 μl reaction containing 4 μl of 100 mM MgCl, 1.0 μl of RNase inhibitor, 4 μl of 10× PCR buffer, 1 μl of 10 mM dNTP mix, 1.0 μl of random hexamers, 1.0 μl of reverse transcriptase (Invitrogen). Mixture was incubated for 50 min at 42°C and the reaction was stopped by heating to 70°C for 15 min. For each of the Genes studied and 18S, 3 μl of 400 ng/mL diluted cDNA was added to 12.5 μl of Taqman Universal Master Mix (Applied Biosystems), 1.25 μl of 20× primer/Probe combination, and 8.25 μl of DEPC-treated water per each reaction. RT-PCR quantification and determination of expression levels were performed on ABI Prism 7900 and Sequence Detection Software 1.6 (Applied Biosystems, Foster City, CA, USA). All PCR reactions were performed in triplicates in our Gene Expression Core Facility and the results were analyzed using comparative method following normalization of expression values to 18S rRNA expression.

The Effects of ALDH Inhibitors on ALDH Activity, Cell Growth and Gene Expression

For comparison to the siRNA induced knock-down of ALDH activity and as a way of further validation of the microarray results, we also used known inhibitors of ALDH, Tetraethylthiuram disulfide (TT, disulfiram) or diethylaminobenzaldehyde (DEAB), to study the effects of ALDH inhibition on cell growth and expression of selected genes. These inhibitors were purchased from Sigma Biochemicals (Milwaukee, WI) dissolved in DMSO and used in a final DMSO concentration of ≤ 0.1%.

In these experiments, we used the A549 and H522 parent cell lines as well as H1299 cells, a cell line known not to have any measurable ALDH activity [32]. Dose response curves (0, 0.1, 1, 10, and 100 μM) were performed for each inhibitor measuring ALDH activity, cell growth and viability. We used either manual cell counts or the calorimetric (MTT) assay, as described above. In subsequent experiments only 1 and 100 μM concentrations were used. Cells were plated at 2 × 105/well in 6-well plates for 24 hr before adding the inhibitor. After 48 hr incubation with and without the inhibitors (untreated controls), the cells were harvested, counted and then used for ALDH activity assay as described above. The results were expressed as % decrease in cell proliferation or ALDH activity in comparison to untreated control.

In order to confirm that changes in the expression of selected genes are indeed specific to the siRNA induced knock-down of ALDH activity, we performed real-time RT-PCR on 3 of the 7 genes mentioned above (CXXC5, LHX8, RAB20) using RNA obtained from WT A549 cells after incubation with 5 μM of either TT or DEAB for 48 hr. Obtaining similar changes in the expression of these 3 genes after treatment with ALDH inhibitors would strongly suggest direct relationship between ALDH1A1 and ALDH3A1 activity and the changes in the gene profile demonstrated by the microarray analysis.

Results

ALDH Knock-Down Effect on ALDH Activity and Proteins

The expression of the specific siRNAs in A549 cells resulted in a significant decrease in ALDH activity as demonstrated by spectrophotometric enzyme activity assay (Figure 2A). The effect of siRNA against each ALDH isozyme was also demonstrated by Western blot analysis (Figure 2A). With such significant decrease in ALDH activity and proteins in the Lenti 1+3 cell line, we proceeded with the microarray gene profiling experiment comparing gene expression in this cell line to that in GFP and WT control cell lines. Similar decrease in ALDH activity was obtained in H522 lung cancer cell line using the same lentiviral constructs (Figure 2B). However, we have observed the loss of siRNA effect after continuous culture of either cell line for more than 4 weeks. Because of that, fresh cells were thawed every two weeks, and ALDH activity verified before any experiments were performed.

Figure 2
figure 2

A. The effect of siRNA against ALDH1A1 and ALDH3A1 on ALDH activity and protein. The mean percent ± SD of remaining ALDH activity (in comparison to control GFP group) as measured by spectrophotometry enzyme activity assay is shown in A549 cells transduced by lentiviral vector containing siRNA against ALDH1A1 (Lenti 1), ALDH3A1 (Lenti 3), both vectors (Lenti 1+3), green fluorescent protein (GFP), and untreated parent cell line (WT). The results represent at least 3 different experiments. Furthermore, Western blot analysis shows the decrease in ALDH1A1 and ALDH3A1 proteins in Lenti 1+3 group in comparison to WT and GFP control groups. Actin is shown as a loading control. The results are typical of an experiment performed at least twice. B. Lentivirally mediated expression of siRNA constructs (shown in A.) results in similar effects on ALDH activity in H522 cell line.

ALDH Knock-Down Effects on Cell Growth and Motility

A549 and H522 cells (2 × 105/group) were plated for 72 hours and then counted and viability determined. Figure 3A shows the results which reflect mean ± SD of cell number from at least 4 experiments. The growth rate was significantly lower (P < 0.024) in Lenti-1, Lenti-3, and Lenti-1+3 in comparison to GFP, except for A549 Lenti-1 with P value of 0.052. No significant difference was found between the GFP and the WT. Using the MTT cell proliferation assay, we were able to show similar effects (data not shown).

Figure 3
figure 3

The effect of siRNA expression on cell proliferation and motility. A. Knock-down of ALDH1A1, ALDH3A1 or both results in significant decrease in cell proliferation in comparison to WT and GFP controls from A549 and H522 cells (P < 0.024, except for A549 Lenti-1 with P value of 0.052). The horizontal line indicates the starting cell number (2 × 105/group) incubated, while the bars represent mean ± SD of viable cell counts after 72 hr incubation obtained from at least 2 experiments. Scratch wound migration assay results are shown. B. The bars represent mean ± SD of 4 measurements of wound healing in the 5 different cell lines of A549 cell. Significant decrease in cell migration into scratch wound is shown after 20 hr in cells with ALDH1A1 (Lenti-1), ALDH3A1 (Lenti-3) and both isozymes (Lenti 1+3) knock-down in comparison to WT cells or cells transduced with GFP siRNA (P < 0.005). C. An example of the delayed wound healing seen with the knock-down of both isozymes (Lenti 1+3 cells) on wound healing in comparison to WT and GFP cells.

A549 cells (200/petri dish in triplicates per each group) were plated for 5 days and colonies adhered to the bottom of the plates were counted after this time period. Results show that the mean number of colonies decreased significantly (range, 27–71%) for the Lenti-1, Lenti-3 and Lenti-1+3 cells in comparison to GFP and WT cells of both A549 and H522 (P < 0.016 and 0.003, respectively).

An example of the scratch wound migration test results for A549 cells are shown in Figure 3C which demonstrates delayed wound healing of the Lenti 1+3 cells in comparison to control cells. Summary of at least 4 readings are given for all the five A549 cell lines (Figure 3B). These results show significant inhibition (P < 0.005) of cell migration by knock-down of each enzyme separately as well as both enzymes (Lenti 1+3).

Significant Changes in Gene Expression Related to "Knock-down" of ALDH Isozymes

To construct the transcriptional profiles associated with the "knock-down" of ALDH1A1 and ALDH3A1 in A549 cells (Lenti 1+3), two other cell lines, WT and GFP, were used to interrogate the Affymetrix U133 Plus human genome GeneChips containing 38500 well-known genes. The statistical analyses are detailed above. We performed a modified F test to identify probe sets whose hybridization signal intensities varied significantly as a function of the ALDH "Knock-down". Overall about 3735 probe sets out of total of 54675 sets were detected above background and statistically different (P < 0.001) in the Lenti 1+3 experimental group when compared to the WT group (Figure 4) (supplemental material is available on http://www.ncbi.nlm.nih.gov/geo using accession number GSE8045). Many more changes in gene expression (17269 probe sets) were seen when 3 way comparisons were done among the 3 experimental groups mainly due to non-specific changes seen in the GFP group. In order to arrange the significant genes into functional framework, we chose only those genes that were significant (P < 0.001) and differed by ≥ 2 fold in the Lenti 1+3 group when compared to the WT and GFP groups. There were 26 up-regulated and 77 down-regulated genes including ALDH1A1 and ALDH3A1. We summarized the chromosomal location for these genes in Figure 5. We also researched their biological and biochemical properties using multiple internet bioinformatics search engines, web sites, and NCBI resources including BLAST, PubMed, etc. The collective results of those searches were synthesized into color-coded, mechanism-based pies of gene expression (Figure 6). These illustrations provide the relative dominance of each of the biological-biochemical processes most affected by the "Knock-down" of the ALDH isozymes. The names of these genes are tabulated in Tables 2 and 3. In summary, the affected genes are found on multiple chromosomes with the largest cluster (8 genes) observed on chromosome 12. The biological pathway most affected was signal transduction, organ development and Transcription among both the induced and repressed genes, although more pathways were found to be involved with the repressed genes. On the other hand, the biochemical pathways most affected in induced genes were protein binding, cell adhesion and lipid metabolism, while with the suppressed genes protein binding, intracellular signaling and DNA related processes were most affected (see Figure 6). As shown in Table 2, eight of the repressed genes (ID4, DDX3Y, RPS4Y1, EIF1AY, CCL20, GPC6, GPR37, and HMGA2) had the highest decrease (> 6 folds) in expression. At least of these genes, RPS4Y1 and CCL20 are highly expressed in WT A549 cells.

Table 2 Down regulated Genes (P < 0.001 and ≥ 2 folds) in Lenti 1+3 cells in comparison to WT and GFP control cells.
Table 3 Up regulated genes (P < 0.001 and ≥ 2 folds) in Lenti 1+3 cells when compared to WT and GFP control cells.
Figure 4
figure 4

Hybridization signal intensities varied significantly as a function of the ALDH "Knock-down". Four different Affymetrix U133 Plus human genome GeneChips were interrogated using labeled cRNA from WT or Lenti 1+3 cells from 4 different established clones (Vertical columns). Overall about 3735 probe sets (horizontal rows) out of total of 54675 sets were detected above background and statistically different (P < 0.001) in the Lenti 1+3 experimental group when compared to the WT group.

Figure 5
figure 5

Chromosomal clustering of affected genes. Chromosomal clustering of 101 genes that were found to be significantly different (P < 0.001) and by ≥ 2 folds in Lenti 1+3 cells in comparison to WT and GFP control groups. The black bars represent the number of up-regulated genes, while the grey bars represent the number of down-regulated genes on each chromosome.

Figure 6
figure 6

Color-coded pies showing the distribution of induced and repressed genes into functional categories. There were 26 significantly induced genes and 75 significantly repressed genes as a result of ALDH activity knock down in Lenti 1+3 cells compared to control groups WT and GFP. Both groups of genes were categorized once into various known biochemical pathways and again into the different biological pathways.

The Overall Effect on Other ALDH Isozymes

Since the human genome contains 19 ALDH functional genes and three pseudogenes [5], we sought to analyze the effect of the siRNA sequences used here on the other members of the ALDH family. We interrogated the data set of probes that showed changes in gene expression and pulled out the ALDH genes and compared changes among the 3 different experimental groups. The data is summarized in Table 4, and it shows that the siRNA sequences used are quite specific which validates our data further. ALDH1A1 and ALDH3A1 transcripts were decreased by about 7 and 4 folds, respectively. Total of another 12 ALDH genes showed some expression variability, but it was either non-specific (similar in GFP and Lenti 1+3) or minimal.

Table 4 Changes in the expression of the different ALDH genes profiled in the microarray including ALDH1A1 and ALDH3A1.

Comparison of Apparent Gene Expression among DNA Array and Real-Time RT PCR

In order to assess the validity of the microarray results, we used real-time RT-PCR to compare the expression of 8 transcripts with significant expression changes (P < 0.001 and ≥ 2 folds in Lenti 1+3) per the microarray in all 5 cell lines mentioned in Figure 2A include WT, GFP, Lenti 1, Lenti 3 and Lenti 1+3. Figure 7 shows results for 3 genes CXXC5, RAB20, and LHX8 that confirm the results seen in the 3 experimental groups of the microarray. The increase in CXXC5 and RAB20 expression was synergistic with the knock-down of both ALDH1A1 and ALDH3A1, while the effect seen on LHX8 seems to be related to the lentiviral vector integration/effect because it is seen equally in all cell lines except WT. HMOX1 and MGST1 were increased in Lenti 1+3 and to a lesser degree with either Lenti 1 or Lenti 3. On the other hand, MAP3K8 expression was decreased in Lenti 1+3 only, while STRA6 was decreased in Lenti 1+3 and Lenti 1. Thus, the real-time RT PCR measurements of all 7 transcripts showed agreement of 100% with the microarray results. LHX8 was the only gene out of this group exhibiting non specific change in expression, both by microarray as well real-time RT PCR. On the other hand, MGST1 and STRA6 change in expression did not meet the cutoff of ≥ 2 folds.

Figure 7
figure 7

Quantitative real-time RT-PCR to validate the results of microarray analysis. Results of quantitative real-time RT-PCR for 3 (CXXC5, RAB20, LHX8) out of 7 seven genes studied are shown here. These genes were found to be significantly different (P > 0.001) by microarray analysis in Lenti 1+3 experimental group, in which both ALDH1A1 and ALDH3A1 were knocked down by siRNA expression, in comparison to WT control group. The real-time RT-PCR was performed on RNA obtained from 5 cell lines including cells that express siRNA against ALDH1A1 only (Lenti 1), against ALDH3A1 only (Lenti 3), or against green fluorescent protein (GFP) as another control, in addition to the Lenti 1+3 and WT groups. The results are expressed as the log 10 of relative change in each gene expression with the denominator being the 18S as a reference control and using the fold increase (+) or decrease (-) compared to the WT/18S ratio as the baseline expression value. The direction (up- or down-regulation) for all 7 genes confirmed the results of microarray analysis even in the Lenti 1 and Lenti 3 cells. LHX8 was the only gene out of the 7 tested that was non-specifically down-regulated since it showed similar changes in the GFP control group as well.

ALDH Inhibitors Effects on Cell Proliferation and Gene Expression

To further validate our results, cells (1000/well) were incubated with ALDH inhibitors, DEAB or TT, for 72 hrs and cell proliferation measured by the calorimetric MTT based assay (Figure 8A). Absorbance was significantly lower (P < 0.0001) after adding any concentration of TT to A549 cells as well as to H522 cells. Absorbance was significantly lower (P < 0.023) after adding 0.1 μM DEAB to A549 cells, while for H522 cells, the absorbance was significantly lower (P < 0.007) only after adding ≥ 1 μM of DEAB. The slight recovery of cell proliferation at 10 μM concentration of TT has been described before [19]. These proliferation inhibitory effects seem to be specific to the inhibition of ALDH activity since the use of these inhibitors in H1299 cells lacking ALDH1A1 and ALDH3A1 expression or any significant ALDH activity [32] resulted in significantly less effects (P < 0.029) of either TT or DEAB on H1299 cell growth in comparison to H522 and A549 cell lines (Figure 8B). The dose-response effect of both inhibitors (1 and 100 μM) on ALDH activity in A549 and H522 cells is shown in Figure 8C.

Figure 8
figure 8

The effects of ALDH inhibitors DEAB and TT on ALDH activity, cell proliferation and gene expression. A. Dose-response effect of DEAB and TT on cell proliferation of A549 and H522 lung cancer cell lines. A colorimetric non-radioactive (MTT based) assay was used as described in Methods to quantify cell proliferation. The results are expressed as the mean ± SD of absorbance and reflect 8 measurements from 2 experiments. Overall, significant decrease (P < 0.023) in cell proliferation was seen with both inhibitors. B &C. H1299 lung cancer cell line with no measurable ALDH activity and A549 and H522 cells were cultured with or without ALDH inhibitors and cell proliferation and ALDH activity were measured as described in Methods. The results are expressed as % decrease (mean ± SD of at least 2 experiments) in cell proliferation (B) or ALDH activity (C) in comparison to untreated controls. The results demonstrate that the effect of ALDH inhibitors (1 and 100 μM) on H1299 cell proliferation is significantly (P < 0.029) less than that on A549 and H522. These results suggest that the effect on cell proliferation is mostly due to the inhibition of ALDH activity. D. After 48 hr incubation of A549 cells with 5 μM of DEAB (dotted bars) or TT (clear bars), RNA was extracted and real-time RT-PCR was performed, as described in Methods, for the 3 genes shown in Figure 7. The results are expressed as the log 10 of relative change in gene expression when compared to untreated WT cells (see Figure 7 legend). These results confirm that the increase in CXXC5 and RAB20 expression is most likely related to the decrease in ALDH activity. The effect on LHX8 was different by the two inhibitors and therefore non specific to ALDH activity inhibition.

Furthermore, 3 of the genes (CXXC5, RAB20 and LHX8, shown in Figure 7) were examined using real-time RT-PCR after 48 hr incubation of A549 cells with 5 μM DEAB or TT. The results show similar changes in the expression of CXXC5 and RAB20 to what was seen in the microarray results and real-time RT-PCR (Figure 8D). Once again change in LHX8 expression was affected differently by the two different inhibitors, which is suggestive of non-specific effect.

Discussion

In this study, we have used microarray analysis to assess global changes in gene expression that result from the knock-down of ALDH1A1 and ALDH3A1 in lung cancer cell line in order to gain insight into the underlying mechanisms that lead to up regulation of these two isozymes in some cancers. The goal of this study was also to identify genes or gene networks that interact with these isozymes and may be responsible for their biological effects. We have been successful in achieving significant knock-down of these two isozymes that resulted in abolishing up to 95% of ALDH activity in A549 cells used in these experiments. We have seen significant changes in the rate of growth of H522 and A549 cell lines as well as significant decrease in cell motility and migration of A549 cells with the knock-down of ALDH1A1 and/or ALDH3A1. Similar effects have been described before with the in vitro use of low doses of disulfiram (TT) [41, 19], and we were able to show similar effects on cell growth by the addition of DEAB or TT, both are known ALDH inhibitors, to cell cultures in experiments presented above. As we mentioned in the introduction, multiple other mechanisms have been elucidated for the TT effects on cell growth and apoptosis, however our results suggest that it may have to do mostly with ALDH activity either directly or indirectly through the multiple metabolic effects of ALDH isozymes. This is mainly based on two facts that are reported here for the first time: another ALDH inhibitor such as DEAB causes similar effects on cell growth like TT and the microarray gene analysis presented here reveals that multiple genes are affected with the knock-down of ALDH1A1 and ALDH3A1.

Our results also show that the siRNA sequences used here are specific to these two isozymes and do not significantly affect 12 other ALDH genes according to the microarray analysis. Thus, we feel confident that the findings presented here are indeed true representation of the changes in ALDH1A1 and ALDH3A1 activity.

The role of ALDH isozymes in many endobiotic and xenobiotic metabolic pathways has been described in details before. ALDH1A1 and/or ALDH3A1 have been implicated in retinoic acid metabolism, metabolism of aldehydes produced during lipid metabolism, alcohol metabolism, cyclophosphamide metabolism, and oxidative stress response [5]. We have reported the upregualtion of ALDH1A1 mRNA and protein in hematopoietic progenitors by interleukin-1 and tumor necrosis factor alpha [42]. Furthermore, some mutations in ALDH isozymes have been described that lead to specific clinical phenotypes or diseases, and many of these diseases are characterized by neurologic abnormalities [5]. Disulfiram (Antabuse, TT), an ALDH inhibitor used to treat alcoholism, has been reported to cause significant neurologic toxicity in patients [43, 44]. Thus, it is not surprising to detect such wide spread changes in the levels of multiple genes as revealed in our microarray analysis. These include down and up regulated genes (≥ 2 fold) that are involved in neural biological pathways, such as GPR37, PMP22, NTS, KCTD12 (PFET1), SMAD1 (all down regulated ≥ 2 fold), as well as PAK3, CNTNAP2, and OR51E2 (up regulated ≥ 2 fold). Furthermore, several biochemical pathways are affected in the A549 cells expressing siRNA against both ALDH1A1 and ALDH3A1 such as protein binding, cell adhesion, intracellular signaling, and lipid metabolism. These pathways are inter-related since protein binding could be in the center of intracellular signaling; while lipid metabolism could lead to changes in membrane structure and, as a result, changes in adhesion properties.

The most important data we hoped to obtain from the microarray analysis was identifying genes (transcription factors) that may have a role in the regulation of expression of these two enzymes and may give clues to the mechanisms involved in their up regulation in the different types of cancer. Previous studies have shown that ALDH1A1 can be up-regulated by HOX11 [45, 46], and that peroxisome proliferator activated receptor gamma (PPAR-γ) affect the expression of ALDH3A1 via other unknown transcription factors [28]. Furthermore, the involvement of the Ah receptor nuclear translocator (ARNT) in xenobiotic induction of ALDH3A1 has been established [47]. Our results confirm the involvement of PPAR-γ in ALDH regulation which was found to be significantly up regulated by ALDH knock down (Table 3). Also in this study, we have identified significant changes in some transcription factors such as TFEC which was up regulated (≥ 2 folds), and TFAM, LISCH7, HBP1 and HOXA5 which were down regulated (≥ 2 folds). HBP1 is highly expressed in bronchial epithelium [48], while HOXA5 was reported to be abnormally expressed in NSCLC [49]. Again, these observations will open the door for further studies in order to define mechanisms involved in the regulation of ALDH1A1 and ALDH3A1.

Most importantly, our finding showing significant decrease in the levels of 8 genes by > 6 folds implies the importance of ALDH activity in cell homeostasis and cancer transformation. These include CCL20, GPR37, DDX3Y, ID4, GPC6, RPS4Y1, EIF1AY, and HMGA2. According to the microarray results, RPS4Y1 and CCL20 are highly expressed in A549 cells, while the other 5 genes are expressed at lower levels. Information about the biological functions of these 8 genes was extracted form reviewing relevant publications and is summarized in Table 5. Some of these genes are involved in transcription regulation, cell growth, differentiation and apoptosis [50–56]. Furthermore, seven of these genes (not including RPS4Y1) were implicated and studied in various cancers, but only HMGA2 was recently reported to be over expressed in lung cancer and inversely associated with survival [56]. HMGA2 is found on chromosome 12q13-15, and indeed clinically relevant chromosome 12 abnormalities have been reported among the frequent chromosomal abnormalities described in non-small cell lung cancer [57–59]. Furthermore, our data show that the largest cluster of genes affected by ALDH knock-down was on chromosome 12. Furthermore, 3 of the highly affected genes are located on Y chromosome (Table 5). Deletion of Y chromosome is one of the frequent abnormalities reported in NSCLC [60–63] and is associated with malignant transformation and development of lung cancer [64–66]. We have recently reported our immunohistochemistry results on archived pathological specimens of patients with lung cancer which showed that both isozymes are highly expressed in non-small cell primary lung cancers [67]. Overall, these findings indicate that some of the genes identified in this study are good candidates for future investigation aiming at defining gene networks in which ALDH isozymes play an important role in cancer biology and may, in turn, enhance our knowledge of overall disease process.

Table 5 Genes that were highly affected by the knock-down of ALDH1A1 and ALDH3A1 in A549 lung cancer cell line.

Conclusion

We have been able to show that reduction in ALDH activity by siRNA or ALDH inhibitors may have significant effects on cell growth and proliferation possibly through effects on wide spectrum of genes with different biological roles in lung cancer cells. We predict that these results will have significant implications on defining the biological significance and mechanisms involved in high expression of ALDH isozymes in NSCLC, and therefore, impact the treatment of lung cancer.

Abbreviations

ALDH:

aldehdye dehydrogenase

ALDH1A1 and ALDH3A1:

aldehyde dehydrogenase class 1A1 and 3A1

DEAB:

diethylaminobezaldehyde

TT:

tetraethylthiuram disulfide or disulfiram

NSCLC:

non small cell lung cancer

Lenti 1 or Lenti 3:

cells expressing siRNA against ALDH1A1 or ALDH3A1 using lentivirus

Lenti 1+3:

cells expressing siRNA against ALDH1A1 and ALDH3A1

WT:

refers to wild type cells

GFP:

refers to cells expressing siRNA against the green fluorescence protein.

References

  1. Sladek NE: Human aldehyde dehydrogenases: potential pathological, pharmacological, and toxicological impact. J Biochem Mol Toxicol. 2003, 17: 7-23. 10.1002/jbt.10057

    Article  CAS  PubMed  Google Scholar 

  2. Vasiliou V, Pappa A, Estey T: Role of human aldehyde dehydrogenases in endobiotic and xenobiotic metabolism. Drug Metab Rev. 2004, 36: 279-299. 10.1081/DMR-120034001

    Article  CAS  PubMed  Google Scholar 

  3. Vasiliou V, Bairoch A, Tipton KF, Nebert DW: Eukaryotic aldehyde dehydrogenase (ALDH) genes: human polymorphisms, and recommended nomenclature based on divergent evolution and chromosomal mapping. Pharmacogenetics. 1999, 9: 421-434. 10.1097/00008571-199910000-00004

    Article  CAS  PubMed  Google Scholar 

  4. Vasiliou V, Nebert DW: Analysis and update of the human aldehyde dehydrogenase (ALDH) gene family. Hum Genomics. 2005, 2: 138-143.

    PubMed Central  CAS  PubMed  Google Scholar 

  5. Vasiliou V, Pappa A: Polymorphisms of human aldehyde dehydrogenases. Consequences for drug metabolism and disease. Pharmacology. 2000, 6: 192-198. 10.1159/000028400.

    Article  Google Scholar 

  6. Niederreither K, McCaffery P, Drager UC, Chambon P, Dolle P: Restricted expression and retinoic acid-induced downregulation of the retinaldehyde dehydrogenase type 2 (RALDH-2) gene during mouse development. Mech Dev. 1997, 62: 67-78. 10.1016/S0925-4773(96)00653-3

    Article  CAS  PubMed  Google Scholar 

  7. Niederreither K, Subbarayan V, Dolle P, Chambon P: Embryonic retinoic acid synthesis is essential for early mouse post-implantation development. Nat Genet. 1999, 21: 444-448. 10.1038/7788

    Article  CAS  PubMed  Google Scholar 

  8. Dupe V, Matt N, Garnier JM, Chambon P, Mark M, Ghyselinck NB: A newborn lethal defect due to inactivation of retinaldehyde dehydrogenase type 3 is prevented by maternal retinoic acid treatment. Proc Natl Acad Sci USA. 2003, 100: 14036-14041. 10.1073/pnas.2336223100

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Magni M, Shammah S, Schiro R, Mellado W, Dalla-Favera R, Gianni AM: Induction of cyclophosphamide-resistance by aldehyde-dehydrogenase gene transfer. Blood. 1996, 87: 1097-1103.

    CAS  PubMed  Google Scholar 

  10. Bunting KD, Lindahl R, Townsend AJ: Oxazaphosphorine-specific resistance in human MCF-7 breast carcinoma cell lines expressing transfected rat class 3 aldehyde dehydrogenase. J Biol Chem. 1994, 269: 23197-23203.

    CAS  PubMed  Google Scholar 

  11. Bunting KD, Townsend AJ: De novo expression of transfected human class 1 aldehyde dehydrogenase (ALDH) causes resistance to oxazaphosphorine anti-cancer alkylating agents in hamster V79 cell lines. Elevated class 1 ALDH activity is closely correlated with reduction in DNA interstrand cross-linking and lethality. J Biol Chem. 1996, 271: 11884-11890. 10.1074/jbc.271.17.10237

    Article  CAS  PubMed  Google Scholar 

  12. Sladek NE: Aldehyde dehydrogenase-mediated cellular relative insensitivity to the oxazaphosphorines. Curr Pharm Des. 1999, 5: 607-625.

    CAS  PubMed  Google Scholar 

  13. Moreb J, Schweder M, Suresh A, Zucali JR: Overexpression of the human aldehyde dehydrogenase class I results in increased resistance to 4-hydroperoxycyclophosphamide. Cancer Gene Ther. 1996, 3: 24-30.

    CAS  PubMed  Google Scholar 

  14. Moreb JS, Schweder M, Gray B, Zucali J, Zori R: In vitro selection for K562 cells with higher retrovirally mediated copy number of aldehyde dehydrogenase class-1 and higher resistance to 4-hydroperoxycyclophosphamide. Hum Gene Ther. 1998, 9: 611-619. 10.1089/hum.1998.9.5-611

    Article  CAS  PubMed  Google Scholar 

  15. Pappa A, Brown D, Koutalos Y, DeGregori J, White C, Vasiliou V: Human aldehyde dehydrogenase 3A1 inhibits proliferation and promotes survival of human corneal epithelial cells. J Biol Chem. 2005, 280: 27998-28006. 10.1074/jbc.M503698200

    Article  CAS  PubMed  Google Scholar 

  16. Moreb JS, Maccow C, Schweder M, Hecomovich J: Expression of antisense RNA to aldehyde dehydrogenase class-1 sensitizes tumor cells to 4-hydroperoxycyclophosphamide in vitro. J Pharmacol Exp Ther. 2000, 293: 390-396.

    CAS  PubMed  Google Scholar 

  17. Moreb JS, Gabr A, Vartikar GR, Gowda S, Zucali JR, Mohuczy D: Retinoic acid down-regulates aldehyde dehydrogenase and increases cytotoxicity of 4-hydroperoxycyclophosphamide and acetaldehyde. J Pharmacol Exp The. 2005, 312: 339-345. 10.1124/jpet.104.072496.

    Article  CAS  Google Scholar 

  18. Moreb JS, Muhoczy D, Ostmark B, Zucali JR: RNAi-mediated knockdown of aldehyde dehydrogenase class-1A1 and class-3A1 is specific and reveals that each contributes equally to the resistance against 4-hydroperoxycyclophosphamide. Cancer Chemother Pharmaco. 2006, 59: 127-136. 10.1007/s00280-006-0233-6.

    Article  Google Scholar 

  19. Wickström M, Danielsson K, Rickardson L, Gullbo J, Nygren P, Isaksson A, Larsson R, Lövborg H: Pharmacological profiling of disulfiram using human tumor cell lines and human tumor cells from patients. Biochem Pharmacol. 2007, 73: 25-33. 10.1016/j.bcp.2006.08.016

    Article  PubMed  Google Scholar 

  20. Chen D, Cui QC, Yang H, Dou QP: Disulfiram, a clinically used anti-alcoholism drug and copper-binding agent, induces apoptotic cell death in breast cancer cultures and xenografts via inhibition of the proteasome activity. Cancer Res. 2006, 66: 10425-10433. 10.1158/0008-5472.CAN-06-2126

    Article  CAS  PubMed  Google Scholar 

  21. Cen D, Brayton D, Shahandeh B, Meyskens FL, Farmer PJ: Disulfiram facilitates intracellular Cu uptake and induces apoptosis in human melanoma cells. J Med Chem. 2004, 47: 6914-6920. 10.1021/jm049568z

    Article  CAS  PubMed  Google Scholar 

  22. Brar SS, Grigg C, Wilson KS, Holder WD, Dreau D, Austin C, Foster M, Ghio AJ, Whorton AR, Stowell GW, Whittall LB, Whittle RR, White DP, Kennedy TP: Disulfiram inhibits activating transcription factor/cyclic AMP-responsive element binding protein and human melanoma growth in a metal-dependent manner in vitro, in mice and in a patient with metastatic disease. Mol Cancer Ther. 2004, 3: 1049-1060.

    CAS  PubMed  Google Scholar 

  23. Wang W, McLeod HL, Cassidy J: Disulfiram-mediated inhibition of NF-kappaB activity enhances cytotoxicity of 5-fluorouracil in human colorectal cancer cell lines. Int J Cancer. 2003, 104: 504-511.10.1002/ijc.10972

    Article  CAS  PubMed  Google Scholar 

  24. Liu GY, Frank N, Bartsch H, Lin JK: Induction of apoptosis by thiuramdisulfides, the reactive metabolites of dithiocarbamates, through coordinative modulation of NFkappaB, c-fos/c-jun, and p53 proteins. Mol Carcinog. 1998, 22: 235-246. 10.1002/(SICI)1098-2744(199808)22:4<235::AID-MC5>3.0.CO;2-I

    Article  CAS  PubMed  Google Scholar 

  25. Kim CH, Kim JH, Moon SJ, Hsu CY, Seo JT, Ahn YS: Biphasic effects of dithiocarbamates on the activity of nuclear factor-kappaB. Eur J Pharmacol. 2000, 392: 133-136. 10.1016/S0014-2999(00)00109-6

    Article  CAS  PubMed  Google Scholar 

  26. Yakisich JS, Sidén A, Eneroth P, Cruz M: Disulfiram is a potent in vitro inhibitor of DNA topoisomerases. Biochem Biophys Res Commun. 2001, 289: 586-590. 10.1006/bbrc.2001.6027

    Article  CAS  PubMed  Google Scholar 

  27. Nobel CS, Kimland M, Nicholson DW, Orrenius S, Slater AF: Disulfiram is a potent inhibitor of proteases of the caspase family. Chem Res Toxicol. 1997, 10: 1319-1324. 10.1021/tx970131m

    Article  CAS  PubMed  Google Scholar 

  28. Canuto RA, Maggiora M, Trombetta A, Martinasso G, Muzio G: Aldehyde dehydrogenase 3 expression is decreased by clofibrate via PPAR gamma induction in JM2 rat hepatoma cell line. Chem Biol Interact. 2003, 143–144: 29-35. 10.1016/S0009-2797(02)00169-2.

    Article  PubMed  Google Scholar 

  29. Prough RA, Falkner KC, Xiao GH, Lindahl RG: Regulation of rat ALDH-3 by hepatic protein kinases and glucocorticoids. Adv Exp Med Biol. 1997, 414: 29-36.

    Article  CAS  PubMed  Google Scholar 

  30. Vasiliou V, Puga A, Nebert DW: Negative regulation of the murine cytosolic aldehyde dehydrogenase-3 (Aldh-3c) gene by functional CYP1A1 and CYP1A2 proteins. Biochem Biophys Res Commun. 1992, 187: 413-419. 10.1016/S0006-291X(05)81508-6

    Article  CAS  PubMed  Google Scholar 

  31. Xiao G, Falkner KC, Xie Y, Lindahl RG, Prough RA: cAMP-dependent negative regulation of rat aldehyde dehydrogenase class 3 gene expression. J Biol Chem. 1997, 272: 3238-3245. 10.1074/jbc.272.6.3238

    Article  CAS  PubMed  Google Scholar 

  32. Moreb JS, Zucali JR, Ostmark B, Benson NA: Heterogeneity of aldehyde dehydrogenase expression in lung cancer cell lines is revealed by Aldefluor flow cytometry-based assay. Cytometry B Clin Cytom. 2007, 72: 281-289.

    Article  PubMed  Google Scholar 

  33. Suyama E, Kawasaki H, Nakajima M, Taira K: Identification of genes involved in cell invasion by using a library of randomized hybrid ribozymes. Proc Natl Acad Sci USA. 2003, 100: 5616-5621. 10.1073/pnas.1035850100

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Xiong C, Tang DQ, Xie CQ, Zhang L, Xu KF, Thompson WE, Chou W, Gibbons GH, Chang LJ, Yang LJ, Chen YE: Genetic engineering of human embryonic stem cells with lentiviral vectors. Stem Cells Dev. 2005, 14: 367-377. 10.1089/scd.2005.14.367

    Article  PubMed  Google Scholar 

  35. Zaiss A-K, Son S, Chang L-J: RNA 3' readthrough of oncoretrovirus and lentivirus: implications for vector safety and efficacy. J Virol. 2002, 76: 7209-7219. 10.1128/JVI.76.14.7209-7219.2002

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Dockham PA, Sreerama L, Sladek NE: Relative cotribution of human erythrocyte aldehyde dehydrogenase to the systemic detoxification of the oxazaphosphorines. Drug Metab Dispos. 1997, 25: 1436-1441.

    CAS  PubMed  Google Scholar 

  37. Sreerama L, Sladek NE: Class 1 and class 3 aldehyde dehydrogenase levels in the human tumor cell lines currently used by the National Cancer Institute to screen for potentially useful antitumor agents. Adv Exp Med Biol. 1997, 414: 81-94.

    Article  CAS  PubMed  Google Scholar 

  38. Li C, Hung Wong W: Model-based analysis of oligonucleotide arrays: model validation, design issues and standard error application. Genome Biol. 2001, 2: 1-0032.

    Google Scholar 

  39. Velardo MJ, Burger C, Williams PR, Baker HV, López MC, Mareci TH, White TE, Muzyczka N, Reier PJ: Patterns of gene expression reveal a temporally orchestrated wound healing response in the injured spinal cord. J Neurosci. 2004, 24: 8562-8576. 10.1523/JNEUROSCI.3316-04.2004

    Article  CAS  PubMed  Google Scholar 

  40. Lee LG, Connell CR, Bloch W: Allelic discrimination by nick-translation PCR with fluorogenic probes. Nucleic Acids Res. 1993, 21: 3761-3766. 10.1093/nar/21.16.3761

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Brar SS, Grigg C, Wilson KS, Holder WD, Dreau D, Austin C, Foster M, Ghio AJ, Whorton AR, Stowell GW, Whittall LB, Whittle RR, White DP, Kennedy TP: Disulfiram inhibits activating transcription factor/cyclic AMP-responsive element binding protein and human melanoma growth in a metal-dependent manner in vitro, in mice and in a patient with metastatic disease. Mol Cancer Ther. 2004, 3: 1049-1060.

    CAS  PubMed  Google Scholar 

  42. Moreb JS, Turner C, Sreerama L, Zucali JR, Sladek NE, Schweder M: Interleukin-1 and tumor necrosis factor alpha induce class 1 aldehyde dehydrogenase mRNA and protein in bone marrow cells. Leuk Lymphoma. 1995, 20: 77-84. 10.3109/10428199509054756

    Article  CAS  PubMed  Google Scholar 

  43. Chick J: Safety issues concerning the use of disulfiram in treating alcohol dependence. Drug Saf. 1999, 20: 427-435. 10.2165/00002018-199920050-00003

    Article  CAS  PubMed  Google Scholar 

  44. Mirsal H, Yalug I, Tan D, Stern TA, Kalyoncu A, Pektas O, Erdogan G, Beyazyürek M: Delirium-associated disulfiram and ethanol interactions. Prim Care Companion J Clin Psychiatry. 2005, 7: 235-237.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Greene WK, Bahn S, Masson N, Rabbitts M: The T-cell oncogenic protein HOX11 activates Aldh1 expression in NIH 3T3 cells but represses its expression in mouse spleen development. Mol Cell Biol. 1998, 18: 7030-7037.

    PubMed Central  CAS  PubMed  Google Scholar 

  46. Owens BM, Zhu YX, Suen TC, Wang PX, Greenblatt JF, Goss PE, Hawley RG: Specific homeodomain-DNA interactions are required for HOX11-mediated transformation. Blood. 2003, 101: 4966-4974. 10.1182/blood-2002-09-2857

    Article  CAS  PubMed  Google Scholar 

  47. Reisdorph R, Lindahl R: Aldehyde dehydrogenase 3 gene regulation: studies on constitutive and hypoxia-modulated expression. Chem Biol Interact. 2001, 130–132: 227-233. 10.1016/S0009-2797(00)00267-2.

    Article  PubMed  Google Scholar 

  48. Ioka RX, Kang MJ, Kamiyama S, Kim DH, Magoori K, Kamataki A, Ito Y, Takei YA, Sasaki M, Suzuki T, Sasano H, Takahashi S, Sakai J, Fujino T, Yamamoto TT: Expression cloning and characterization of a novel glycosylphosphatidylinositol-anchored high density lipoprotein-binding protein, GPI-HBP1. J Biol Chem. 2003, 278: 7344-7349. 10.1074/jbc.M211932200

    Article  CAS  PubMed  Google Scholar 

  49. Abe M, Hamada J, Takahashi O, Takahashi Y, Tada M, Miyamoto M, Morikawa T, Kondo S, Moriuchi T: Disordered expression of HOX genes in human non-small cell lung cancer. Oncol Rep. 2006, 15: 797-802.

    CAS  PubMed  Google Scholar 

  50. Abdelhaleem M: RNA helicases: regulators of differentiation. Clin Biochem. 2005, 38: 499-503. 10.1016/j.clinbiochem.2005.01.010

    Article  CAS  PubMed  Google Scholar 

  51. Filmus J: Glypicans in growth control and cancer. Glycobiology. 2001, 11: 19R-23R. 10.1093/glycob/11.3.19R

    Article  CAS  PubMed  Google Scholar 

  52. Kleiman SE, Yogev L, Hauser R, Botchan A, Maymon BB, Paz G, Yavetz H: Expression profile of AZF genes in testicular biopsies of azoospermic men. Hum Reprod. 2007, 22: 151-158.

    Article  CAS  PubMed  Google Scholar 

  53. Peddada S, Yasui DH, LaSalle JM: Inhibitors of differentiation (ID1, ID2, ID3 and ID4) genes are neuronal targets of MeCP2 that are elevated in Rett syndrome. Hum Mol Genet. 2006, 15: 2003-2014. 10.1093/hmg/ddl124

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Rezgaoui M, Susens U, Ignatov A, Gelderblom M, Glassmeier G, Franke I, Urny J, Imai Y, Takahashi R, Schaller HC: The neuropeptide head activator is a high-affinity ligand for the orphan G-protein-coupled receptor GPR37. J Cell Sci. 2006, 119: 542-549. 10.1242/jcs.02766

    Article  CAS  PubMed  Google Scholar 

  55. Ruzinova MB, Benezra R: Id proteins in development, cell cycle and cancer. Trends Cell Biol. 2003, 13: 410-418. 10.1016/S0962-8924(03)00147-8

    Article  CAS  PubMed  Google Scholar 

  56. Sarhadi VK, Wikman H, Salmenkivi K, Kuosma E, Sioris T, Salo J, Karjalainen A, Knuutila S, Anttila S: Increased expression of high mobility group A proteins in lung cancer. J Pathol. 2006, 209: 206-212. 10.1002/path.1960

    Article  CAS  PubMed  Google Scholar 

  57. Grepmeier U, Dietmaier W, Merk J, Wild PJ, Obermann EC, Pfeifer M, Hofstaedter F, Hartmann A, Woenckhaus M: Deletions at chromosome 2q and 12p are early and frequent molecular alterations in bronchial pithelium and NSCLC of long-term smokers. Int J Oncol. 2005, 27: 481-488.

    CAS  PubMed  Google Scholar 

  58. Hoglund M, Gisselsson D, Hansen GB, Mitelman F: Statistical dissection of cytogenetic patterns in lung cancer reveals multiple modes of karyotypic evolution independent of histological classification. Cancer Genet Cytogenet. 2004, 154: 99-109. 10.1016/j.cancergencyto.2004.01.030

    Article  PubMed  Google Scholar 

  59. Feder M, Siegfried JM, Balshem A, Litwin S, Keller SM, Liu Z, Testa JR: Clinical relevance of chromosome abnormalities in non-small cell lung cancer. Cancer Genet Cytogenet. 1998, 102: 25-31. 10.1016/S0165-4608(97)00274-4

    Article  CAS  PubMed  Google Scholar 

  60. Johansson M, Jin Y, Mandahl N, Hambraeus G, Johansson L, Mitelman F, Heim S: Cytogenetic analysis of short-term cultured squamous cell carcinomas of the lung. Cancer Genet Cytogenet. 1995, 81: 46-55. 10.1016/0165-4608(94)00207-X

    Article  CAS  PubMed  Google Scholar 

  61. Dave BJ, Hopwood VL, Spitz MR, Pathak S: Sex chromosome abnormalities in lung cancer patients. Cancer Genet Cytogenet. 1996, 87: 24-28. 10.1016/0165-4608(95)00245-6

    Article  CAS  PubMed  Google Scholar 

  62. Berker-Karaüzüm S, Lüleci G, Ozbilim G, Erdoǧan A, Kuzucu A, Demircan A: Cytogenetic findings in thirty lung carcinoma patients. Cancer Genet Cytogenet. 1998, 100: 114-123. 10.1016/S0165-4608(96)00422-0

    Article  PubMed  Google Scholar 

  63. Berrieman HK, Ashman JN, Cowen ME, Greenman J, Lind MJ, Cawkwell L: Chromosomal analysis of non-small-cell lung cancer by multicolour fluorescent in situ hybridisation. Br J Cancer. 2004, 90: 900-905. 10.1038/sj.bjc.6601569

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  64. Weaver DA, Hei TK, Hukku B, McRaven JA, Willey JC: Cytogenetic and molecular genetic analysis of tumorigenic human bronchial epithelial cells induced by radon alpha particles. Carcinogenesis. 1997, 18: 1251-1257. 10.1093/carcin/18.6.1251

    Article  CAS  PubMed  Google Scholar 

  65. Weaver DA, Hei TK, Hukku B, Demuth JP, Crawford EL, McRaven JA, Girgis S, Willey JC: Localization of tumor suppressor gene candidates by cytogenetic and short tandem repeat analyses in tumorigenic human bronchial epithelial cells. Carcinogenesis. 2000, 21: 205-211. 10.1093/carcin/21.2.205

    Article  CAS  PubMed  Google Scholar 

  66. Kayser K, Kosjerina Z, Goldmann T, Kayser G, Kazmierczak B, Vollmer E: Lung carcinoma-associated atypical adenomatoid hyperplasia, squamous cell dysplasia, and chromosome alterations in non-neoplastic bronchial mucosa. Lung Cancer. 2005, 47: 205-214. 10.1016/j.lungcan.2004.07.042

    Article  PubMed  Google Scholar 

  67. Patel M, Lu L, Zander DS, Sreerama L, Coco D, Moreb JS: ALDH1A1 and ALDH3A1 expression in lung cancers: correlation with histologic type and potential precursors. Lung Cancer. 2008, 59: 340-349. 10.1016/j.lungcan.2007.08.033

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Financial support for this work was provided by a grant (to JSM) from the Flight Attendant Medical Research Institute (Miami, FL).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jan S Moreb.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JSM was behind the hypothesis tested here, designed and analyzed the studies, his laboratory (MA and BO) performed the cell transduction, RNA extraction and cell culture experiments; WC prepared the lentiviral constructs with the specific siRNA in Dr Chang's Laboratory; HVB and his laboratory technician (MCL) performed the microarray hypridization and analysis. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Moreb, J.S., Baker, H.V., Chang, LJ. et al. ALDH isozymes downregulation affects cell growth, cell motility and gene expression in lung cancer cells. Mol Cancer 7, 87 (2008). https://doi.org/10.1186/1476-4598-7-87

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-4598-7-87

Keywords