Carbon nanotubes (CNT) are among the most widely used nanomaterials for various industrial and biomedical applications [1-3]. Their unique properties of high tensile strength, flexibility, thermal and electrical conductivity, and their light weight have contributed to their widespread use. However, their rapid growth in utility has raised a major concern about their safety, especially on human health and environment. CNT share several properties with asbestos fibers such as high aspect ratio, durability, biopersistence, and mode of exposure (e.g., inhalation). Because asbestos fibers are classified as a Group I human carcinogen by the International Agency of Research in Cancer , there is a concern about the potential carcinogenicity of CNT . Exposure to asbestos has been shown to cause various lung diseases including lung cancer and mesothelioma . Thus, there is a clear need for carcinogenicity studies of CNT to develop control and prevention strategies .
Recent studies have shown that both single-walled (SW) and multi-walled (MW) CNT deposit in the deep lung tissue of mice with low clearance after pulmonary administration [7-9]. Occupationally relevant in vivo exposures (10 to 80 μg/mouse) and pulmonary fibroblast in vitro exposure studies, using in vivo dose equivalents (0.02 to 0.2 μg/cm2), resulted in dose-dependent transient pulmonary inflammation followed by fibroblast cell proliferation, alveolar wall thickening, and collagen I production culminating in persistent pulmonary fibrosis [9,10]. Several in vitro and in vivo studies have suggested the potential carcinogenicity of CNT [11-13]. The possible carcinogenic mechanisms include DNA damage [14-16], mitotic disruption , impaired apoptosis, and activation of oncogenic signaling events, which recapitulate asbestos-induced lung cancer and mesothelioma [17-19]. It was suggested that nanotube bundles are similar in size to microtubules that form the mitotic spindle and may be incorporated into the mitotic spindle apparatus, resulting in multipolar mitosis and aneuploidy  which were observed in asbestos-treated cells . Furthermore, CNT exposure induced asbestos-like granulomas in mice in a short-term abdominal instillation study .
Multiple genes may be involved in single-walled carbon nanotubes (SWCNT) pathogenesis. For examples, p38 MAPK was reported to regulate SWCNT-induced fibrogenic and angiogenic responses . Other short-term exposure studies report increased NF-kβ, TNFα, NRLP3, iNOS, p53, and TGFβ activity all culminating in lung inflammation and the onset of pulmonary interstitial fibrosis [22,23]. Long-term in vitro and in vivo exposure studies found that SWCNT and MWCNT exposure to lung tissues resulted in KRAS activation, enhanced micronuclei and nuclear protrusions, mitosis disruption, mutated p16, cMyc proto-oncogene signaling, and enhanced metalloproteinase remodeling of the extracellular matrix [9,24-27], which all indicate tumorigenesis potential following CNT exposure. Recent transcriptomic and sequencing studies suggest that long-term CNT exposure results in similar signaling compared to other biopersistent fibers (i.e., asbestos) but also display unique signaling pathways [27-30] raising the possibility of a distinct CNT pathogenesis paradigm . Although a majority of these pro-inflammation and pro-fibrotic genes have known roles in cancer, a consensus on those key molecules and signaling pathways that potentially drive and contribute to CNT-associated carcinogenesis remain elusive.
Our recent studies showed that low-dose (0.02 μg/cm2) SWCNT exposure caused malignant transformation and tumorigenesis of human lung bronchial epithelial cells . We also found that phosphorylation at various positions on p53 in these malignant cells significantly decreased after sub-chronic SWCNT exposure . These phosphorylation sites are crucial for stabilization and activation of p53-dependent functions of tumor suppression, DNA damage repair, and apoptosis resistance [32-35]. At present, the upstream and downstream signal transduction pathways associated with altered p53 status contributing to a SWCNT-induced epithelial tumorigenic phenotype are unknown.
Since carcinogenesis is a long-term, multi-step process and involves multiple genes, it is of critical importance to study the chronic exposure and genome-wide expression changes to define the underlying molecular mechanisms of carcinogenesis, which is largely unknown. To accomplish this goal, we have developed a sub-chronic exposure model in which human lung bronchial epithelial cells, the primary target of CNT inhalation exposure, were continuously exposed to low-dose SWCNT in culture over a long period (6 months) . Because isolated primary lung cells do not survive in long-term culture, we employed an immortalized human lung epithelial BEAS-2B cell line which is non-tumorigenic and exhibits similar cellular responses [7,20,36]. The long-term SWCNT-exposed BEAS-2B (B-SWCNT) cells were evaluated, for the first time, for genome-wide expression profiling and functional analysis compared to unexposed, passage-matched control cells. Our study revealed multiple genes and signaling networks that are affected by SWCNT-induced oncogenic transformation which might serve as potential SWCNT-specific exposure or disease markers. This in vitro approach also supports prudent adoption of exposure control strategies protection of workers, consumers, and the environment.
Preparation of single-walled carbon nanotubes
SWCNT (CNI, Houston, TX, USA) were produced by the high-pressure CO disproportionation (HiPco) technique, employing CO in a continuous-flow gas phase as the carbon feedstock and Fe(CO)5 as the iron-containing catalyst precursor and used in our previous study . Briefly, the SWCNT were purified by acid treatment to remove metal contaminates. Elemental analysis of the supplied SWCNT showed that the SWCNT were 99% elemental carbon and 0.23% iron. The specific surface area measured at −196°C by the nitrogen absorption-desorption technique (Brunauer-Emmet-Teller method) was 400 to 1,000 m2/g. The diameter and length distribution of the SWCNT measured by field emission scanning electron microscopy which were 0.8 to 1.2 nm and 0.1 to 1 μm, respectively. SWCNT were dispersed by acetone/sonication method as previously described . Briefly, SWCNT were treated with acetone and placed in an ultrasonic bath for 24 h. The dispersed CNT were then filtered from the solution using a 20-μm nylon mesh screen followed by a 0.2-μm polytetrafluoroethylene filter. After filter collection, the dispersed CNT were washed thoroughly with distilled water and suspended in phosphate-buffered saline (PBS) with 2- to 3-min sonication (Sonic Vibra Cell Sonicator, Sonic & Material Inc, Newtown, CT, USA).
In vitro treatment and cell sample preparation for microarrays
Microarray sample collection and analyses were performed on unexposed and SWCNT-exposed cells from our previous study . Briefly, human lung bronchial epithelial BEAS-2B cells were continuously exposed for 6 months to a sub-cytotoxic concentration (0.02 μg/cm2; equivalent to 0.1 μg/ml) of dispersed SWCNT or PBS (non-treatment control) in a six-well culture plate in triplicate. Each replicate was exposed and independently assayed throughout the study. This relatively low concentration was chosen due to its relevance to in vivo SWCNT exposure dose of 10 μg/mouse previously reported [8,37,38]. Repeated, long-term dosing of particulate in in vitro systems may result in accumulation of particle burden over time in the assay system. To reduce and safeguard against potential SWCNT accumulation over time, every 3 to 4 days, the exposure media was aspirated and the cells were triplicate washed with PBS to remove suspended and unbound SWCNT. Cells were then resupplied with new exposure media containing SWCNT. The cells were passaged weekly at pre-confluent densities using a solution containing 0.05% w/v trypsin and 0.5 mM EDTA (Invitrogen, Carlsbad, CA, USA). Greater than 50% of cells were removed during passaging which resulted in a relatively small proportion of intracellular SWCNT remaining in the system. Parallel cell cultures grown in the same media without SWCNT provided passage-matched control. Following the 6-month exposure, SWCNT-treated (B-SWCNT) and control (B-Control) cells were evaluated for genome-wide gene expression. Three B-SWCNT cell biological replicates showing similar phenotypic behavior and p53 expression patterns  were pooled and analyzed as one. Additional filters such as increasing the fold change to >2 will limit this disadvantage. Since the objective of this study was to explore the candidate genes involved in the carcinogenic process for further studies, this study is considered a preliminary step.
RNA sample preparation and hybridization were performed as previously described . Briefly, total RNA was harvested and isolated using TRIzol (Invitrogen) according to manufacturer's instructions. Isolated RNA in DEPC-treated water was then immediately frozen at −80°C for 24 h and shipped on dry ice to ArrayStar (Rockville, MD, USA) for mRNA processing, microarray hybridization, and probe expression normalization.
DNA microarray and data analysis
Whole genome oligo microarrays (#014850, Agilent Technology, Santa Clara, CA, USA) were used. Each array represents more than 41,000 unique human genes and transcripts sourced from RefSeq, Goldenpath Ensembl Unigene Human Genome (Build 33) and GenBank databases. Fluorescent RNA targets were prepared using Agilent Quick Amp Labeling Kit. Microarray hybridization was performed at 65°C for 17 h in Agilent's SureHyb hybridization chambers. After washing in an ozone-free environment, the slides were scanned using the Agilent DNA microarray scanner (model G2505B).
Raw data were extracted using Agilent Feature Extraction Software (Santa Clara, CA, USA). Expression raw data from passage control cells were used for other exposure studies. The resulting text files were imported into the Partek Genomic Suite (PGS; St. Louis, MO, USA) for preprocessing, normalization, and statistical analysis. Briefly, normalization was performed using the Agilent FE one-color scenario (mainly median normalization). Genes marked present or marginal in all samples (‘All Targets Value’) were chosen for further data analysis. Differentially expressed genes (DEGs) were identified through fold-change screening and t-test assuming unequal variances. Genes were differentially expressed in B-SWCNT cells if they exhibited ≥ ±twofold expression and p ≤ 0.05 compared to B-Control cells. Unsupervised hierarchical clustering, gene ontology (GO) analysis, and pathway analysis were performed using the PGS (St. Louis, MO). False discovery rate (FDR) assessment by Benjimini Hochberg methods was applied for multiple hypothesis testing purpose . The filtered gene lists were used for GO, pathways, and gene-gene interaction analyses by using Ingenuity Pathway Analysis (IPA) software package (Redwood City, CA, USA).
Apoptosis protein array
To substantiate the genomic data and further identify proteins contributing to the apoptotic resistance phenotype, B-SWCNT and B-Control cell protein lysates were subjected to a human apoptosis protein expression array (see Additional file 1: Table S1 for array layout) in duplicate following manufacturer's instructions (R & D Systems, Minneapolis, MN, USA). Briefly, 2 × 106 cells from each treatment were seeded into a 25-cm2 flask in particle-free growth medium and cultured for 24 h. Next, plated cells were placed on ice, washed in ice cold PBS, and lysed using the manufacturer's lysis buffer. Cell lysate was collected via rubber policemen, rocked on ice in suspension of 30 min, and spun at 12.5 K × g to collect supernatant. Protein levels were determined via BCA absorbance method on a spectrophotometer at 562 nm. Five hundred micrograms of protein was incubated with each protein array overnight at 4°C. Each array was triplicate washed in buffer, incubated with detection antibody for 1 h, washed again, and incubated with streptavidin-HRP antibody for 30 min. Protein expression was captured using HRP chemiluminescence on X-ray film.
Bcl-2, pAkt, and Akt protein expression was determined using previously described Western blot methods . Briefly, B-SWCNT and B-Control cells were seeded (5 × 105 cells/well) and cultured in six-well plates for 24 h in particle-free medium or re-exposed to dispersed SWCNT (1 to 50 μg/ml or 0.1 to 5.2 μg/cm2). Cells were placed on ice, washed with ice cold PBS buffer, and lysed in protein buffer containing PMSF and protease inhibitor (Invitrogen). Cells were scraped, collected, and rocked on ice for 30 min. Following centrifugation at 12.5 K × g, BCA assay was used to measure protein concentration in collected supernatants. Equal amounts of protein were loaded into SDS PAGE gel, allowed to separate under current, and transferred to PVDF membranes using semi-dry transfer apparatus (Fisher Scientific, Pittsburgh, PA, USA). Membranes were incubated with antibodies to either Bcl-2, pAkt, Akt (Cell Signaling, Beverly, MA, USA), or B-actin (Sigma, St. Louis, MO, USA) overnight at 4°C, washed thrice in Tris buffer with 0.1% Tween, and incubated with anti-rabbit or anti-mouse HRP antibody for 1 h (Santa Cruz, Inc, Paso Robles, CA, USA). Protein levels were determined using Millipore Immobilon HRP chemiluminescent substrate (Millipore, Billerica, MA, USA) and exposure to X-ray film. Western blot experiments were independently performed three times. Quantification of Western blot results were performed by using ImageJ [41,42]. Statistical analyses were performed by using two group Student t-tests after quantification and normalization to corresponding controls.
Results and discussion
As we reported previously , transformed B-SWCNT cells displayed aberrant p53 signaling, increased aggressiveness, apoptosis resistance, anchorage-independent growth, angiogenicity, and in vivo xenograft tumor formation compared to unexposed passage control cells. These findings suggested the potential carcinogenicity of B-SWCNT cells. The present transcriptomic analysis and supporting data also indicate aberrant p53 signaling, pAkt, Ras, Notch 1, and altered intrinsic mitochondrial pathway signaling as potential mechanisms for SWCNT carcinogenicity. This is consistent with the previous observations that most human cancers possess p53 tumor suppressor gene inactivation resulting in dysfunctional signaling contributing to tumorigenesis [43,44].
Gene expression profile by hierarchical clustering
Hierarchical clustering is a convenient way to visualize the overall similarity among a large pool of samples. The heatmap pattern also offers clues of possible co-expression of genes. As shown in Figure 1, three control samples showed similar gene expression patterns to cluster on the left (control, control-1, and control-2). The B-SWCNT cells possessed a different expression pattern on the right (SWCNT1_3). Here, three treated samples were pooled and analyzed as one. This is a cost-effective way to get the most likely candidate genes but with the disadvantage of lowering the power and introducing more false positives. Since the objective of this study was to explore the candidate genes involved in the carcinogenic process for further studies, this study is considered a preliminary step.
Gene ontology (GO) analysis for biological processes and molecular functions
We analyzed top GO biological processes that are affected by SWCNT exposure based on the list of DEGs in B-SWCNT cells. Genes that are involved in immune response were among the top biological processes (Table 1). The top biological processes include the responses to metal ion, antigen process and presenting, regulation of phagocytosis, and vacuolar transport. Apparently, lung epithelial cells may recognize SWCNT as foreign which can alter the innate immune response and influence the immune system as reported [27,45,46]. In addition, the positive regulation of microtubule depolymerization suggests that long-term SWCNT exposure altered microtubule dynamics as previously reported [12,47]. Similarly, genome wide changes of GO molecular functions after SWCNT exposure included genes responsible for immune response (i.e., MHC class I) protein binding, post-transcriptional modifications (e.g., ubiquitin protease activity), and other enzymatic activities were among the top molecular function changes (Table 2). Previous ‘omics’ studies using long-term SWCNT exposure in vivo or in vitro models reported similar findings that protein post-translational modification constitutes a large GO signal [27,28] and is likely to play a large role in disease development and pathology.
Functional analysis of gene-gene interactions
Multiple gene-gene interactions are important for gene functions and signal transduction. Therefore, gene-gene interaction networks were plotted, scored, and ranked to identify interactions potentially driving the observed B-SWCNT malignant phenotype. Consistent with the observed carcinogenic and apoptosis-resistant phenotype of B-SWCNT cells , immune response, cell growth, cell death/survival, and cell cycle control signaling network were among the major functions identified in the top-ranked networks in the B-SWCNT toxicogenomic profile (Table 3). The functional involvement of individual genes or pathways in the carcinogenic process needs to be further studied individually to determine their actual involvement. Regardless, the pathways of Cellular growth and proliferation, hematological system development and function, humoral immune response, cell death and survival, embryonic development, organismal development, and cancer-related gene-gene interaction networks possessed gene-gene interactions known to participate in tumorigenic and CNT-induced neoplastic signaling mechanisms .
Aberrant p53 signaling
Our initial study reported B-SWCNT malignant transformed cells possessing an apoptotic resistant phenotype due in part to aberrant p53 signaling . Apoptosis plays an essential role in the removal of mutated or transformed cells and its disruption contributes to abnormal cell growth and malignancy [48,49]. To expand on this finding, we first examined p53 (TP53) and p53-associated protein mRNA expression in B-SWCNT compared to B-Control cells. Moderate decreased expression of TP53 mRNA was observed, but no statistical significance was found (Table 4). Although 13 other p53-associated transcripts showed no statistical difference, three transcripts experienced increased or decreased expression (>twofold) which included TP53I3, TP53INP1, and TP53INP2. This data suggested a moderate alteration of TP53-related signaling potentially due to post-transcriptional modifications. Overexpressed TP53I3, an oxidoreductase, is typically observed as a response to reactive oxidative species (ROS) stress which has been shown to occur in B-SWCNT cells [10,50]. More interestingly, decreased TP53INP1 and TP53INP2 expression suggests a decreased ability in B-SWCNT cells to respond to cell stress. Both genes encode for antioxidant proteins that respond to cell stress by positively regulating autophagy and directly phosphorylating serine 46 on p53 protein which is known to increase p53's apoptotic signaling ability . BEAS-2B cells exhibit increased autophagy in response to ROS stress, including CNT and other inhalable xenobiotic exposure, and dysregulation of autophagy may contribute to tumor promotion potential [50,51]. Decreased expression of these two transcripts suggests that B-SWCNT cells lost two major cell protection mechanisms including a major p53 pro-apoptotic signaling pathway and positive regulation of autophagy.
To identify pathways that potentially altered p53 signaling and fate, several canonical pathways with DEGs with known p53 functions were plotted in IPA. Upregulation of PTEN and lowered expression of AKT2 suggested decreased Akt signaling upstream of MDM2 regulation of p53 (Figure 2). Also, decreased GSK3β expression suggests release of β-catenin from APC complex resulting in activation of Wnt signaling  and a potential block to MDM2-mediated p53 degradation via the proteosome.
To substantiate post-translational phosphorylation of p53 and identify other proteins contributing to B-SWCNT cell apoptotic resistance ability, we conducted an apoptosis protein expression array. Figure 3 shows, in addition to substantial decreased phosphorylation of Ser46, a target of TP53INP1/2, Ser15 and Ser392, targets of MDM2 and casein kinase II [13,33,53], which indicated an unstable and dysfunctional p53 response in contributing to B-SWCNT malignant phenotype and enhanced cancer signaling (Table 3, network 7). Other apoptotic-associated proteins were found to alter expression. A substantial decrease in claspin expression indicated loss of checkpoint-mediated ATR control of the cell cycle following DNA damage. Overexpressed heat shock protein 60 (HSPD1) in B-SWCNT cells may contribute to cell proliferation and apoptotic resistance in response to genotoxic stress (Table 3, networks 1 and 9) by stabilizing, surviving, and inhibiting p53 . Moderate decreased TNFR1, the major receptor for TNFα, indicated a reduced ability to respond to pro-death TNFα signal which was previously confirmed . Conversely, a fivefold overexpressed Fas receptor in the extrinsic apoptotic pathway suggested a sensitization to apoptosis; however, it is possible that downstream modulators of Caspase 8 activity (i.e., cFLIP) may block this effect. Bcl-xl and Xiap, anti-apoptotic Bcl-2 family proteins, experienced a moderate decrease (twofold) in B-SWCNT compared to B-Control cells. Cells undergoing tumorigenic transformation typically exhibit changes in both pro-survival and death proteins . Investigation of both extrinsic and intrinsic signaling pathways and their role in B-SWCNT apoptotic resistance is currently under investigation in our group.
A recent demonstration of MWCNT-induced mesothelioma in p53 knockout mice [19,55] further supports our finding and the role of p53 in CNT carcinogenesis. A portion of this response may be due to a weakened p53 in the BEAS-2B cell model since the SV40 large T antigen immortalization procedure does affect p53 function . Numerous studies acknowledge that p53 status is a contributing factor to adverse effects following many different types of inhalable and respirable particle exposure to the lung [34,36,38,39,43]. In addition, other studies note that a CNT's physicochemical properties and toxicokinetics in combination with p53 status of exposed tissue determine disease response [56,57]. Regardless, SWCNT exposure did alter p53 signaling in these cells compared to B-Control cells. These findings add to the growing literature that pre-existing p53 condition along with physicochemical properties of CNTs contribute to disease outcome following respirable fiber exposure .
Since the above evidence suggested MDM2 activity and AKT1 is a major signaling protein in most cells compared to AKT2, we investigated whether B-SWCNT cells compared to B-Control cells exhibited increased AKT activation via phosphorylation, both in the presence and absence of SWCNT. Western blot analysis indicated increased ratio (1.7-fold, p < 0.05) of phosphorylated Akt (pAkt) to Akt protein expression in 0 to 5 μg/ml SWCNT in B-SWCNT compared to B-Control cells (Figure 4A,B). Exposure to SWCNT caused a moderate 20% increase in pAkt levels in B-Control cells at ≥10 μg/ml after 12 h exposure while re-exposure to B-SWCNT exhibited a 27% decrease at high dose levels (≥25 μg/ml). These results show that unexposed lung epithelial cells experience enhanced Akt activation following SWCNT exposure . In addition, it suggested that prolonged, low-dose SWCNT exposure results in a dysregulated, constitutively active pAKT and that Akt signaling may not be regulated at the mRNA level. Akt activation, in turn, potentially activated MDM2, leading to enhanced p53 ubiquitination and potential blockage of GSK3β (Figure 2). Thus, PI3K-AKT signaling may be involved in SWCNT oncogenesis by promoting cell survival (Table 3, network 3). Since pAkt is known to suppress intrinsic mitochondrial pathway and p53 blocks Bcl-2 expression, protein analysis showed that Bcl-2 protein expression was significantly elevated 2.5-fold in B-SWCNT compared to B-Control cells at 0 μg/ml (p < 0.05; Figure 4A,C). At 12 h post-exposure to SWCNT, neither cell type displayed a significant change in Bcl-2 expression. Following 24-h SWCNT exposure, B-Control cells displayed a significant increase in Bcl-2 expression at ≥5 μg/ml. Conversely, B-SWCNT displayed a significant decrease in Bcl-2 expression at ≥10 μg/ml. These dose-dependent Bcl-2 expression changes correlated with pAkt levels in each cell type. This suggested that SWCNT exposure to lung epithelial cells first results in Akt activation at 12 h followed by increased Bcl-2 at 24 h, which potentially occurs via either diminished p53 signaling to Noxa  or blockage of BAD activity . Persistent stimulation of this pathway could result in an apoptosis resistance phenotype via Bcl-2 overexpression. Interestingly, the dose-dependent drop in pAkt expression mirrored the drop in Bcl-2 expression (−28%) in re-exposed B-SWCNT cells, which resulted in equivalent pAKt and Bcl-2 expression levels in both cell types at high SWCNT doses (≥10 μg/ml). This suggests that pAkt expression sufficiently regulates Bcl-2-associated anti-apoptotic signaling in B-SWCNT cells  at low doses that typically occur in occupational exposures. Under unrealistic high doses, however, B-SWCNT cells may partially succumb to pro-apoptotic intrinsic signaling and may rely on other apoptotic resistance mechanisms (e.g., dysfunctional extrinsic pathway). In summary, prolonged SWCNT exposure coupled with p53 instability resulted in overexpressed pAkt and Bcl-2 which potentially plays a role in the apoptotic resistance and oncogenic development of B-SWCNT cells.
Genes involved in the molecular mechanism of cancer
Among the various possible molecular mechanisms of cancer, as shown in Additional file 2: Figure S1, both increased (RASGRP, ARHGEF9, SHC, and RHOB) and decreased gene expressions (RBPJK, JNK, DVL1, GSK3A, GSK3β, PAK, AKT2, PMAIP1 [Noxa], BAX) were observed in B-SWCNT cells. Based on overexpressed pAkt, decreased stabilization of p53 via several serine dephosphorylations and increased Bcl-2 protein, it appears that B-SWCNT cells gained apoptotic resistance in the intrinsic apoptotic pathway via the pAkt/MDM2/p53/Bcl-2 signaling axis. With the decreased pro-apoptotic BAX and PMAIP1 (Table 3, networks 3 and 9) expression downstream of pAkt, and our elevated levels of Bcl-2 protein (Western blot data), we postulated this anti-apoptotic signaling mechanism in the intrinsic apoptosis pathway of B-SWCNT cells.
Several DEGs in other signaling pathways potentially contributed to B-SWCNT malignant phenotype. Increased expression of RASGRP1 and RasGEF (ARHGEF9) may contribute to Ras-ERK signaling that regulates cyclin expression and thus cell cycle progression (Table 3, network 1). Increased SHC and Rho (RHOB) expression may also activate Cdc42, a key cell cycle regulator, and contribute to PI3K/Akt activation. It was reported that smoke-concentrated medium exposure induced Cdc42 translocation in human bronchial epithelial cells that may contribute to lung carcinogenesis . Downregulation of Dsh (DVL1; Table 3, network 3) may remove its blockage of Notch1, which could trigger notch signaling-mediated tumorigenesis  and p53-dependent carcinogenesis . Decreased GSK3A and GSK3β expression (Table 3, networks 3 and 9) potentially enhanced β-catenin tumor promotion and increased Cyclin D protein levels leading to enhanced G1 phase transition and accelerated cell proliferation. It is believed that GSK-3β regulates Wnt signaling pathway and its aberrant activation often results in tumor formation . In this study, downregulated GSK3β in APC complex (a prominent anti-lung cancer mediator), at the transcript or post-transcriptional level via pAkt activation, may lead to dysfunctional β-catenin and thus potential elevated Wnt signaling (Figure 4). Lastly, downregulation of JNK, a MAPK family kinase, does not seem to support Fos/Jun-mediated transcription in SWCNT carcinogenic potential. However, decreased JNK can result in decreased phosphorylation of p53 Ser15, resulting in MDM2-mediated p53 degradation .
In conclusion, we observed that long-term exposure in vitro (6 months) to SWCNT caused malignant transformation of human lung epithelial cells as suggested by our genome-wide expression data. Several genes involved in the apoptosis, cell cycle control, and oncogenic development were altered by the SWCNT treatment with the immune response representing the most altered biological process. The oncogenic phenotype of transformed B-SWCNT cells  may be explained, among other possible mechanisms, by active pAkt signaling, destabilized p53, increased expression of Ras family genes, Dsh-mediated Notch1, overexpressed Bcl-2 and downregulation of the anti-apoptotic genes BAX and PMAIP1. The development of the sub-chronic in vitro exposure model coupled with global transcriptome analysis described in this study could aid in the investigation of the potential mechanisms of CNT carcinogenesis, while the genome-wide expression approach can provide new insights into the genes involved in the carcinogenic process .
single-walled carbon nanotubes
transformed BEAS-2B cells following SWCNT exposure
passage matched, unexposed BEAS-2B cells
differentially expressed genes
Ingenuity Pathway Analysis
reactive oxygen species
Competing interests The authors declare that they have no competing interests.
Authors’ contributions DC performed data analysis, annotation, and manuscript preparation. TS did the experimental design and manuscript preparation. SL and YL performed the experiment and data collection. YR did the experimental design, interpretation, and manuscript preparation. LW was the study principle investigator for experimental design, data collection, analysis, interpretation and did the manuscript preparation. All authors contributed to preparation of this manuscript, comments and proof reading. All authors read and approved the final manuscript.
This work was partially supported by the National Institute for Occupational Safety and Health, NORA/NTRC and by grants from National Institutes of Health (R01-HL095579 and R01-ES022968) and National Science Foundation (EPS-1003907 and CBET-1434503). Authors thank Dr. Robert Oster for helpful discussion and Kai Wu at Wayne State University for the Western blot data analysis assistance.
The findings and conclusions in this report are those of the authors and do not necessarily represent the views of the National Institute for Occupational Safety and Health. Mention of brand name does not constitute product endorsement.
Zhang, Z; Yang, X; Zhang, Y; Zeng, B; Wang, S; Zhu, T: Delivery of telomerase reverse transcriptase small interfering RNA in complex with positively charged single-walled carbon nanotubes suppresses tumor growth, Clin cancer Res., 2006, 12, S. 4933–9CrossrefPubMed
Roco, MC: Environmentally responsible development of nanotechnology, Environ Sci Technol., 2005, 39, S. 106A–12Crossref
Mercer, RR; Hubbs, AF; Scabilloni, JF; Wang, L; Battelli, LA; Schwegler-Berry, D: Distribution and persistence of pleural penetrations by multi-walled carbon nanotubes, Part Fibre Toxicol., 2010, 7, S. 28CrossrefPubMed
Mercer, RR; Scabilloni, J; Wang, L; Kisin, E; Murray, AR; Schwegler-Berry, D: Alteration of deposition pattern and pulmonary response as a result of improved dispersion of aspirated single-walled carbon nanotubes in a mouse model, Am J Physiol Lung Cell Mol Physiol., 2008, 294, S. L87–97CrossrefPubMed
Shvedova, AA; Kisin, E; Murray, AR; Johnson, VJ; Gorelik, O; Arepalli, S: Inhalation vs. aspiration of single-walled carbon nanotubes in C57BL/6 mice: inflammation, fibrosis, oxidative stress, and mutagenesis, Am J Physiol Lung Cell Mol Physiol, 2008, 295, S. L552–65CrossrefPubMed
Wang, L; Mercer, RR; Rojanasakul, Y; Qiu, A; Lu, Y; Scabilloni, JF: Direct fibrogenic effects of dispersed single-walled carbon nanotubes on human lung fibroblasts, J Toxic Environ Health A., 2010, 73, S. 410–22Crossref
Sargent, LM; Porter, DW; Staska, LM; Hubbs, AF; Lowry, DT; Battelli, L: Promotion of lung adenocarcinoma following inhalation exposure to multi-walled carbon nanotubes, Part Fibre Toxicol., 2014, 11, S. 3CrossrefPubMed
Wang, L; Luanpitpong, S; Castranova, V; Tse, W; Lu, Y; Pongrakhananon, V: Carbon nanotubes induce malignant transformation and tumorigenesis of human lung epithelial cells, Nano Lett., 2011, 11, S. 2796–803CrossrefPubMed
Folkmann, JK; Risom, L; Jacobsen, NR; Wallin, H; Loft, S; Moller, P: Oxidatively damaged DNA in rats exposed by oral gavage to C60 fullerenes and single-walled carbon nanotubes, Environ Health Perspect., 2009, 117, S. 703–8CrossrefPubMed
Lindberg, HK; Falck, GC; Suhonen, S; Vippola, M; Vanhala, E; Catalan, J: Genotoxicity of nanomaterials: DNA damage and micronuclei induced by carbon nanotubes and graphite nanofibres in human bronchial epithelial cells in vitro, Toxicol Lett., 2009, 186, S. 166–73CrossrefPubMed
Yamashita, K; Yoshioka, Y; Higashisaka, K; Morishita, Y; Yoshida, T; Fujimura, M: Carbon nanotubes elicit DNA damage and inflammatory response relative to their size and shape, Inflammation., 2010, 33, S. 276–80CrossrefPubMed
Pacurari, M; Yin, XJ; Zhao, J; Ding, M; Leonard, SS; Schwegler-Berry, D: Raw single-wall carbon nanotubes induce oxidative stress and activate MAPKs, AP-1, NF-kappaB, and Akt in normal and malignant human mesothelial cells, Environ Health Perspect., 2008, 116, S. 1211–7CrossrefPubMed
Poland, CA; Duffin, R; Kinloch, I; Maynard, A; Wallace, WA; Seaton, A: Carbon nanotubes introduced into the abdominal cavity of mice show asbestos-like pathogenicity in a pilot study, Nat Nanotechnol., 2008, 3, S. 423–8CrossrefPubMed
Takagi, A; Hirose, A; Nishimura, T; Fukumori, N; Ogata, A; Ohashi, N: Induction of mesothelioma in p53+/− mouse by intraperitoneal application of multi-wall carbon nanotube, J Toxicol Sci., 2008, 33, S. 105–16CrossrefPubMed
Azad, N; Iyer, AK; Wang, L; Liu, Y; Lu, Y; Rojanasakul, Y: Reactive oxygen species-mediated p38 MAPK regulates carbon nanotube-induced fibrogenic and angiogenic responses, Nanotoxicology., 2013, 7, S. 157–68CrossrefPubMed
Snyder-Talkington, BN; Pacurari, M; Dong, C; Leonard, SS; Schwegler-Berry, D; Castranova, V: Systematic analysis of multiwalled carbon nanotube-induced cellular signaling and gene expression in human small airway epithelial cells, Toxicol Sci., 2013, 133, S. 79–89CrossrefPubMed
Lohcharoenkal, W; Wang, L; Stueckle, TA; Dinu, CZ; Castranova, V; Liu, Y: Chronic exposure to carbon nanotubes induces invasion of human mesothelial cells through matrix metalloproteinase-2, ACS Nano., 2013, 7, S. 7711–23CrossrefPubMed
Nagai, H; Okazaki, Y; Chew, SH; Misawa, N; Yamashita, Y; Akatsuka, S: Diameter and rigidity of multiwalled carbon nanotubes are critical factors in mesothelial injury and carcinogenesis, Proc Natl Acad Sci U S A., 2011, 108, S. E1330–8CrossrefPubMed
Shvedova, AA; Yanamala, N; Kisin, ER; Tkach, AV; Murray, AR; Hubbs, A: Long-term effects of carbon containing engineered nanomaterials and asbestos in the lung: one year postexposure comparisons, Am J Physiol Lung Cell Mol Physiol., 2014, 306, S. L170–82CrossrefPubMed
Wang, L; Stueckle, TA; Mishra, A; Derk, R; Meighan, T; Castranova, V: Neoplastic-like transformation effect of single-walled and multi-walled carbon nanotubes compared to asbestos on human lung small airway epithelial cells, Nanotoxicology., 2014, 8, S. 485–507CrossrefPubMed
Teeguarden, JG; Webb-Robertson, BJ; Waters, KM; Murray, AR; Kisin, ER; Varnum, SM: Comparative proteomics and pulmonary toxicity of instilled single-walled carbon nanotubes, crocidolite asbestos, and ultrafine carbon black in mice, Toxicol Sci., 2011, 120, S. 123–35CrossrefPubMed
Pacurari, M; Qian, Y; Porter, DW; Wolfarth, M; Wan, Y; Luo, D: Multi-walled carbon nanotube-induced gene expression in the mouse lung: association with lung pathology, Toxicol Appl Pharmacol., 2011, 255, S. 18–31CrossrefPubMed
Wu, P; Yuan, SS; Ho, CC; Hsieh, WY; Hong, QS; Yu, SL: Focal amplification of HOXD-harboring chromosome region is implicated in multiple-walled carbon nanotubes-induced carcinogenicity, Nano Lett., 2013, 13, S. 4632–41CrossrefPubMed
Donaldson, K; Murphy, FA; Duffin, R; Poland, CA: Asbestos, carbon nanotubes and the pleural mesothelium: a review of the hypothesis regarding the role of long fibre retention in the parietal pleura, inflammation and mesothelioma, Toxicol Sci., 2010, 7, S. 5
Gerwin, BI; Spillare, E; Forrester, K; Lehman, TA; Kispert, J; Welsh, JA: Mutant p53 can induce tumorigenic conversion of human bronchial epithelial cells and reduce their responsiveness to a negative growth factor, transforming growth factor beta 1, Proc Natl Acad Sci U S A., 1992, 89, S. 2759–63CrossrefPubMed
Hussain, SP; Amstad, P; Raja, K; Sawyer, M; Hofseth, L; Shields, PG: Mutability of p53 hotspot codons to benzo(a)pyrene diol epoxide (BPDE) and the frequency of p53 mutations in nontumorous human lung, Cancer Res., 2001, 61, S. 6350–5PubMed
Liao, WT; Lin, P; Cheng, TS; Yu, HS; Chang, LW: Arsenic promotes centrosome abnormalities and cell colony formation in p53 compromised human lung cells, Toxicol Appl Pharmacol., 2007, 225, S. 162–70CrossrefPubMed
Stueckle, TA; Lu, Y; Davis, ME; Wang, L; Jiang, BH; Holaskova, I: Chronic occupational exposure to arsenic induces carcinogenic gene signaling networks and neoplastic transformation in human lung epithelial cells, Toxicol Appl Pharmacol., 2012, 261, S. 204–16CrossrefPubMed
Harris, CC: p53 tumor suppressor gene: at the crossroads of molecular carcinogenesis, molecular epidemiology, and cancer risk assessment, Environ Health Perspect., 1996, 104, Suppl 3, S. 435–9CrossrefPubMed
Herzog, E; Byrne, HJ; Casey, A; Davoren, M; Lenz, AG; Maier, KL: SWCNT suppress inflammatory mediator responses in human lung epithelium in vitro, Toxicol Appl Pharmacol., 2009, 234, S. 378–90CrossrefPubMed
Turabekova, M; Rasulev, B; Theodore, M; Jackman, J; Leszczynska, D; Leszczynski, J: Immunotoxicity of nanoparticles: a computational study suggests that CNTs and C60 fullerenes might be recognized as pathogens by Toll-like receptors, Nanoscale., 2014, 6, S. 3488–95CrossrefPubMed
Donaldson, K; Stone, V; Seaton, A; Tran, L; Aitken, R; Poland, C: Re: Induction of mesothelioma in p53+/− mouse by intraperitoneal application of multi-wall carbon nanotube, J Toxicol Sci., 2008, 33, S. 385CrossrefPubMed
Nagai, H; Toyokuni, S: Differences and similarities between carbon nanotubes and asbestos fibers during mesothelial carcinogenesis: shedding light on fiber entry mechanism, Cancer Sci., 2012, 103, S. 1378–90CrossrefPubMed
Huang, X; Zhang, F; Sun, X; Choi, KY; Niu, G; Zhang, G: The genotype-dependent influence of functionalized multiwalled carbon nanotubes on fetal development, Biomaterials., 2014, 35, S. 856–65CrossrefPubMed
Zhang, L; Gallup, M; Zlock, L; Finkbeiner, WE; McNamara, NA: Rac1 and Cdc42 differentially modulate cigarette smoke-induced airway cell migration through p120-catenin-dependent and -independent pathways, Am J Pathol., 2013, 182, S. 1986–95CrossrefPubMed
Pez, F; Lopez, A; Kim, M; Wands, JR; Caron de Fromentel, C; Merle, P: Wnt signaling and hepatocarcinogenesis: molecular targets for the development of innovative anticancer drugs, J Hepatol, 2013, 59, S. 1107–17CrossrefPubMed
Licciulli, S; Avila, JL; Hanlon, L; Troutman, S; Cesaroni, M; Kota, S: Notch1 is required for Kras-induced lung adenocarcinoma and controls tumor cell survival via p53, Cancer Res., 2013, 73, S. 5974–84CrossrefPubMed