Skip to main content

Global gene expression profiling in human lung cells exposed to cobalt



It has been estimated that more than 1 million workers in the United States are exposed to cobalt. Occupational exposure to 59 Co occurs mainly via inhalation and leads to various lung diseases. Cobalt is classified by the IARC as a possible human carcinogen (group 2B). Although there is evidence for in vivo and in vitro toxicity, the mechanisms of cobalt-induced lung toxicity are not fully known. The purpose of this work was to identify potential signatures of acute cobalt exposure using a toxicogenomic approach. Data analysis focused on some cellular processes and protein targets that are thought to be relevant for carcinogenesis, transport and biomarker research.


A time course transcriptome analysis was performed on A549 human pulmonary cells, leading to the identification of 85 genes which are repressed or induced in response to soluble 59 Co. A group of 29 of these genes, representing the main biological functions, was assessed by quantitative RT-PCR. The expression profiles of six of them were then tested by quantitative RT-PCR in a time-dependent manner and three modulations were confirmed by Western blotting. The 85 modulated genes include potential cobalt carriers (FBXL2, ZNT1, SLC12A5), tumor suppressors or transcription factors (MAZ, DLG1, MYC, AXL) and genes linked to the stress response (UBC, HSPCB, BNIP3L). We also identified nine genes coding for secreted proteins as candidates for biomarker research. Of those, TIMP 2 was found to be down-regulated and this modulation was confirmed, in a dose-dependent manner, at protein level in the supernatant of exposed cells.


Most of these genes have never been described as related to cobalt stress and provide original hypotheses for further study of the effects of this metal ion on human lung epithelial cells. A putative biomarker of cobalt toxicity was identified.


In the United States, more than a million workers are potentially exposed to cobalt or its compounds [1]. Cobalt is massively used in the steel industry, being a major constituent of hard metal alloys, in combination with tungsten carbides. Other industrial uses include diamond polishing with Co-containing disks and the production of drying agents, pigments, and catalysts [2]. Radioactive isotopes of cobalt are used in industry, medicine and nuclear research. In nuclear power plants, 59Co-containing alloys can be activated into radioactive 60Co oxides, dispersed in the cooling water and then contaminate workers [3, 4]. A study measuring the ambient air in cobalt powder production reported concentrations of cobalt ranging from 0.675 to 10 mg/m3 [5]. Airborne concentrations measured in the working environment from a factory producing hard-metal inserts ranged from 14.6 to 37.4 mg/m3 [6]. Occupational exposure to Co occurs mainly via inhalation leading to various lung diseases, such as pneumonitis, fibrosis and asthma [7, 8]. As with human exposure, animal exposure to cobalt-containing aerosols causes pronounced respiratory effects. A single 30-minute exposure of rats to relatively high levels (26–236 mg cobalt hydrocarbonyl/m3), resulted in lung congestion, oedema, and haemorrhage [9]. Necrosis and inflammation of the respiratory tract epithelium were reported in rats exposed to 19 mg cobalt/m3 and mice exposed to 1.9 mg cobalt sulfate/m3 over 16 days [10, 11]. Some acute effects have been observed concerning general public exposure. Lethal cardiomyopathy was reported in people who consumed large quantities of beer containing cobalt as a foam stabilizer (0.04–0.14 mg cobalt/kg/day), and acute mortality accounted for 18% of the deaths [12]. A 19-month-old boy who swallowed an unknown amount of cobalt chloride solution died 6.5 hours after ingestion [13].

Following absorbtion by inhalation, cobalt is eliminated in the urine. Biological monitoring of accidental exposure mainly involves measuring the concentration of metal in the urine. This might be inadequate for several reasons. Firstly, the quantity of metal excreted (exposure marker) does not necessarily reflect organ damage, which varies from one person to another. Secondly, depending on the chemical form, excretion does not necessarily reflect the level of metal in the body because some forms are retained in the lungs. Thirdly, depending on its solubility, clearance can be very rapid and the cobalt may have left the body by the time samples can be taken. Therefore, a key issue in monitoring occupational exposure is the availability of adequate biomarkers.

Although the chemical toxicity of cobalt has been proven, the molecular mechanisms of its toxicity are not all known. Cobalt is genotoxic [14, 15], and an oxidizing stress inducer [16]. It also induces apoptosis [17]. Cobalt is used as a hypoxia-simulating agent [18], leading to increased apoptosis, glycolysis, angiogenesis and erythropoiesis[19].

Since the lung is the main target organ of cobalt toxicity, the human A549 lung cell line was chosen as a model for this study, to evaluate cobalt toxicity. Noteworthy, this cell line has been widely documented in molecular toxicology, including hypoxia mimicked by cobalt [16, 20].

Microarrays are currently used for large scale gene profiling, measuring sensitive cell changes in response to xenobiotic exposure. Such investigative studies may help identify new molecular targets for toxicants or provide new hypotheses about their mechanisms of action [21]. We used toxicogenomic tools to detect biomolecular targets of acute cobalt exposure and identify candidates as biomarkers of cobalt toxicity.

Results and discussion

Cobalt cytotoxicity

A549 is a stable tumor cell line, obtained from human lung carcinoma, with properties of type II alveolar epithelial cells [22]. The response of A549 cells to increasing concentrations of cobalt (CoCl2) after 24 h exposure was first analysed using the intracellular ATP measurement that reflects early metabolic modification [23]. We determined that the concentration needed to decrease ATP concentration to 50% was reached at 2 mM cobalt (figure 1). At this concentration, we measured an average load of 2.4 pg of cobalt per cell, using flame atomic absorption spectroscopy (data not shown). Some toxicologists, including one reviewer of this paper, assume that cytotoxic doses are too high to study specific cellular response to a toxicant [31]. The consequence, on gene modulation, could be the activation of large numbers of nonspecific pathways of toxicity, particularly a predominance of apoptosis and stress related genes. However, as the purpose of this study was to simulate accidental acute exposure, we nevertheless chose to expose A549 cells to an acute dose of cobalt (2 mM).

Figure 1
figure 1

ATP viability test on A549 cells exposed to cobalt. The cells were cultured for 24 h with increasing concentrations of cobalt. Cell viability was determined by ATP Celltiter-Glow Luminescence Cell Viability Assay (Promega). Measurements were made on the Lumistar Galaxy (BMG) luminometer. The viability (V) was determined as the ATP ratio of treated cells versus control cells in %. The results are presented as mean ± SD (n = 4 to 8).

Differential gene expression and functional classification

The global transcriptional response was monitored using CEA microarrays (GPL4263) [24]. A549 cells were treated with 2 mM cobalt in a time course experiment. For each experiment, control cells were grown in parallel and collected simultaneously. At 30 min. and 2 h, six microarrays (three dyeswaps) were hybridized with RNA from two different cell exposures. At 4 h, four microarrays (two dyeswaps) were hybridized with RNA from two exposures and at 24 h, 14 microarrays (seven dyeswaps) were hybridized with RNA from three cell exposures. Details of the experimental design and microarray data were submitted to the Gene Expression Omnibus repository at the National Center for Biotechnology Information [24] and are accessible via number GSE5892. Data analysis was performed as described in the methods section. After applying our selection criteria to the data (ratio > 1.5 and p < 0.05), we obtained a list of 173 modulated clones corresponding to 85 known genes. The modulation trend was up-regulation since 53 genes were up-regulated and 32 were down-regulated (table 1 and 2). The results showed that the main effect was obtained after 24 h of cobalt treatment (50 genes), and to a lesser extent after 30 min (16 genes). At 2 h and 4 h, 11 and 8 genes respectively were found to be modulated. Since 24 h was the main category in terms of modulated gene numbers, we reinforced the study for this time point by performing a total of 14 microarrays.

Table 1 Functional classification and ratios of genes up-regulated (> 1.5-fold) by exposure to cobalt across all time points
Table 2 Functional classification and ratios of genes down-regulated (< 1.5 fold) by exposure to cobalt across all time points

We used DAVID [25] web-accessible software, to obtain information from the Gene Ontology Database. We then sorted the modulated genes manually into 8 major functional classes (table 1 and 2), combining information from DAVID, SWISS-PROT and additional elements taken from the literature. The classifications of genes modulated after 30 min and 24 h of cobalt treatment were compared. At 30 min, the main category regulated was signal transduction and trafficking (38%) and half of the modulated genes coded for membrane proteins. This suggests strong involvement of membrane transport in the early response to cobalt exposure. The later response (24 h) mainly concerned protein expression and turnover (34%), metabolism and energy (20%), and finally signal transduction and trafficking (16%). This different pattern indicates that after 24 h, the cell response is established, with metabolic adaptation and strong involvement of protein translation. It can be noticed that in spite of the cytotoxic dose we used, the genes modulation trend was rather increase than decrease, indicating an active metabolic response of the cell, and also that apoptotic and stress related genes are not over represented.

Quantitative RT-PCR

We selected genes to be validated by quantitative RT-PCR (qRT-PCR) in some specific functional classes: development, differentiation and proliferation, signal transduction and trafficking, cell defense and finally gene transcription and modification. We also selected genes coding for secreted proteins as potential hits for biomarker research. Tables 1 and 2 show the results of quantitative RT-PCR validation (lines in bold) and the primers used are described in table 3. Of the 29 genes tested, most confirmed the variation observed in the microarrays (83%). Five genes appeared to be modulated in the opposite direction by qRT-PCR (AXL, IFNA4, DLG1, TNFRSF9 and SIN3A). Some discrepancies between the two techniques may be due to technical artefacts even if the primers were designed within the microarray cDNA sequence. For further analysis of gene modulation, we performed time-course quantitative RT-PCR on several genes: AXL, UBC, FBXL2 and SLC12A5 (figure 2). The results show that modulation varies over time following various profiles: either diminishing at first, then increasing up to 24 h (AXL, SLC12A5), or being induced from 2 h (UBC), or finally, early and continuous diminishing (FBXL2). For genes displaying temporal modulation, down-fluctuation quickly followed by up-regulation, as with AXL (figure 2), a slight difference in the kinetics of different exposures could explain the opposing results we found between microarray and qRT-PCR ratios, because RNAs from different exposures were used to perform the microarray experiments and qRT-PCR. These results, revealing an oscillatory gene expression profile, suggest complex regulation pathways in response to cobalt.

Table 3 List of primers used in qRT- PCR
Figure 2
figure 2

Time-course qRT-PCR analysis. The temporal expression pattern of genes regulated by cobalt was analysed by qRT-PCR in sextuplets for each time point. Mean results are given in Log2 ratios ± SD. *: significant modulation calculated by pairwise testing with p < 0.05


Cobalt is known as a "hypoxia-simulating" agent[26] and most of the literature therefore relates to this property. Hypoxia-inducible factor (HIF-1) is a transcription factor which acts as the main regulator of oxygen homeostasis and the genes it induces have mostly been described [27]. We compared this list with our results to identify modulations due to hypoxia. Strikingly, very poor overlap was observed, as only 7 of the 85 modulated genes (8 %) in our study matched HIF1 target genes. These 7 genes are: aldolase A (ALDOA), glucose transporter 1 (SLC2A1), glyceraldehyde 3-phosphate dehydrogenase (GAPD), lactate dehydrogenase A (LDHA), proapoptotic factor (BNIP3L), phosphoglycerate kinase 1 (PGK1) and transferrin receptor (TFRC), (table 1 and 2). Like HIF1 regulated genes, they were all induced at 24 h by cobalt (except TFRC which was repressed at 30 min), therefore our results agree with those already published and validate our system relative to previous reports, since they represent an appropriate internal positive control. Recently, two transcriptomic studies of hypoxia compared the action of low oxygen with that of metal ions such as cobalt or nickel on embryonic mouse fibrobasts [28] and human liver carcinoma [29]. 21% and 28% of the genes respectively were found to be modulated by both hypoxia and cobalt. This shows that the toxic effect of cobalt is not limited to hypoxia. It should also be noted that nickel had an action closer to that of low oxygen than cobalt, with 65% common modulated genes [29]. Comparing these results with ours showed that, apart from the hypoxia genes, we had no overlap. This discrepancy might be due to differences between cell models, species and conditions of exposure (one non-lethal dose of 0.1 mM cobalt for 24 h).

Stress response – apoptosis

The HSPCB gene coding for heat shock protein 90, was induced by a factor of 2.2 at 24 h (table 1). HSP90 protein plays a major role in the stress response by preventing irreversible protein aggregation. In our previous study of the toxicogenomics of uranium [30], we noted that both HSPCB gene and HSP90 protein were strongly repressed. This protein can even be thought of as a metal "sensor", since the HSPCB gene was also reported to be modulated by other metals (Ni, As, Cr, Cd) [31]. However, we did not see any change in the quantity of HSP90 in Western blotting experiments (figure 3). This proteome/transcriptome discordance can be explained by an increase in protein turnover associated with gene induction, resulting in a constant amount of protein being maintained in the cell. Another explanation may be that, as Gerner observed for apoptosis, HSP90 is aggregated and accumulates in the insoluble fraction of the cytosol and is not collected during protein extraction [32]. Two other genes described as HSP mediators are also induced by cobalt: ST13 and DNAJA2.

Figure 3
figure 3

Western blot confirmation. Western blot analyses were performed on A549 extracts following exposure to 2 mM cobalt for 24 h. Treated (Co+) and control (Co-) extracts were probed with anti-ubiquitin, anti-DLG1, or anti-HSP90 antibodies. An anti-actin antibody was used as a control.

Ubiquitin is a key element in eliminating proteins via proteasome because polyubiquitinylation is the recognition signal for protein elimination [33, 34]. The results show that ubiquitin genes were up-regulated after 24 h of toxic exposure (see table 1, UBB, ratio 1.94, UBC, ratio 5.6) while UBC gene induction was observed in the kinetics from 4 h with qRT-PCR (figure 2), showing that it is not an early event in cell response. The Western blot revealed protein ubiquitinylation, with the anti-ubiquitin antibody, shown by a smear at high molecular weight in the cobalt extract, indicating an accumulation of ubiquitinylated protein (figure 3). This demonstrates increased turnover of proteins involved in defense against stress, or the elimination of proteins partly damaged by the metal, or the start of apoptosis.

Cobalt has been described as an apoptosis inducer [17] and this is confirmed here by BAG1 repression and BNIP3L induction (table 1 and 2).

Tumor suppressors – transcription factors

The protein encoded by MAZ (Myc-associated zinc finger protein), is a transcription factor involved in both the initiation and termination of target gene transcription. It binds zinc and can act as a tumor suppressor [35]. This gene is repressed early at 2 h. Other cobalt-modulated genes code for zinc-binding proteins involved in transcription or DNA metabolism: BHC80 (BHC80 protein), NT5E (CD73) and SDCCAG33 (Teashirt homolog 1). Cobalt ions have been shown to substitute for zinc in zinc finger protein domains which control the transcription of several genes and also in zinc-DNA repair proteins, inhibiting DNA repair [14, 36, 37].

DLG1 (disk large homolog 1) has been identified as a tumor suppressor gene in drosophila, with its mutation inducing the loss of cell polarity and neoplastic tissue growth. This gene is highly conserved. The molecular mechanism which regulates cell proliferation by DLG1 is not well known but there are arguments for implicating epithelial cell polarization in this regulation [38]. The DLG1 (disk large homolog 1) gene was strongly repressed in qRT-PCR. The anti-DLG1 antibody revealed clear depletion of DLG1 protein in Co+ samples (figure 3). This confirmed the results of the qRT-PCR.

Two tumor suppressors (MAZ and DLG1) were repressed by cobalt, and in one case (DLG1), this repression was also observed at protein level.

We noted the modulation of three genes related to MYC, a key element in oncogenic response (table 1 and 2). The first, SIN3A, known as an MYC suppressor, is induced by cobalt. The two others: NDRG4, a gene negatively regulated by MYC and MAZ, a MYC transcription inhibitor, are down-regulated by cobalt. These results are not clear regarding MYC up- or down-regulation, but MYC regulation pathways are complex. These results led us to hypothesize that MYC might be modulated by cobalt, so we tested MYC temporal modulation with qRT-PCR. Interestingly, we noted significant MYC inhibition from 30 min to 24 h (figure 2) indicating that cobalt modulates regulation pathways controlling or controlled by MYC. Liu has recently shown that MYC is involved in the carcinogenic response of mouse liver to arsenic, another toxic metal [39].

The AXL gene is positively modulated. Members of the AXL/UFO family of tyrosine kinases are prone to transcriptional regulation and perform various functions including regulation of cell adhesion, migration, phagocytosis, and survival. The biological consequences of AXL activation are complex. AXL was initially identified as a transforming gene product, and AXL expression is indeed up-regulated in human tumors [40, 41].

Concerning the potential effects of cobalt on cancer development, it will be interesting to carry out further studies of the mechanisms by which cobalt induces these responses. These genes have never been linked to cobalt effects before.


Two cobalt transporters are described in yeast: COT1 and COT2 [42]. Eleven cobalt-tolerant mutants were obtained with a common COT2 mutation, suggesting that it is involved in cobalt uptake [43]. The mutants were resistant to cobalt toxicity and cobalt was not incorporated in yeast cells. The human homolog of COT2, FBXL2, was strongly repressed in our study as early as 2 h and with a ratio of 0.1 at 24 h; this is the strongest inhibition observed in time-course qRT-PCR analysis (figure 2). FBXL2 function is not very well described, but it is potentially involved in ubiquitinylation. Conklin's observations on the yeast homolog of FBXL2 [43] and our results, lead us to suggest that FBXL2 gene repression could be a cellular defense mechanism and this protein might be associated with cobalt uptake in humans.

The product of the COT1 gene is involved in the uptake of cobalt ions [44]. COT1, has a human homolog in ZNT1, coding for Zinc transporter 1 protein, probably involved in zinc transport out of the cell. ZNT1, not spotted on the microarray, was also tested in kinetics with qRT-PCR and was found to be slightly modulated with a biphasic response; first down- then up-regulated (figure 2).

The SLC12A5 gene, coding for an integral membrane K-Cl cotransporter, was repressed at 4 h then strongly increased at 24 h (qRT-PCR ratio : 4.3). The three genes, FBXL2, ZNT1 and SLC12A5, are therefore appropriate candidates for research into cobalt transport proteins in human cells and will be studied further using targeted biological approaches.

Potential biomarkers and secreted proteins

Biological monitoring of exposure to cobalt is mainly based on its concentration in urine. This method may be inadequate for two reasons. Firstly, measuring the quantity of excreted metal (exposure marker) does not necessarily reflect organ damage, which varies from one person to another. Furthermore, depending on the chemical form, excretion does not necessarily reflect the concentration of metal in the body; indeed, certain forms are eliminated in several phases; only 40% of cobalt oxides inhaled are excreted within 72 h [45]. In this context, only an effective biomarker could evaluate the damage and provide a valuable tool for monitoring occupational contamination in the event of accidental acute exposure or chronic exposure. We selected nine modulated genes coding for secreted proteins: ICAM1 (CD54), SERPIN C1 (antithrombin-III), IFNA4 (alpha-4 interferon), TFRC (transferrin receptor), PLA2R1 (phospholipase receptor A2), TNFSF6 (FAS ligand), WFDC2 (WAP four-disulfide core domain protein 2), LTBP3 (latent transforming growth factor beta binding protein 3) and TIMP2 (tissue inhibitor of metalloproteinase 2). To further evaluate the interest of these genes, we checked the secretion modulation of the corresponding proteins. A549 cells were exposed to cobalt in a dose-response manner, to 0.2, 1 and 2 mM, and the supernatants analysed using available immunoassays. The levels of FAS ligand, WFDC2 [46], ICAM1, antithrombin-III, alpha-4 interferon, and soluble transferrin receptor were found to be invariant or below the kit detection limit. These results indicate that most of these immunoassays were not sensitive enough to detect faint modulations. TIMP2 protein (tissue inhibitor of metalloproteinase 2) was measured using ELISA and we confirmed the modulation recorded at gene level (microarray ratio: 0.63, qRT-PCR ratio: 0.69). A strong decrease in TIMP2 supernatant concentration was noted at 1 or 2 mM of cobalt (figure 4). The changes in TIMP2 levels in cell supernatants between cobalt and the control were -5% for 0.2 mM, -57% for 1 mM and -72% for 2 mM. The reason for this reduction in TIMP2 can be explained as a consequence of cobalt exposure. Matrix metalloproteinases (MMPs) are a family of endopeptidases that can degrade all the components of the extracellular matrix [47]. Endogenous protease inhibitors, known as tissue inhibitors of metalloproteinase (TIMP), provide critical extracellular regulation of MMP proteolytic activities [48, 49], regulating tumor growth, progression, and angiogenesis in a variety of experimental cancer models and human malignancies. TIMP2 mRNA is described as being suppressed by hypoxia in endothelial cells [50]. Exposure of cells to cobalt up-regulates indirectly hypoxia-inducible genes by a mechanism which has recently been explained [51]. Therefore, TIMP-2 down-regulation might be a consequence of cobalt-simulated hypoxia. Down-regulation of TIMP2 mediated by ROS production has been observed by Rezzani in rat cardiac tissue exposed to cyclosporin [52]. Since it is known that cobalt triggers ROS production [16], TIMP2 modulation by cobalt could also be explained by this pathway.

Figure 4
figure 4

TIMP2 ELISA. Concentrations of TIMP2 in supernatants were determined by ELISA following treatment of A549 cells with cobalt over 24 h. The changes in level of TIMP2 in cell supernatants between cobalt and the control were -5% for 0.2 mM, -57% for 1 mM and -72% for 2 mM (n = 2).

TIMP2 protein modulation in cobalt-exposed cell supernatant corroborated the gene down-regulation observed on the microarrays and with qRT-PCR. TIMP2 modulation has never been associated with metal stress. This result is innovative and will be studied further in the biological fluids of rats exposed to different chemical forms of cobalt by inhalation.


This study provides the first toxicogenomic analysis of human lung cell response to acute cobalt exposure. We have confirmed that genes involved in the cobalt hypoxia response and apoptosis are modulated. We have also revealed genes linked to heat shock response and proteasome function that have already been described in other metal stress responses. Newly identified genes linked to cobalt acute toxicity include potential cobalt carriers (FBXL2, ZNT1, SLC12A5) and tumor suppressors or transcription factors (MAZ, DLG1, MYC, AXL). Some of these genes provide new hypotheses for elucidating the mechanisms of cobalt intracellular chemical toxicity. Targeted biological approaches might confirm their biochemical role in the cobalt response.

Regarding biomarkers, we have highlighted the down-regulation of TIMP2, a gene coding for a secreted protein. TIMP2 modulation was confirmed at protein level, in a dose-dependent manner, in the supernatant of exposed A549 cells. TIMP2 provides a putative biomarker of cobalt toxicity that will be further studied on animal models.


Cell culture

The human type II epithelial cell line A549 from ATCC, was cultured as described previously [53]. At midlog phase, a medium without FCS, containing cobalt (CoCl2, Sigma) or not, was added for 24 hours. A CoCl2 stock solution (1 M in MilliQ water) was prepared extemporaneously, filtered through a sterile 0.22 μm membrane and then diluted in the cell culture medium (final concentration 2 mM). After 30 mn, 2 h, 4 h or 24 h, the cells were harvested with trypsin and washed in PBS containing 1 mM EDTA. Culture supernatants for ELISA and cells for Western blot confirmations were collected after 24 h of exposure.

Cytotoxicity studies

Cobalt cytotoxicity was determined by measuring the intracellular ATP after 24 h using the Promega Celltiter-GlowTM Luminescence Cell Viability Assay (four to eight replicates per concentration). The viability rate was determined as the ratio between the ATP in treated cells and control cells. Cobalt concentration in cell pellets was measured using flame atomic absorption spectroscopy (FAAS, CERECO laboratory, Nimes, France).

RNA extraction

Total RNA was isolated using the Quiagen RNeasy miniprep kit according to the manufacturer's instructions. mRNAs were prepared using the Oligotex mRNA Quiagen kit, the purification process being repeated once to eliminate any contamination. RNA concentration was determined by OD measurement (260 nm/280 nm), purity and integrity were assessed using an Agilent 2100 Bioanalyser.

Microarrays and data analysis

Microarrays were obtained from the Service de Genomique Fonctionnelle (CEA, Evry, France). Two types of DNA collections were used for preparing the cDNA arrays: a collection of 5760 full length cDNA clones from the human infant brain 1NIB library, kindly provided by Genethon, and a collection of 2304 human PCR products (400–600 bp) corresponding to specific key-words selected by direct query in the Unigene database. A complete description of the microarrays used in this study, including the protocols for slide production, has been submitted to the GEO database [24] under accession number GPL4263.

For each microarray experiment, 20 μg of total RNA were reverse transcripted and indirectly labelled using the FairPlay Microarray Labelling Kit (Stratagene). Amino reactive Cy3- and Cy5- dyes (Cy™ Dye Post-Labelling Reactive Dye Pack, Amersham) were then chemically bound to cobalt and control cDNA respectively. Reverse labelling (or dyeswap) was performed for each experiment. Microarrays were pre-treated for 20 min with an N-methyl-pyrrolidinone solution containing 20 g/L succinic anhydride and 20 mM sodium borate, pH 8, immersed quickly in ethanol and dried by centrifugation for 6 min at 500 rpm. Following a quick rinse in RNase-free water, the slides were prehybridized in 30 ml 5X SSC, 1% SDS, 1% BSA (w/v) at 50°C for 40 min. The microarrays were then rinsed in RNAse-free water, immersed quickly in isopropanol and dried by centrifugation for 6 min at 500 rpm. The hybridization solution was prepared using 100 μl formamide, 10 μl SDS 10% and 30 μl RNase-free water. Labelled cDNA was solubilized in 17.5 μl of this solution plus 7.5 μl of 20 × SSPE and heat-denatured. Hybridization was performed overnight at 42°C. The microarrays were quickly rinsed in 0.1 × SSC, 0.01% SDS, washed twice for 10 min in 0.01 × SSC, then dried quickly in a stream of nitrogen and scanned with GenePix 4000B (Axon Instrument Inc., Forster City, CA). Cy3 and Cy5 spot fluorescence intensities were quantified after local subtraction of background using Genepix Pro 4.0 software (Axon Instrument Inc.). For each time point the result files were submitted to GeneSpring software 6.2 (Agilent Technologies) as follows. The data were first converted to take into account the results of the dyeswap reverse labelling, then normalized using the Lowess method, applying robust locally-weighted regression to smooth the intensity-dependence of the log ratios. The normalized data were then filtered on a quality test basis. This involved selecting spots detected on at least half of the microarrays with at least 70% pixels above threshold intensity (set to the median background plus two standard deviations).

From these remaining spots, we selected those with fluorescence ratios (representing cobalt-treated cells versus control) above 1.5 fold with p-value < 0.05 using a t-test statistical analysis on Genespring software and performing a Benjamini and Hochberg false discovery rate multiple testing correction.

Quantitative real-time polymerase chain reaction (qRT-PCR)

Specific primers were designed with Primer3 [54] using cDNA sequences spotted on microarrays, and amplicons were controlled using Mfold [55]. The list of primers is given in table 3. Before differential analysis, good primer efficiency (80% -120 %) was checked as a mandatory test. One μg of total RNA (for 30 min, 2 h and 4 h) or 50 ng of mRNA (for 24 h) were reverse transcribed with oligo dT (12–18) using Omniscript reverse transcriptase (Qiagen) following the manufacturer's instructions. cDNA was diluted at least 10-fold with DNase-free water to 10 ng/μL. PCR was performed on 15 μl using a DyNAmoTM HS SYBR Green qPCR Kit (Finnzymes) on a VWR DNA Engine Opticon® 2 system. The amplification program consisted of 1 cycle at 95°C with 10 min-hold followed by 40 cycles at 95°C with 15 sec-hold, 60°C with 1 min-hold, and a reading step at 60°C for 1 sec. Amplification was followed by melting curve analysis between 65°C and 95°C. RNA was used as a negative template for the absence of residual genomic DNA and a negative control without cDNA was used to control overall specificity. House- keeping genes, which are generally used as reference genes, are often modulated by stress. For example, GAPDH is strongly induced by cobalt (ratio 3). To find an invariant gene at each time point, several candidates were selected for their invariance from all the microarray data, and tested as invariant genes in qRT-PCR. FIGF (c- fos induced growth factor) was selected for time points 0.5 h, 2 h and 4 h. For samples at 24 h, mRNA had to be purified from total RNA to detect TUBA3 as an invariant reference gene. Differential analysis was performed on cDNA templates obtained from cobalt-treated or untreated cells in sextuplets for each gene, in the same way as for the reference gene (5 ng for cDNA from total RNA or 0.1 ng for cDNA from mRNA). The results were processed using REST-MCS software [56] and tested for significance using a Pair Wise Fixed Reallocation Randomized Test by calculating a p value. Modulation was considered to be significant when p calculated by REST software was < 0.05.

Western blot

Western blots were performed as previously described [53] after 24 h of cobalt exposure. Cells were lysed with 7 M urea, 2 M thiourea, 4% CHAPS, 20 mM spermine base, anti proteases (Roche cocktail), 40 mM DTT or 50 mM Tris-HCl, NaCl 150 mM ph7.5, 20 mM spermine base, anti proteases, 0.5% NP40 and 1% DOC. 10 to 50 μg of proteins were loaded onto a 4–12% or 12% NuPAGE gel in a MOPS or MES buffer (Invitrogen). Primary antibodies were: mouse monoclonal anti-DLG1 (sc-9961, 1/100, Santa Cruz), rabbit polyclonal anti-hsp90β (Ab1, 1/1666, NeoMarkers), mouse monoclonal anti-βactin (AC-15, 1/2000, Sigma), or mouse monoclonal anti-ubiquitin (P4D1, 1/400, Santa Cruz). A VECTASTAIN ABC kit (Vector laboratories) or a goat anti-mouse IgG coupled to HRP (1/10000, Novagen) was used for detection.


Immunoassays were performed on culture supernatants from cells exposed to various concentrations of cobalt over a 24 h period. The supernatants were tested crude or at 10 fold concentration. TIMP2 and FASL proteins were detected using ELISA kits (Raybiotech and R&D systems respectively) according to the manufacturer's instructions. The WFDC2 protein, was analysed in Dr Hellström's laboratory because they developed the WFDC2 ELISA [46]. Other immunoassays were performed by the CERBA laboratory (Cergy Pontoise, France).


  1. ATSDR: toxicological profile for cobalt. 2004, []

    Google Scholar 

  2. Barceloux DG: Cobalt. J Toxicol Clin Toxicol. 1999, 37 (2): 201-206. 10.1081/CLT-100102420.

    Article  CAS  PubMed  Google Scholar 

  3. Hartmann P: Inventaire des radioexpositions [in French]. EDF - Service de radioprotection. 2000, 17 (Radiotoxicologie): 2-5.

    Google Scholar 

  4. Le Guen B, Ansoborlo E: Le cobalt et ses isotopes [in French]. Toxicologie-Pathologie professionnelle. 2005, Elsevier, 1-11.

    Google Scholar 

  5. Chlorinated drinking-water; chlorination by-products; some other halogenated compounds; cobalt and cobalt compounds. International Agency for Research on Cancer (IARC) Working Group, Lyon, 12-19 June 1990. 1991, 52: 1-544.

  6. Goldoni M, Catalani S, De Palma G, Manini P, Acampa O, Corradi M, Bergonzi R, Apostoli P, Mutti A: Exhaled breath condensate as a suitable matrix to assess lung dose and effects in workers exposed to cobalt and tungsten. Environ Health Perspect. 2004, 112 (13): 1293-1298.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  7. Rosenberg N: Allergie respiratoire professionnelle au cobalt. INRS. 1989, 40TR15:

    Google Scholar 

  8. INERIS: Cobalt et ses dérivés [in French]. fiche de données toxicologiques et environnementales des substances chimiques. 2004, version1 (3- décembre 04): 1-52.

    Google Scholar 

  9. Palmes ED, Nelson N, Laskin S, Kuschner M: Inhalation toxicity of cobalt hydrocarbonyl. Am Ind Hyg Assoc J. 1959, 20: 453-468.

    Article  CAS  PubMed  Google Scholar 

  10. Bucher J: NTP technical report on the toxicity studies of Cobalt Sulfate Heptahydrate in F344/N Rats and B6C3F1 Mice (Inhalation Studies) (CAS No. 10026-24-1). Toxic Rep Ser. 1991, 5: 1-38. 10.1016/0890-6238(91)90105-O.

    Article  PubMed  Google Scholar 

  11. Bucher JR, Elwell MR, Thompson MB, Chou BJ, Renne R, Ragan HA: Inhalation toxicity studies of cobalt sulfate in F344/N rats and B6C3F1 mice. Fundam Appl Toxicol. 1990, 15 (2): 357-372. 10.1016/0272-0590(90)90061-N.

    Article  CAS  PubMed  Google Scholar 

  12. Alexander CS: Cobalt-beer cardiomyopathy. A clinical and pathologic study of twenty-eight cases. Am J Med. 1972, 53 (4): 395-417. 10.1016/0002-9343(72)90136-2.

    Article  CAS  PubMed  Google Scholar 

  13. Jacobziner H, Raybin HW: Poison control...accidental cobalt poisoning. Arch Pediatr. 1961, 78: 200-205.

    CAS  PubMed  Google Scholar 

  14. De Boeck M, Kirsch-Volders M, Lison D: Cobalt and antimony: genotoxicity and carcinogenicity. Mutat Res. 2003, 533 (1-2): 135-152.

    Article  CAS  PubMed  Google Scholar 

  15. Lison D, De Boeck M, Verougstraete V, Kirsch-Volders M: Update on the genotoxicity and carcinogenicity of cobalt compounds. Occup Environ Med. 2001, 58 (10): 619-625. 10.1136/oem.58.10.619.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  16. Salnikow K, Su W, Blagosklonny MV, Costa M: Carcinogenic metals induce hypoxia-inducible factor-stimulated transcription by reactive oxygen species-independent mechanism. Cancer Res. 2000, 60 (13): 3375-3378.

    CAS  PubMed  Google Scholar 

  17. Pulido MD, Parrish AR: Metal-induced apoptosis: mechanisms. Mutat Res. 2003, 533 (1-2): 227-241.

    Article  CAS  PubMed  Google Scholar 

  18. Bruick RK: Oxygen sensing in the hypoxic response pathway: regulation of the hypoxia-inducible transcription factor. Genes Dev. 2003, 17 (21): 2614-2623. 10.1101/gad.1145503.

    Article  CAS  PubMed  Google Scholar 

  19. Harris AL: Hypoxia--a key regulatory factor in tumour growth. Nat Rev Cancer. 2002, 2 (1): 38-47. 10.1038/nrc704.

    Article  CAS  PubMed  Google Scholar 

  20. Li Q, Chen H, Huang X, Costa M: Effects of 12 metal ions on iron regulatory protein 1 (IRP-1) and hypoxia-inducible factor-1 alpha (HIF-1alpha) and HIF-regulated genes. Toxicol Appl Pharmacol. 2006, 213 (3): 245-255. 10.1016/j.taap.2005.11.006.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Lettieri T: Recent applications of DNA microarray technology to toxicology and ecotoxicology. Environ Health Perspect. 2006, 114 (1): 4-9.

    CAS  PubMed Central  PubMed  Google Scholar 

  22. Lieber M, Smith B, Szakal A, Nelson-Rees W, Todaro G: A continuous tumor-cell line from a human lung carcinoma with properties of type II alveolar epithelial cells. Int J Cancer. 1976, 17 (1): 62-70. 10.1002/ijc.2910170110.

    Article  CAS  PubMed  Google Scholar 

  23. Crouch SP, Kozlowski R, Slater KJ, Fletcher J: The use of ATP bioluminescence as a measure of cell proliferation and cytotoxicity. J Immunol Methods. 1993, 160 (1): 81-88. 10.1016/0022-1759(93)90011-U.

    Article  CAS  PubMed  Google Scholar 

  24. Andrew AS, Warren AJ, Barchowsky A, Temple KA, Klei L, Soucy NV, O'Hara KA, Hamilton JW: Genomic and proteomic profiling of responses to toxic metals in human lung cells. Environ Health Perspect. 2003, 111 (6): 825-835.

    CAS  PubMed Central  PubMed  Google Scholar 

  25. National Center for Biotechnology Information. []

  26. DAVID: DAVID. []

  27. Maxwell P, Salnikow K: HIF-1: an oxygen and metal responsive transcription factor. Cancer Biol Ther. 2004, 3 (1): 29-35.

    Article  CAS  PubMed  Google Scholar 

  28. Semenza GL: Targeting HIF-1 for cancer therapy. Nat Rev Cancer. 2003, 3 (10): 721-732. 10.1038/nrc1187.

    Article  CAS  PubMed  Google Scholar 

  29. Vengellur A, Woods BG, Ryan HE, Johnson RS, LaPres JJ: Gene expression profiling of the hypoxia signaling pathway in hypoxia-inducible factor 1alpha null mouse embryonic fibroblasts. Gene Expr. 2003, 11 (3-4): 181-197.

    Article  CAS  PubMed  Google Scholar 

  30. Vengellur A, Phillips JM, Hogenesch JB, LaPres JJ: Gene expression profiling of hypoxia signaling in human hepatocellular carcinoma cells. Physiol Genomics. 2005, 22 (3): 308-318. 10.1152/physiolgenomics.00045.2004.

    Article  CAS  PubMed  Google Scholar 

  31. Prat O, Berenguer F, Malard V, Tavan E, Sage N, Steinmetz G, Quemeneur E: Transcriptomic and proteomic responses of human renal HEK293 cells to uranium toxicity. Proteomics. 2005, 5 (1): 297-306. 10.1002/pmic.200400896.

    Article  CAS  PubMed  Google Scholar 

  32. Gerner C, Frohwein U, Gotzmann J, Bayer E, Gelbmann D, Bursch W, Schulte-Hermann R: The Fas-induced apoptosis analyzed by high throughput proteome analysis. J Biol Chem. 2000, 275 (50): 39018-39026. 10.1074/jbc.M006495200.

    Article  CAS  PubMed  Google Scholar 

  33. Bachmair A, Varshavsky A: The degradation signal in a short-lived protein. Cell. 1989, 56 (6): 1019-1032. 10.1016/0092-8674(89)90635-1.

    Article  CAS  PubMed  Google Scholar 

  34. Hershko A, Leshinsky E, Ganoth D, Heller H: ATP-dependent degradation of ubiquitin-protein conjugates. Proc Natl Acad Sci U S A. 1984, 81 (6): 1619-1623. 10.1073/pnas.81.6.1619.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Stubbs MC, Min I, Izzo MW, Rallapalli R, Derfoul A, Hall DJ: The ZF87/MAZ transcription factor functions as a growth suppressor in fibroblasts. Biochem Cell Biol. 2000, 78 (4): 477-485. 10.1139/bcb-78-4-477.

    Article  CAS  PubMed  Google Scholar 

  36. Kopera E, Schwerdtle T, Hartwig A, Bal W: Co(II) and Cd(II) substitute for Zn(II) in the zinc finger derived from the DNA repair protein XPA, demonstrating a variety of potential mechanisms of toxicity. Chem Res Toxicol. 2004, 17 (11): 1452-1458. 10.1021/tx049842s.

    Article  CAS  PubMed  Google Scholar 

  37. Witkiewicz-Kucharczyk A, Bal W: Damage of zinc fingers in DNA repair proteins, a novel molecular mechanism in carcinogenesis. Toxicol Lett. 2006, 162 (1): 29-42. 10.1016/j.toxlet.2005.10.018.

    Article  CAS  PubMed  Google Scholar 

  38. Humbert P, Russell S, Richardson H: Dlg, Scribble and Lgl in cell polarity, cell proliferation and cancer. Bioessays. 2003, 25 (6): 542-553. 10.1002/bies.10286.

    Article  CAS  PubMed  Google Scholar 

  39. Liu J, Xie Y, Ducharme DM, Shen J, Diwan BA, Merrick BA, Grissom SF, Tucker CJ, Paules RS, Tennant R, Waalkes MP: Global gene expression associated with hepatocarcinogenesis in adult male mice induced by in utero arsenic exposure. Environ Health Perspect. 2006, 114 (3): 404-411.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Nakano T, Tani M, Ishibashi Y, Kimura K, Park YB, Imaizumi N, Tsuda H, Aoyagi K, Sasaki H, Ohwada S, Yokota J: Biological properties and gene expression associated with metastatic potential of human osteosarcoma. Clin Exp Metastasis. 2003, 20 (7): 665-674. 10.1023/A:1027355610603.

    Article  CAS  PubMed  Google Scholar 

  41. Janssen JW, Schulz AS, Steenvoorden AC, Schmidberger M, Strehl S, Ambros PF, Bartram CR: A novel putative tyrosine kinase receptor with oncogenic potential. Oncogene. 1991, 6 (11): 2113-2120.

    CAS  PubMed  Google Scholar 

  42. Kobayashi M, Shimizu S: Cobalt proteins. Eur J Biochem. 1999, 261 (1): 1-9. 10.1046/j.1432-1327.1999.00186.x.

    Article  CAS  PubMed  Google Scholar 

  43. Conklin DS, Kung C, Culbertson MR: The COT2 gene is required for glucose-dependent divalent cation transport in Saccharomyces cerevisiae. Mol Cell Biol. 1993, 13 (4): 2041-2049.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. Conklin DS, McMaster JA, Culbertson MR, Kung C: COT1, a gene involved in cobalt accumulation in Saccharomyces cerevisiae. Mol Cell Biol. 1992, 12 (9): 3678-3688.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Kreyling WG, Cox C, Ferron GA, Oberdorster G: Lung clearance in Long-Evans rats after inhalation of porous, monodisperse cobalt oxide particles. Exp Lung Res. 1993, 19 (4): 445-467.

    Article  CAS  PubMed  Google Scholar 

  46. Hellstrom I, Raycraft J, Hayden-Ledbetter M, Ledbetter JA, Schummer M, McIntosh M, Drescher C, Urban N, Hellstrom KE: The HE4 (WFDC2) protein is a biomarker for ovarian carcinoma. Cancer Res. 2003, 63 (13): 3695-3700.

    PubMed  Google Scholar 

  47. Coussens LM, Fingleton B, Matrisian LM: Matrix metalloproteinase inhibitors and cancer: trials and tribulations. Science. 2002, 295 (5564): 2387-2392. 10.1126/science.1067100.

    Article  CAS  PubMed  Google Scholar 

  48. Jiang Y, Goldberg ID, Shi YE: Complex roles of tissue inhibitors of metalloproteinases in cancer. Oncogene. 2002, 21 (14): 2245-2252. 10.1038/sj.onc.1205291.

    Article  CAS  PubMed  Google Scholar 

  49. Baker AH, Edwards DR, Murphy G: Metalloproteinase inhibitors: biological actions and therapeutic opportunities. J Cell Sci. 2002, 115 (Pt 19): 3719-3727. 10.1242/jcs.00063.

    Article  CAS  PubMed  Google Scholar 

  50. Ben-Yosef Y, Lahat N, Shapiro S, Bitterman H, Miller A: Regulation of endothelial matrix metalloproteinase-2 by hypoxia/reoxygenation. Circ Res. 2002, 90 (7): 784-791. 10.1161/01.RES.0000015588.70132.DC.

    Article  CAS  PubMed  Google Scholar 

  51. Salnikow K, Donald SP, Bruick RK, Zhitkovich A, Phang JM, Kasprzak KS: Depletion of intracellular ascorbate by the carcinogenic metals nickel and cobalt results in the induction of hypoxic stress. J Biol Chem. 2004, 279 (39): 40337-40344. 10.1074/jbc.M403057200.

    Article  CAS  PubMed  Google Scholar 

  52. Rezzani R, Giugno L, Buffoli B, Bonomini F, Bianchi R: The protective effect of caffeic acid phenethyl ester against cyclosporine A-induced cardiotoxicity in rats. Toxicology. 2005, 212 (2-3): 155-164. 10.1016/j.tox.2005.04.020.

    Article  CAS  PubMed  Google Scholar 

  53. Malard V, Prat O, Darrouzet E, Berenguer F, Sage N, Quemeneur E: Proteomic analysis of the response of human lung cells to uranium. Proteomics. 2005, 5 (17): 4568-4580. 10.1002/pmic.200402038.

    Article  CAS  PubMed  Google Scholar 

  54. Primer3: Primer3. Primer3. []

  55. Mfold: Mfold. []

  56. Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR - Nucleic Acids Res. 2002, 30 (9): e36-10.1093/nar/30.9.e36.

    Article  PubMed Central  PubMed  Google Scholar 

Download references


This research was supported in part by "Electricité De France". We thank the Service de Genomique Fonctionnelle (CEA Evry, France) for providing us with microarrays and Dr I.Hellström (Washington University) for performing the WFDC2 ELISA test. We thank Philippe Guérin and Nicole Sage for their technical support. We thank Elisabeth Darrouzet for critical review of this manuscript and for helpful discussions.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Veronique Malard.

Additional information

Authors' contributions

VM designed the study and drafted the manuscript.

FB carried out data acquisition and analysis and prepared the microarray results for tables and figures.

OP was involved in the design stage, qRT-PCR analysis and microarray experiments.

SR carried out qRT-PCR.

GS carried out the RNA extraction, labelling and microarray hybridization.

EQ was involved in the design stage and revised the manuscript

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 (, 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

Malard, V., Berenguer, F., Prat, O. et al. Global gene expression profiling in human lung cells exposed to cobalt. BMC Genomics 8, 147 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: