Abstract
Abstract
Chronic exposure to cigarette smoke markedly increases the risk for lung cancer. Regulation of gene expression at the post-transcriptional level by miRNAs influences a variety of cancer-related interactomes. Yet, relatively little is known on the effects of long-term cigarette smoke exposure on miRNA expression and gene regulation. NCI-H292 (H292) is a cell line sensitive to cigarette smoke with mucoepidermoid characteristics in culture. We report, in this study, original observations on long-term (12 months) cigarette smoke effects in the H292 cell line, using microarray-based miRNA expression profiling, and stable isotopic labeling with amino acids in cell culture-based quantitative proteomic analysis. We identified 112 upregulated and 147 downregulated miRNAs (by twofold) in cigarette smoke-treated H292 cells. The liquid chromatography–tandem mass spectrometry analysis identified 3,959 proteins, of which, 303 proteins were overexpressed and 112 proteins downregulated (by twofold). We observed 39 miRNA target pairs (proven targets) that were differentially expressed in response to chronic cigarette smoke exposure. Gene ontology analysis of the target proteins revealed enrichment of proteins in biological processes driving metabolism, cell communication, and nucleic acid metabolism. Pathway analysis revealed the enrichment of phagosome maturation, antigen presentation pathway, nuclear factor erythroid 2-related factor 2-mediated oxidative stress response, and cholesterol biosynthesis pathways in cigarette smoke-exposed cells. In conclusion, this report makes an important contribution to knowledge on molecular changes in a lung cell line in response to long term cigarette smoke exposure. The findings might inform future strategies for drug target, biomarker and diagnostics innovation in lung cancer, and clinical oncology. These observations also call for further research on the extent to which continuing or stopping cigarette smoking in patients diagnosed with lung cancer translates into molecular and clinical outcomes.
Introduction
S
Patients with COPD have increased incidence of lung cancer (Tockman et al., 1987). Even the low-grade COPD without noticeable airflow obstruction is associated with elevated risk of lung cancer (de Torres et al., 2007). Cigarette smoking has also been found to induce a number of genetic and molecular changes in the respiratory tract, including cellular atypia (Franklin et al., 1997) and loss of heterozygosity (Powell et al., 1999; Wistuba et al., 1997), which can occur in cytologically normal airway epithelium (Thiberville et al., 1995; Wistuba et al., 1997). Yet, molecular alterations induced by smoking and the factors that impact the survival rates remain elusive.
Accumulation of both genetic and epigenetic alterations results in sequential morphological changes leading to tumorigenesis (Osada and Takahashi, 2002). Regulation of gene expression at the posttranscriptional level by miRNAs influences a variety of pathways, including developmental and oncogenic pathways (Calin and Croce, 2006; He et al., 2007; Hu et al., 2010; Johnson et al., 2005).
The crucial role of miRNAs in various cellular processes such as differentiation, cell growth, proliferation, and apoptosis has been established by several groups (Kent and Mendell, 2006; Xu et al., 2011). Dysregulation of miRNAs and their targets is associated with hallmarks of cancer, including cytoskeletal remodeling, proliferation, and migration (Huang et al., 2008; Kobayashi et al., 2008; Lee et al., 2005). Distinct miRNA expression profiles have been observed for different cancers (Barker et al., 2009). miRNAs ultimately affect protein expression either by mediating degradation of mRNA or by inhibiting translation (Jacobsen et al., 2013; Yang et al., 2009). Since dysregulation of miRNAs can affect multiple targets leading to diverse functional consequences, it is imperative to study the targets of the dysregulated miRNAs in any disease.
miRNA expression profile of lung cancer has been previously characterized using microarrays (Yanaihara et al., 2006) as well as next-generation sequencing technologies (Ma et al., 2014; Wang et al., 2015). Cigarette smoke has been shown to alter the expression of several genes in various biological models, including mice (Izzotti et al., 2004), rats (Izzotti et al., 2005), and airway epithelial cells from smokers (Spira et al., 2004). In addition, miRNA dysregulation has also been reported in bronchial epithelium of smokers (Schembri et al., 2009) and animal models exposed to cigarette smoke (Izzotti et al., 2009). Recently, lung cancer cells exposed to cigarette smoke revealed the dysregulation of several genes in response to particulate matter and vapor phase components of cigarette smoke (Sekine et al., 2015).
One of the genes that showed enhanced expression in response to cigarette smoke includes HMOX1, an antioxidant gene regulated by the transcription factor nuclear factor erythroid 2-related factor 2 (NRF2). In addition, Nrf2 signaling and DNA damage response pathways were found to be enriched in response to cigarette smoke. Integration of miRNA expression and proteomic data in cigarette smoke-exposed lung cells could lead to the identification of potential signaling mechanisms affected by cigarette smoke. NCI-H292 (H292) is a bronchial epithelial cell line derived from a mucoepidermoid carcinoma established from a nonsmoker. It is one of the widely used cell lines used to study the effect of cigarette smoke as it is sensitive to cigarette smoke and retain their mucoepidermoid characteristics in culture (Newland and Richter, 2008). Several studies have reported the effect of cigarette smoke on NCI-H292 cells to be similar to that of primary human bronchial epithelial cells and the airway epithelium in vivo (Baginski et al., 2006; Luppi et al., 2005; Newland and Richter, 2008; Shao et al., 2004; Thorne et al., 2009).
In vitro cellular models where human lung cells have been exposed to varied doses of cigarette smoke have served as useful tool to interpret disease progression and understand the molecular mechanisms associated with them. Even though these studies have identified few molecular mechanisms by which cigarette smoke may exert its effects, there has been no study on the chronic effect of cigarette smoke in lung cells. It is known that the carcinogenic effect of cigarette smoke in lung cancer is through chronic exposure and not by acute exposure.
In this study, we investigated chronic effects of cigarette smoke on miRNA expression as well as proteomic changes in H292 cells chronically (12 months) exposed to cigarette smoke using microarray and mass spectrometry (MS) approaches.
Materials and Methods
Cell culture and stable isotopic labeling with amino acids in cell culture labeling
H292 cells were obtained from the American Type Culture Collection and grown in Dulbecco's modified Eagle's medium (DMEM) medium (Invitrogen) supplemented with 10% fetal bovine serum (FBS), and 100 U/mL penicillin and 100 μg/mL streptomycin in a humidified incubator at 37°C with 5% CO2. These cells were and adapted to stable isotopic labeling with amino acids in cell culture (SILAC) media as previously described (Harsha et al., 2008). Briefly, the cells were maintained in DMEM medium without lysine and arginine supplemented with 10% FBS, 2 mM
To study the chronic effect of cigarette smoke condensate (CSC; Murty Pharmaceuticals, Inc.), H292 cells were subjected to chronic treatment with 0.1% CSC for 12 months (Chang et al., 2011). H292 cells were grown in a smoke-dedicated incubator. Parental cells were grown in regular CO2 incubator. Cells that were grown and passaged in the smoke-dedicated incubator exposed to CSC were labeled as H292-smoke or termed as cigarette smoke-treated cells. Henceforth, the H292 cells exposed to CSC will be referred to as H292-smoke and untreated cells as H292-parental.
Trypsin digestion and basic pH reverse phase liquid chromatography
The CSC-treated and CSC-untreated H292 cells were washed with cold phosphate-buffered saline and lysed in lysis buffer (20 mM HEPES, pH 8.0, 9 M urea, 1 mM sodium orthovanadate, 2.5 mM sodium pyrophosphate, and 1 mM β-glycerophosphate), sonicated, and centrifuged at 16,000 g at 15°C for 20 min. The protein concentration was determined using the BCA assay (Pierce). Equal amounts of protein from CSC-treated and CSC-untreated H292 cells were reduced with 5 mM dithiothreitol and incubated at 60°C for 20 min. These were further alkylated using 10 mM iodoacetamide for 10 min at room temperature in the dark. The samples were diluted such that urea was <2 M with 20 mM HEPES, pH 8.0.
The samples were subjected to digestion with TPCK-treated trypsin (Worthington Biochemical Corporation) overnight at room temperature. The peptide mixture was acidified using 1% trifluoroacetic acid and desalted using C18 Sep-Pak cartridge (Waters; Cat#WAT051910). The extracted peptides were lyophilized and subjected to basic pH reverse phase chromatography (bRPLC). The samples were reconstituted in bRPLC solvent A [7mM triethyl ammonium bicarbonate (TEABC), pH 9] and separated on a XBridge BEH C18 column (Waters) using a linear increase in gradient from 5% to 100% of 7 mM TEABC with 90% acetonitrile (pH 9) over 30 min. A total of 96 fractions were initially collected and later concatenated to 24 fractions and dried using SpeedVac.
bRPLC fractionation and liquid chromatography–tandem mass spectrometry analysis
The samples were analyzed on an LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) interfaced with an Easy-nLC II nanoflow liquid chromatography system (Thermo Scientific). The peptides were reconstituted in 0.1% formic acid and loaded onto a 75 μm × 2 cm Magic C18 AQ (Michrom Bioresources, Inc.), 5 μm, 120 Å trap column. Peptides were resolved on a 75 μm × 20 cm analytical column at a flow rate of 350 nL/min using a linear gradient of 8–30% solvent B (0.1% formic acid in 95% acetonitrile) over 75 min. The total run time, including sample loading and column reconditioning, was 95 min. Data-dependent acquisition was carried out with full scans in 350–1800 m/z range with a mass resolution of 60,000 at 400 m/z. Fifteen most intense precursor ions from a survey scan were selected for MS/MS, fragmented using higher energy collisional dissociation (HCD) with 39% normalized collision energy, and detected at a mass resolution of 15,000 at 400 m/z. Dynamic exclusion was set for 40 sec with a 7 ppm mass window. Internal calibration was carried out using lock mass option (m/z 445.1200025) from ambient air.
Proteomics data analysis
MS-derived data were searched against human RefSeq database (version 65 containing 34,454 entries along with common contaminants) using SequestHT and Mascot (version 2.2) search algorithms through Proteome Discoverer (Thermo Scientific; version 1.4.1.14). The search parameters included carbamidomethylation of cysteine as a fixed modification and N-terminal acetylation, oxidation of methionine, and SILAC labeling (13C6) at lysine and arginine (6.02013 Da) as variable modifications. MS/MS spectra were searched with a precursor mass tolerance of 20 ppm and a fragment mass tolerance of 0.1 Da. Trypsin was specified as the protease and a maximum of two missed cleavages were allowed. The data were searched against decoy database and a 1% false discovery rate was set at the peptide level. Proteins with ratios ≥2-fold were considered to be upregulated where as those ≤0.5-fold were considered to be downregulated.
miRNA microarray analysis
The miRNA microarray profiling was carried out using Affymetrix GeneChip miRNA Arrays (Affymetrix) as per the manufacturer's instructions. Total RNA was extracted from H292 cells and quality assessment of the samples was done using the TaqMan assay. The miRNA targets were biotinylated using an Asuragen developed direct labeling procedure. Once labeled, the miRNA targets were hybridized overnight onto the microarrays following which the arrays were washed and stained using Streptavidin-Phycoerythrin. Arrays were processed using the gene chip operating software (GCOS). Data extraction was from the Affymetrix miRNA QC tool using robust multi-array average (RMA) background correction and median polish summarization. Candidate miRNAs that were increased or decreased in expression with fold changes ≥2-fold following treatment with CSC were considered to be differentially regulated.
Bioinformatics analysis
The identified proteins were sorted with a gene ontology-based functional classification by Human Protein Reference Database (HPRD 1 ) (Keshava Prasad et al., 2009a, 2009b). Molecule class, molecular function, biological process, subcellular localization, tissue expression, and domain information were fetched from HPRD for all proteins. FunRich tool was used for functional enrichment of the 415 dysregulated proteins using the HPRD annotation as backend database (Pathan et al., 2015). The potential targets of differentially expressed miRNAs were predicted using TargetScan algorithm 2 (Lewis et al., 2005). These were then compared with the proteomic data to infer the expression of target proteins under chronic smoke exposure. Pathway analysis of the differentially expressed proteins was carried out using QIAGEN's Ingenuity Pathway Analysis (IPA®; QIAGEN Redwood City 3 ).
Identification of dysregulated miRNAs from The Cancer Genome Atlas
miRNA expression of NSCLC (squamous cell carcinoma and adenocarcinoma) and its corresponding clinical information for Illumina GA and HiSeq platform were obtained from The Cancer Genome Atlas (TCGA 4 ). Data used were obtained from the study entitled “miRNA Analysis of TCGA LUSC and LUAD Samples using Illumina.” Patient barcodes were segregated into “smokers” and “nonsmokers” based on the clinical information provided. Cases where no information was available were not used.
Based on the clinical information provided by TCGA, the expression profile of miRNA for each case was obtained using both Illumina platforms. In addition, an in-house PERL script was used to create a matrix of the miRNA expression of each patient across each miRNA for tumor smoker versus nonsmoker from the matrix; a median value of the expression for each miRNA was calculated across patient samples. Differential expression of miRNA is further calculated by taking the ratios of the expression levels of smoker versus nonsmoker.
Data availability
The MS proteomics data have been deposited to the ProteomeXchange Consortium 5 (Vizcaino et al., 2014) by the PRIDE partner repository with the dataset identifier PXD003086.
Results
Dysregulated miRNAs
We employed a microarray-based approach to carry out miRNA expression analysis of H292 cells chronically treated with cigarette smoke. A schematic representation of the work flow is presented in Figure 1 (left panel). Our analysis led to the identification of 834 miRNAs among which 259 were dysregulated (≥2-fold) in the lung cells treated with cigarette smoke. A total of 112 miRNAs were upregulated and 147 were downregulated in the cigarette smoke-treated cells. A partial list of dysregulated miRNAs is provided in Table 1. A complete list of miRNAs along with corresponding fold changes in smoke-treated cells is provided in Supplementary Table S1.

Experimental workflow employed to profile miRNA expression in smoke-exposed H292 cells. RNA was isolated from H292 cells, chronically treated with CSC. Samples were processed further depending on the integrity of the isolated RNA and labeling was carried out. The samples were hybridized on miRNA microarrays. Data were processed further to identify dysregulated miRNAs. Experimental workflow employed to analyze the proteome of smoke-exposed H292 cells. H292 cells were chronically treated with CSC. The cells were adapted to stable isotopic labeling with amino acids in cell culture media. The cells treated with CSC were grown in a medium containing normal arginine and lysine. The untreated cells were grown in a medium containing heavy arginine and lysine. Equal amounts of proteins from both treated and untreated samples were digested. The peptides were fractionated by strong cation exchange chromatography. The fractionated samples were subjected to LC-MS/MS analysis on a LTQ-Orbitrap Velos mass spectrometry. The data were searched and quantitated using the Sequest and Mascot search algorithms. CSC, cigarette smoke condensate; LC-MS/MS, liquid chromatography–tandem mass spectrometry.
Proteins differentially expressed
To understand how the dysregulated miRNAs affect the expression of their targets, we employed SILAC-based proteomic strategy. Liquid chromatography–tandem mass spectrometry-based quantitative proteomics of H292 parental and H292 smoke cells led to the identification of 3959 proteins. The proteomic workflow used for this study is outlined in Figure 1 (right panel). Among the differentially expressed proteins, 302 proteins were overexpressed (≥2-fold) and 112 proteins were found to be downregulated (≤0.5-fold) in the cigarette smoke-treated H292 cells. A partial list of proteins is provided in Table 2. The complete list of identified proteins and peptides is provided in Supplementary Table S2. MS/MS spectra of a subset of these candidates are provided in Figure 2 and Supplementary Figure S1.

Representative MS/MS spectra depicting the overexpression of
Correlation of miRNA with protein expression
TargetScan (Friedman et al., 2009) and mirTarBase (Chou et al., 2016) were used to determine the targets of dysregulated miRNAs identified for H292 cells treated with cigarette smoke. Only experimentally proven targets were considered for further downstream analysis. As mentioned earlier, 259 miRNAs were dysregulated in the smoke-exposed H292 cells. Among the dysregulated miRNAs, 202/259 miRNAs were dysregulated and identified in our data set only (Fig. 3). Of the 202 miRNAs that were dysregulated and unique to our data set, 132 were present in the TargetScan and mirTarBase databases and were utilized for proven target prediction (Supplementary Table S3). This led to the identification of 6853 targets for H292 data sets. 2431/6853 were identified in our proteomic data.

Bioinformatics workflow. Pipeline depicting the analysis of dysregulated miRNA in H292 smoke-exposed cells and non-small-cell lung cancer The Cancer Genome Atlas data sets.
Supplementary Table S3 lists the dysregulated miRNAs in the smoke-exposed H292 cells and their targets. The overall significant directional changes in miRNA and protein expression were found for 39 miRNAs (Fig. 4), which were identified in our data set. A partial list of miRNAs and their target proteins with directional changes is provided in Table 3. A complete list of miRNAs and their targets showing reciprocal expression is provided as Supplementary Table S4A and B.

Bubbleplot representing miRNAs and their protein targets, which shows reciprocal expression in H292 smoke-exposed cells. Significant dysregulation in inverse direction among miRNAs and their respective target proteins was observed in 39 miRNAs. Size of the bubble represents fold change.
Bioinformatics analysis
The identified proteins were classified based on their molecule class, molecular function, biological process, subcellular localization, and tissue expression using HPRD (Keshava Prasad et al., 2009a, 2009b) 6 . In addition, information pertaining to the presence of signal peptides and domains in identified proteins was also obtained from HPRD. Pathway analysis of proteins that were differentially expressed in response to cigarette smoke exposure was carried out using IPA. Pathway analysis of upregulated proteins led to the identification of several enriched pathways, including phagosome maturation, antigen presentation pathway, NRF2-mediated oxidative stress response, and cholesterol biosynthesis pathways. Similar analysis using downregulated proteins identified the enrichment of retionic acid receptor (RAR) activation pathway.
In addition, we carried out the analysis of differentially expressed proteins using FunRich tool. Subcellular localization-based classification revealed that 47% of the identified proteins localized to cytoplasm, 30% of the identified proteins localized to the lysosome, and 21% of the proteins were present in the endoplasmic reticulum. Further classification based on the biological process showed an enrichment of proteins involved in metabolism (20%), energy pathways (19%), and transport (11%) processes. The gene ontology-based functional annotation of identified proteins is depicted as Figure 5, which shows significant enrichment of lysosomal proteins (30%) in the smoke-exposed lung cells. Other components enriched under gene ontology included endoplasmic reticulum (20%), exosomes (27%), and mitochondria (16%).

Gene ontology-based functional annotation of proteins identified in H292 cell lines.
Comparison to the TCGA datasets of dysregulated miRNAs in NSCLC
Apart from studying the miRNA dysregulation in H292 cells in response to cigarette smoke, we also compared our data set with NSCLC miRNA data set of TCGA. We fetched miRNA-seq raw data from TCGA repository for NSCLC. miRNA expression was calculated in smokers compared to nonsmokers. In total, we identified 933 miRNAs (Fig. 3). We performed an F-test between smokers and nonsmokers using p-value to select all the potential miRNAs expressed. We found that 606 miRNAs had a significant p-value (≤0.05) and 263/606 miRNAs were dysregulated in TCGA data sets (Supplementary Fig. S2A, B). We then compared the miRNAs dysregulated in TCGA dataset with our miRNA data pertaining to the H292 cells treated chronically with cigarette smoke.
Fifty-seven miRNAs were dysregulated and common between H292 (smoke) and TCGA data set (Fig. 3 and Supplementary Fig. S2B); among these, 26 miRNAs showed concordance expression between the two data sets (Supplementary Fig. S3A, B). Supplementary Table S5A and B lists the dysregulated miRNAs in NSCLC (adenocarcinoma and squamous cell carcinoma, respectively) in TCGA data sets with their targets.
Discussion
miRNAs are one of the major regulators of gene expression. They are known to regulate several key cellular processes, including proliferation, differentiation, and apoptosis (Bueno and Malumbres, 2011). Dysregulation of miRNAs have been widely associated with the development of various cancers, including lung cancer (Guz et al., 2014; Hayes et al., 2014). Cigarette smoke has been identified as one of the risk factors associated with the development of lung cancer. Several research groups have shown in both human and rat cell models that exposure of cigarette smoke resulted in the dysregulation of miRNAs (Izzotti et al., 2009, 2010; Schembri et al., 2009; Vucic et al., 2014). However, the exact mechanisms by which cigarette smoke induces the dysregulation of miRNAs are not fully understood.
Against the above overarching scientific context, we carried out microarray-based miRNA expression profiling of cigarette smoke-treated H292 cells. In addition, we also carried out a SILAC-based proteomic analysis to identify differentially expressed proteins in H292 cells after exposure to cigarette smoke. The identification of reciprocal expression of miRNAs and their targets in H292 cells exposed to cigarette smoke provides a unique molecular basis to further explore the biological role of miRNAs in NSCLC. Integrative analyses identified both predicted and experimentally validated miRNA-target pairs, but only experimentally validated targets were considered.
We identified several miRNAs previously reported to be dysregulated in lung cancer corroborating our approach. miR-124 has been reported to be dysregulated in both highly invasive subcell lines and NSCLC (Sun et al., 2015). miR-124 (miR-124-3p) was found to be upregulated in our data set by about sevenfold in response to cigarette smoke. We observed an overexpression of miR-21 (hsa-miR-21-3p) and miR-192 (hsa-miR-192-5p) in the H292 smoke cells. Expression levels of serum miR-21 and miR-192 were previously reported to be significantly higher in stage I NSCLC patients when compared to controls (Qi et al., 2014b).
The reduced expression of miR-125a (hsa-miR-125a-3p; identified as downregulated in our study) has been shown to be strongly associated with NSCLC dedifferentiation in tissue and serum samples (Zhu et al., 2014). It has also been reported to be a key miRNA responsible for the activation of p53 and induction of apoptosis in human lung cancer cells (Jiang et al., 2013). Low expression of miR-374a (hsa-miR-374a-5p; downregulated by 5.4-fold in this study) has been associated with poor patient survival in early-stage NSCLC (Vosa et al., 2011). It has also been known to promote metastasis in breast cancer cells through the Wnt/β-catenin signaling cascade (Cai et al., 2013). Xia et al. (2014) have reported involvement of miR-204 (hsa-miR-204-5p) in NSCLC development as its expression was decreased in tumor samples.
Furthermore, overexpression of miR-204 led to reduced cell proliferation and invasion, increased expression of E-cadherin, and decreased expression of N-cadherin and Vimentin (Xia et al., 2014). In our data set, miR-204 (hsa-miR-204-5p) was found to be downregulated by about 4-fold in response to cigarette smoke. Similarly, miR-101, miR-125a, miR-140, miR-145, and miR-27b, which were found to be downregulated in the H292 cells chronically treated with cigarette smoke, have been reported to be downregulated in lung adenocarcinoma (Yanaihara et al., 2006).
Apart from miRNAs reported to be dysregulated in NSCLC, we also identified several miRNAs reported to be downregulated in other cancers, but identified and reported for first time in NSCLC and cigarette smoke. These include miR-214, miR-215, miR-342, miR-34b, miR-378, and miR-424 in colon cancer (Gaur et al., 2007; Schetter et al., 2008), miR-194 in breast cancer (Iorio et al., 2005), miR-195 in hepatocellular carcinoma (Murakami et al., 2006), miR-152 and miR-155 in pancreatic cancer (Volinia et al., 2006), and miR-32 in uterine leiomyoma (Wang et al., 2007). A total of four differentially expressed miRNAs (miR-146, miR-27, miR-877, and miR-186) have been known to be involved in metastatic melanoma (Qi et al., 2014a). We observed a downregulation of miR-27 and miR-877 in H292 cells treated with cigarette smoke.
Consistent with the cancer biology literature, many miRNAs identified to be dysregulated in our study showed reciprocal relationship with their targets. PPP2R2A has been shown to inhibit Akt phosphorylation in a cohort of NSCLC tumors and is a target of miR-222. HMGA1 is also reported to be regulated by miR-222 (Zhang et al., 2011). We found miR-222 to be downregulated by 0.5-fold and its targets PPP2R2A (twofold) and HMGA1 (1.8-fold) were overexpressed in response to cigarette smoke.
Liu et al. (2009) showed that ectopic transfection of miR-222 reduced the expression of MMP1 and SOD2 in oral tongue squamous cell carcinoma cell lines. Our data indicate a nine and sevenfold increased expression of MMP1 and SOD2, respectively, in response to cigarette smoke. miR-204 (hsa-miR-204-5p; downregulated by fourfold), which acts as a tumor suppressor, has been reported to be downregulated in renal clear cell carcinoma. It has been reported to be regulated by von Hippel-Lindau tumor-suppressor gene through its target, MAP1LC3B (Mikhaylova et al., 2012).
miR-145 is reported to inhibit cell proliferation of human lung adenocarcinoma by negative regulation of EGFR and NUDT1 (Cho et al., 2011). We found miR-145 (hsa-miR-145-3p) downregulated by 4.6-fold and also identified its target EGFR downregulated by 0.4-fold in response to cigarette smoke. In many cancers, including glioblastoma (Chen et al., 2008a), oral squamous cell carcinoma (Reis et al., 2010), and cholangiocarcinoma (Liu et al., 2012), programmed cell death 4 (PDCD4), a tumor-suppressor gene, was shown as a target of miR-21, and its expression is downregulated by miR-21. PDCD4 was downregulated in our study, while miR-21 (hsa-miR-21-5p) was upregulated. miR-183 (hsa-miR-183-5p; upregulated by 4.6-fold) has been shown to inhibit apoptosis and promote proliferation and invasion of gastric cancer cells by targeting PDCD4 (Gu et al., 2014).
Similarly, miR-34a (hsa-miR-34a-5p), which has been previously shown to suppress cell proliferation and metastasis and targets CD44 in human renal carcinoma cells (Yu et al., 2014), was found to be upregulated in our data set. Interestingly, contrary to the existing literature, we also identified several miRNAs that were downregulated in our data set and showed reciprocal correlation with their targets. These include miR-155 (hsa-miR-155-5p), which was downregulated by about 10-fold, and its target CEBPB is upregulated by 2.25-fold in our data. Salemi et al. (2015) reported miR-155 to be upregulated and its target CEBPB as downregulated in acute myeloid leukemia.
Similarly, miR-196a was shown to be upregulated and its target S100A9 was reported as downregulated in Barrett's metaplasia-dysplasia-invasive adenocarcinoma (Maru et al., 2009), while our data revealed downregulation of miR-196a (hsa-miR-196a-5p) by about 1.5-fold and its target S100A9 upregulated by 89.25-fold. Such observations indicate that the expression of miRNAs may be either tissue specific or cigarette smoke has its own unique mechanism of regulation of miRNAs and their corresponding targets, which needs further investigation and is beyond the scope of this article.
Apart from dysregulation of miRNAs, we also focused on the regulation of proteins, which are the key targets of miRNAs. Pathway analysis using IPA indicated an enrichment of pathways, including NRF2-mediated oxidative stress response, cholesterol biosynthesis, antigen presentation, and phagosome maturation pathways. In addition, our gene ontology data indicated a significant enrichment of lysosomal proteins in the smoke-treated cells. Nrf2, a transcription factor, has been known to regulate redox homeostasis to protect the cell from oxidative damage (Fujihara et al., 2008). Chronic cigarette smoke exposure has been previously observed to induce oxidative damage in murine retinal pigmented epithelial cells (Cano et al., 2010).
In addition, Nrf2 has been known to increase sulfiredoxin-1 (Srx1) expression during cigarette smoke exposure-induced oxidative stress in murine lung cells (Singh et al., 2009), which in turn has been reported to protect cells from oxidative injury in the lungs. Nrf2-deficient mice have been found to be highly susceptible to cigarette smoke-induced emphysema (Iizuka et al., 2005). Furthermore, cigarette smoking among heavy smokers has been observed to repress the Nrf2 pathway expression in peripheral mononuclear cells through the recruitment of NFκB (Garbin et al., 2009).
Our analysis found the enrichment of cholesterol biosynthesis pathway in response to cigarette smoke exposure. Cigarette smoke exposure has been reported to impair plasma lipid profile as well as reverse cholesterol transport in murine models (Zong et al., 2015). It was found that cigarette smoke significantly decreased plasma high-density lipoprotein (HDL) cholesterol levels, while increasing plasma low-density lipoprotein (LDL) cholesterol levels in smokers (Meenakshisundaram et al., 2010). A study investigating the lipid profiles among 274 smokers found a significant rise in LDL and a fall in HDL levels. Interestingly, another study has indicated that the production of progesterone, a cholesterol derivative, was inhibited by cigarette smoke alkaloids (Gocze and Freeman, 2000). Smoking cessation has been found to improve total HDL levels (Gepner et al., 2011).
We observed enrichment of the antigen presentation pathway in response to cigarette smoke. Literature evidence indicates that cigarette smoke decreases immunoproteasome function of immune cells isolated from COPD patients and impairs immunoproteasome-specific MHC I antigen presentation (Kammerl et al., 2016). In addition, long-term cigarette smoke exposure has been shown to suppress the development of functional dendritic cells, which are antigen-presenting cells, suggesting that cigarette smoke exposure may lead to impaired immune response to infection (Givi et al., 2015).
Apart from the pathways enriched, we also observed an enrichment of lysosomal proteins in the smoke-exposed lung cells. Cancer cells are known to demonstrate alterations of the lysosomal compartment that have pro-oncogenic effects, including tumor growth, migration, invasion, and angiogenesis (Mohamed and Sloane, 2006). Furthermore, literature evidence indicates increased numbers of autophagosomes in tissues from patients with COPD and in lung epithelial cells exposed to cigarette smoke extract (Chen et al., 2008b; Kim et al., 2008). Although there is an increase in the number of autophagosomes in alveolar macrophages of smokers, there is defect or decrease in the autophagy function. This defect in autophagy function in smokers leads to mitochondrial dysfunction, increased accumulation of cigarette smoke particulates, and decreased clearance of pathogens (Monick et al., 2010).
Among the lysosomal proteins identified in our study are acid ceramidase (ASAH1) and LAMP1 and 2. ASAH1 has been reported to be associated with the cigarette smoke-mediated increase in sphingosine (Sph) synthesis (Titz et al., 2016). Furthermore, acid ceramidase (ASAH1) has been proposed as a chemotherapeutic target in NSCLC to overcome resistance to the antitumoral effect of choline kinase α inhibition as it synergistically sensitizes lung cancer cells (Ramirez de Molina et al., 2012). In addition to acid ceramidase, we also found an increased expression and enrichment of lysosome-associated membrane glycoproteins (LAMP1 and LAMP2). LAMP1 plays a decisive role in lung cancer metastasis by promoting interactions with organ extracellular matrix and basement membrane (Agarwal et al., 2015).
Furthermore, LAMP1 has been reported to be overexpressed in A549 lung cancer cells (Su et al., 2016) and known as a target of miR-373, which is downregulated in lung cancer tissues and cell lines (Seol et al., 2014). Similarly, higher expression of LAMP2 has been reported in bronchoalveolar lavage of lung cancer cells (Li et al., 2013) and its higher expression was related to high histology grade in lung cancer patients (Giatromanolaki et al., 2015). miR-487b-5p is known to regulate temozolomide resistance in lung cancer cells through LAMP2-medicated autophagy (Bao et al., 2016).
Other lysosomal proteins identified in our data include gamma-glutamyl hydrolase, cathepsin, and calnexin. Gamma-glutamyl hydrolase expression is correlated with poor prognosis in lung and breast cancers (He et al., 2004; Shubbar et al., 2013). Similarly, cathepsin B has been proposed as a potential prognostic and therapeutic marker for lung squamous cell carcinoma (Gong et al., 2013) and its expression in pulmonary adenocarcinomas and squamous cell carcinomas predicts poor prognosis (Cordes et al., 2009). Calnexin has been proposed as a diagnostic marker in lung cancer (Kobayashi et al., 2015). The recapitulation of existing literature on lung cancer by our cellular model indicates that our cellular model system is robust to study the molecular changes induced by cigarette smoke and our data indicate that cigarette smoke leads of impaired functioning of lysosomal proteins.
Conclusions and Outlook
This study used an integrated approach to analyze miRNA and protein expression profiles of non-small-cell lung cancer cells chronically exposed to cigarette smoke. Previous studies have shown that continued cigarette smoking by patients with lung cancer receiving chemoradiotherapy is associated with decreased survival (Videtic et al., 2003). We identified an altered expression of several miRNAs as well as proteins in response to cigarette smoke. Pathway analysis using our proteomics data set revealed enrichment of several pathways in response to cigarette smoke, including NRF2-mediated oxidative stress response, cholesterol biosynthesis pathways.
Cigarette smoke has been associated with defects in the autophagy/lysosomal pathway, which in turn lead to accumulation of damaged proteins and particulate matter, mitochondrial damage, and impaired innate immune responses. These results indicate that chronic exposure to cigarette smoke alters miRNA-dependent pathways in NSCLC, potentially aggravating the disease.
In summary, this report makes an important contribution to knowledge on molecular changes in a lung cell line in response to long-term cigarette smoke exposure. The findings might also inform the future strategies for molecular drug target, biomarker and diagnostics innovation in lung cancer, and clinical oncology.
Footnotes
Acknowledgments
We thank the Department of Biotechnology (DBT), Government of India, for research support to the Institute of Bioinformatics (IOB), Bangalore. IOB is supported by DBT Program Support on Neuroproteomics and infrastructure for proteomic data analysis (BT/01/COE/08/05). This work was supported by the Science and Engineering Research Board, Department of Science and Technology, Government of India, grant “Role of miRNAs in promoting lung cancer in response to cigarette smoke (SR/FT/LS-113/2011),” the NCI's Clinical Proteomic Tumor Analysis Consortium initiative (U24CA160036), and FAMRI-funded 072017_YCSA.
J.A. and A.R. are recipients of Senior Research Fellowship from the Council of Scientific and Industrial Research (CSIR), New Delhi, India. A.A.K. is a recipient of a Senior Research Fellowship from Indian Council of Medical Research (ICMR). miRNA expression of NSCLC and its corresponding clinical information for Illumina GA and HiSeq platform were obtained from The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga/dataAccessMatrix.htm) as noted under the “
” section.
Author Disclosure Statement
The authors declare that no conflicting financial interests exist.
Abbreviations Used
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
