Next Article in Journal
Seasonal Distribution and Genotyping of Antibiotic Resistant Strains of Listeria Innocua Isolated from A River Basin Categorized by ERIC-PCR
Previous Article in Journal
Foot Arch Height and Quality of Life in Adults: A Strobe Observational Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differential Effect of Smoking on Gene Expression in Head and Neck Cancer Patients

1
Department of Prosthetic Dentistry and Dental Materials, Division Dental Propaedeutics, Aesthetic, Iuliu Hatieganu University of Medicine and Pharmacy, Cluj-Napoca, 23 Marinescu Street, Cluj-Napoca 40015, Romania
2
Research Center for Functional Genomics and Translational Medicine, Iuliu Hatieganu University of Medicine and Pharmacy, 23 Marinescu Street, Cluj-Napoca 40015, Romania
3
MEDFUTURE-Research Center for Advanced Medicine, University of Medicine and Pharmacy Iuliu Hatieganu, 23 Marinescu Street, Cluj-Napoca 40015, Romania
4
Department of Medical Biology, Faculty of Medicine, Medical University-Plovdiv, 15-А Vassil Aprilov Blvd., Plovdiv 4000, Bulgaria
5
Technological Center for Emergency Medicine, 15-А Vassil Aprilov Blvd., Plovdiv 4000, Bulgaria
6
Prosthetics and Dental Materials, Faculty of Dental Medicine, Iuliu Hatieganu University of Medicine and Pharmacy, Cluj-Napoca, 32 Clinicilor Street, Cluj-Napoca 400006, Romania
7
Department of Functional Genomics and Experimental Pathology, The Oncology Institute Ion Chiricuta, Republicii 34th Street, Cluj-Napoca 400015, Romania
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(7), 1558; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15071558
Submission received: 27 April 2018 / Revised: 11 July 2018 / Accepted: 17 July 2018 / Published: 23 July 2018

Abstract

:
Smoking is a well-known behavior that has an important negative impact on human health, and is considered to be a significant factor related to the development and progression of head and neck squamous cell carcinomas (HNSCCs). Use of high-dimensional datasets to discern novel HNSCC driver genes related to smoking represents an important challenge. The Cancer Genome Atlas (TCGA) analysis was performed in three co-existing groups of HNSCC in order to assess whether gene expression landscape is affected by tobacco smoking, having quit, or non-smoking status. We identified a set of differentially expressed genes that discriminate between smokers and non-smokers or based on human papilloma virus (HPV)16 status, or the co-occurrence of these two exposome components in HNSCC. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways classification shows that most of the genes are specific to cellular metabolism, emphasizing metabolic detoxification pathways, metabolism of chemical carcinogenesis, or drug metabolism. In the case of HPV16-positive patients it has been demonstrated that the altered genes are related to cellular adhesion and inflammation. The correlation between smoking and the survival rate was not statistically significant. This emphasizes the importance of the complex environmental exposure and genetic factors in order to establish prevention assays and personalized care system for HNSCC, with the potential for being extended to other cancer types.

1. Introduction

Head and neck squamous cell carcinomas (HNSCCs) represent a preventable pathology which continues to be an important factor of morbidity with high mortality rates at global level [1,2], with over 600,000 new cases detected each year [3,4], and a mortality rate of around 50% [5]. HNSCCs have as common localizations the oral cavity, oropharynx, hypopharynx, and larynx [1,5,6].
An important aspect in prevention and treatment is related to genetic and environmental components [2,7,8]. While the acquired genetic factors cannot be controlled in their early steps of mutation accumulation, environmental exposure can significantly affect the pathogenesis and the prognosis of these patients [8,9]. The major environmental components are tobacco smoking, betel quid chewing, alcohol consumption, and poor oral hygiene and infections, alongside other specific dietary habits or specific pollutant exposure [6,7,10,11]. The totality of risk factors are integrated into the exposome [12], this being an important step in the evaluation of internal and external exposure, generally related to the co-occurrence of multiple toxic environmental agents [13]. Usually, the effect of co-occurrence is much more dramatic than that of the single exposure [13].
Microarray data are provided by large consortium programs such as The Cancer Genome Atlas (TCGA), offering new possibilities and a better understanding of the role of genes in different cancers [2,14]. Analysis of altered gene expression signatures is used in a wide range of pathologies for achieving relevant information with prognostic value [15]. The altered transcriptomic signatures can be integrated in multiple biological pathways, thus leading to a better comprehension of the fundamental mechanisms that are related to pathological processes influenced by smoking, and can be used for the selection of optimal therapy [16,17].
We anticipate that smoking affects molecular mechanisms, including transcriptomic patterns; therefore, we performed a TGCA data analysis in order to identify specific transcriptomic alterations in smoking versus non-smoking patients with head and neck cancer. These data provide a unique opportunity to study the potential oncogenic role of tobacco smoking; furthermore, based on the human papilloma virus (HPV)16 status and its effects on gene expression patterns, we intend to find whether these patients have different gene expression landscapes for the selected exposome components, or for the co-occurrence of these two important components of the exposome.

2. Materials and Methods

2.1. The Cancer Genome Atlas Gene Expression Data for Head and Neck Squamous Cell Carcinomas

The data were downloaded from the University of California Santa Cruz (USCS) Genome Browser [18] as a gene expression data matrix containing log2 transformed, normalized gene expression data for 519 tumor and 43 normal tumor-adjacent tissue samples (from 13 females and 30 males) and generated by RNA sequencing. Table 1 contains the demographic and clinical characteristics of the 519 patients that the tumor tissues were collected from, including clinical stages, information related to TNM (T: tumor, N: node, Mx: metastases), smoking history, and tumor localization.
The initial differential expression analysis was performed on the entire group of samples, namely 519 HNSCC tumors and 43 tumor-adjacent normal tissues. Further analyses were conducted on groups of patients divided according to their smoking status at the moment of sample collection: current smokers, quitters (also known as reformed smokers or ex-smokers), and never smokers.
Differential expression analysis was performed using the Gene Spring GX v.13.0 software from Agilent Technologies (Santa Clara, CA, USA), using the “volcano plot” module, and applying a fold change cut-off of ±2, moderated t-test and false discovery rate (FDR) correction. The bioinformatics analyses for differential expression were performed in the case of tumor tissue (n = 519) versus normal tissue (n = 43), and different comparisons were made based on smoking status, such as current smoking (n = 174) versus never smoked (n = 118), having quit smoking (n = 209) versus never smoked (n = 118), and finally for currently smoking (n = 174) versus having quit smoking (n = 209).

2.2. Molecular Classification for Gene Expression

Signature was performed using different online tools, such as String version 10.5 [19], Kyoto Encyclopedia of Genes and Genomes (KEEG) pathways [20], PantherDB [21] and miRnet data base [22].

2.3. Survival Analysis

Kaplan–Meier survival analysis was performed to investigate the survival distribution between selected groups based on the smoking status using Graph Pad Prism software (Version 6, Graph Pad software Inc., San Diego, CA, USA). A plot of the Kaplan–Meier analysis with the selected groups based on smoking status was performed.

3. Results

3.1. Differential Gene Expression in Tumor Tissues Versus Normal Tissues for Head and Neck Squamous Cell Carcinomas

Global gene expression was evaluated in tumor tissues (n = 519) versus normal tissues (n = 43), where we identified 1216 upregulated genes and 1751 downregulated genes considering as cut-off the fold change (FC) value of ±2 and p-value ≤0.001 (Benjamini–Hochberg correction). Based on the KEGG classification, most of the upregulated genes belong to the extracellular matrix ECM–receptor interaction, focal adhesion, the PI3K (Phosphoinositide 3-kinase) –Akt (Protein kinase B) signaling pathway, or cell cycle regulation, while downregulated genes are involved in altered pathways belonging to drug metabolism cytochrome P450, chemical carcinogenesis and metabolism of xenobiotic by cytochrome P450. We also generated a map describing interconnections of the altered genes with specific targeting miRNAs (microRNAs) in Figure 1, using the miRnet database.

3.2. Differential Gene Expression Levels in Smokers Compared with Non-Smokers or Ex-Smokers in Head and Neck Squamous Cell Carcinomas

In order to evaluate the transcriptomic alterations related to smoking in the tumor tissues, we performed multiple comparisons. In the first analysis we compared the gene expression pattern between current smokers and patients who had never smoked, identifying 119 altered transcripts (9 downregulated genes and 110 upregulated) (Table 2). The second analysis compared the gene expression pattern between ex-smokers and patients who had never smoked, revealing 24 altered transcripts (22 upregulated genes and 2 downregulated). The String network is presented in Figure 2A, for the 20 common upregulated genes for the group of ex-smokers vs. never smoked, and current smokers vs. never smoked, respectively. The third analysis compared the gene expression pattern between current smokers versus those who had quit, revealing 15 overexpressed genes, with the String network displayed in Figure 2A; the genes are not connected in specific networks. Figure 2B shows downregulated genes, presented as Venn diagrams for the analyzed groups (current smokers, ex-smokers, non-smokers), and emphasizes a signature in the ex-smokers group, revealing a panel of genes with an altered expression level event after quitting smoking in HNSCC patients.
In Figure S1, the String Network was generated for the altered signature in the case of smoking versus non-smoking, showing that most genes are involved in the metabolism of xenobiotics by cytochrome P450.

3.3. Molecular Classification for Altered Gene Expression Signature in Smoking versus Never Smoking Head and Neck Squamous Cell Carcinomas Patients

In order to perform the classification of the 119 altered genes in smoking versus never smoking HNSCC patients we used different online tools, such as String database [19] KEGG pathways [20] and PantherDB [21].
The String network for the modified gene expression is presented in Figure 3A, and the KEGG classification in Figure 3B. In the KEGG pathways classification, most of the genes are related to cellular metabolism, emphasizing the activation of detoxification pathways, chemical carcinogenesis, or drug metabolism. Gene ontology classification based on molecular function and biological processes is presented in Table 3.

3.4. Effect of Smoking on Head and Neck Squamous Cell Carcinomas Stages

To evaluate potential alteration of gene expression, specific for early stages to advanced tumor status, we compared the molecular profiles of tumor tissue in smoking versus non-smoking, according to tumor stages. Results from the computation of specific gene expressions for current smokers versus never smokers on HNSCC stages identified stage-specific gene expression signatures. The data presented as Venn diagrams illustrate the common and different overexpressed and downregulated genes (Figure 4A,B), while the network created using the String tool for the downregulated genes in the case of current smokers versus never smokers on HNSCC stage 1, revealed 32 genes involved in cell cycle regulation (Figure 4C).
To address the probable alteration of gene expression as an effect of advanced tumor status, we restricted the gene ontology classification of smoking versus non-smoking only for stage 1 tumors. Gene ontology based on the altered signature of smoking versus never smoking for stage 1 HNSCC is illustrated in Table 4.

3.5. Evaluation of Gene Expression Signature Based on Human Papilloma Virus 16 Status with/without Correlation with Smoking in Head and Neck Squamous Cell Carcinomas Patients

Using the same TCGA data, we performed a new analysis which allowed the identification of a signature composed of 2087 genes (1143 downregulated and 844 upregulated genes) that discriminates HPV16-induced HNSCC from their HPV-negative counterparts, comprising a patient cohort of 37 patients HPV-positive (HPV+) for subtype 16 and 72 HPV-negative (HPV−) patients. Using the miRnet data base, an analysis of the altered transcripts revealed the most relevant interconnected miRNAs and the most significantly altered pathways (Figure 5).
A differential expression level comparison was performed, taking into consideration as reference group the nonsmoker patients negative for HPV (HPV− Smoke−), represented by 32 cases, while the other three analyzed groups were represented by smoking patients that were HPV-positive (HPV+ Smoke+, 11 cases), nonsmoking patients that were HPV-positive (HPV+ Smoke−), and smoking patients that were HPV-negative (Smoke+ HPV−, 11 patients). The heat map depicted in Figure 6A illustrates a specific signature in each analyzed group and in Table 5 being presented GO classification for the altered expression signature identified based on HPV-16 status.
Regarding the overexpressed genes, we observed a common signature represented by 724 genes in the case of the HPV+ versus HPV−, and the group (HPV+ Smoke−) versus (Smoke− HPV−); based on the KEGG classification, these genes are related to the ECM–receptor interaction, focal adhesion, and PI3K–Akt signaling (Figure 6B). We also identified 374 genes specifics for HPV+ versus HPV−, and 309 genes specific for the HPV+ Smoke− group (Figure 6B). In addition, 507 common genes were identified for the overexpressed genes involved in DNA replication and cell cycle, as obtained by KEGG classification (Figure 6C).

3.6. Survival Prognosis Analysis Related to Smoking Status in Head and Neck Squamous Cell Carcinomas Patients

The overall survival of HNSC patients related to three different groups based on smoking status: current smoker (n = 174), ex-smoker (n = 209), and never-smoking groups (n = 118) are presented in Figure S2. Also, survival analysis was performed in the case of HPV16+ (n = 72) versus HPV16– group (n = 37), observing a slightly increased survival rate in HPV− patients compared to HPV+ cases.

4. Discussion

The HNSCC disease etiology is complex, being related to genetic background and exposome, where smoking and viral infection are two important players in its causality [3,18,23,24,25,26,27]. HPV and smoking converge in more aggressive diseases through complex altered pathways (particularly those related to xenobiotic metabolism [23,28,29]) as observed in the presented data, with potentially important clinical implications. At the same time, smoking patients have a reduced overall survival when compared to non-smoking groups [30,31]; in our case, we can observe a slightly increased survival rate in the non-smoking group, with no statistical significance. The study of Osazuwa–Peters et al. shows that the survival rate is almost double in the non-smoking versus smoking group with HNSCC [32].
The overall variation in gene expression profiles for patients who quit smoking versus those who never smoked, and current smokers versus those who quit, was different when comparing tumors with normal tumor adjacent samples. The most significant differences were observed in the case of smoking versus never smoking. These observations are sustained by similar studies [33,34,35]. A set of 49 differentially expressed genes were detected based on smoking status, targeting NFkB-related pathways [36]. In comprehensive genomic characterization, we showed that most of the altered genes are related to the regulation of mutated TP53 and cell cycle progression [37]. A study similar to ours emphasized the important role of xenobiotic metabolizing enzymes in several cancer types like bladder [38] and lung cancer [39] or leukemia [40]; cytochrome (CY) P450 enzymes such as CYP1A1 are activated in the case of the smoking group as compared to never smokers [41]. Xenobiotic metabolizing enzymes CYP1A1 and CYP1B1 were observed also in a cellular model of oral leukoplakia [42]. The metabolic detoxification pathways have an important role in chemotherapeutics metabolism [42], affecting the response to therapy in smoking groups. The negative effects can be counteracted by chemopreventive agents [43,44,45,46].
Our study demonstrated the complex biologic effects of smoking through the analysis related to the effect of smoking on HNSCC stages, particularly for stage 1, emphasizing the altered pathways leading to carcinogenesis. In the case of gene expression signatures in smoking versus never smoking for stage 1 HNSCC patients, using the PantherDB online tool we observed an important number of representative transcripts that are responsible for biological adhesion, including for early stages (DLL3, CDH17, TINAGL1, STRC, PCDHAC2, PCDHB13, TNR, PCDHGB6, PCDHGA9, PLXNB3). The same analysis identified 32 downregulated genes related to cell cycle regulation. A cell culture-based study on human placental cells using cigarette smoke extracts showed alterations in cell cycle, cell migration, and endocrine activity [44,47], sustaining our findings. Alterations of these vital genes denote a frequent mechanism essential for the susceptibility to a variety of smoking-induced diseases [48]. These adhesion molecules are retrieved in the circulatory system and not only at tumor sites, especially in advanced stages [49,50,51]. These adhesion effectors and angiogenic markers could thus be used as biomarkers of invasion and metastasis [49,50,52,53,54,55]. Immune and inflammation-related genes may provide a better understanding of the mechanisms through which tobacco smoking causes disease [56], as well as the possible benefits of immune agonist therapy. A previous study showed that tumors with a genetic smoking pattern had decreased immune infiltration, connected with an unfavorable survival rate [57]. It was shown that the circulating immune markers of inflammation could mirror the overall immune and inflammatory cancer-promoting microenvironment [58,59], and may suggest probable etiologic mechanisms related to smoking-induced diseases, particularly in HNSCC. In our case, the altered immune and inflammatory response genes for stage I smoking vs. nonsmoking identified a panel of 13 genes obtained by PantherDB classification on biological processes (CCL26, TRAF5, XCL1, BLK, DLL3, MAPK8IP2, CD1B, CD1E, CRYAB, GPX3, CCL22, LHX4, ULBP1). Smoking also affects immunity in the oral cavity and promotes oral cavity diseases, including oral cancer [59]; hence, there is no doubt that immune microenvironment of HNSCC significantly affects the response to therapy [60].
The global gene expression signatures show the interaction between genetics and exposure characteristics, making the subtraction of a single agent effect related to HNSCC very difficult; in spite of this we were able to demonstrate that the transcriptomic pattern is highly influenced by tobacco smoking and HPV status. In the case of HPV+ HNSCCs, the mutational and transcriptomic pattern appeared similar to that of cervical cancer, with a higher mutation incidence in the PI3K pathway and DNA repair genes [61]. For transcriptomic pattern, altered genes implicated in cell cycle, apoptosis, inflammatory response, DNA replication and repair, or other important transcription factors involved in transcription regulation were revealed [47]. Gene expression characterization in HPV+ tumors can be used to predict response to therapy, and this information could be used for better tailored therapies [62]. Our data sustain the idea that HPV-related HNSCC represents a distinct entity, and that current treatment options are not responding to the needs of these patients [47,61], proving the necessity for routine testing of HPV in clinical practice [63] and at the same time underlining the importance of patient stratification based on smoking status or HPV infection.

5. Conclusions

In conclusion, we demonstrated that smoking and HPV infection make important contributions to the HNSCC genomic portrait; this information can be translated into the creation of better-tailored therapies. A part of the gene expression alteration pattern was the reversible signature related mainly to the metabolism of xenobiotics by cytochrome P450, but some genes remained altered even after quitting smoking. This study emphasizes the utility of HNSCC classification based on smoking status in the management of cancer risk and also in establishing therapeutic options based on the many altered cell signaling pathways that we identified (metabolic detoxification pathways, adhesion cell signaling, or immune and inflammation pathways).
Our pathway analysis was able to identify the most relevant gene expression signature for the smoking HNSCC patients related to xenobiotic metabolizing enzymes, which might affect the response to therapy. These data support the fact that smoking is a major risk factor for HNSCC outcomes and that smoking cessation therapy should be a part of standard HNSCC care. From a research viewpoint, these results emphasize the importance of environmental toxic agent exposure in corroboration with genetic background.
Since a wide range of factors affect gene expression, it is very likely that not all relevant gene transcripts were identified, and some genes with altered expression levels may not have been confirmed; however, these data represent an important starting point for new investigations. Future studies could examine other complexities of the transcriptome in relation with other environmental carcinogens. Regarding the exposome, is difficult to analyze single exposures due to the fact that for most cases a co-occurrence of toxic elements is observed.

Supplementary Materials

Author Contributions

A.I.I. wrote the manuscript; C.B., R.C. and N.M. analyzed the data; L.M., A.O., C.C., D.D. and S.B. performed the bibliographic search and participated in writing the discussion section. I.B.-N. designed the study and supervised the final correction of the manuscript. All the authors aided in the preparation of the manuscript and endorsed the final version of the manuscript.

Funding

This work was supported by the follow research grants: PN-III-P4-ID-PCE-2016-0795:“Addressing the complex exposome profile in hormone-dependent cancers of the breast and prostate and its influence on tumoral genome”, “Genomic mapping of population from polluted area with radioactivity and heavy metals to increase national security” and Project No. 164/2017.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abrahao, R.; Anantharaman, D.; Gaborieau, V.; Abedi-Ardekani, B.; Lagiou, P.; Lagiou, A.; Ahrens, W.; Holcatova, I.; Betka, J.; Merletti, F.; et al. The influence of smoking, age and stage at diagnosis on the survival after larynx, hypopharynx and oral cavity cancers in Europe: The ARCAGE study. Int. J. Cancer 2018, 143, 32–44. [Google Scholar] [CrossRef] [PubMed]
  2. Irimie, A.I.; Braicu, C.; Cojocneanu-Petric, R.; Berindan-Neagoe, I.; Campian, R.S. Novel technologies for oral squamous carcinoma biomarkers in diagnostics and prognostics. Acta Odontol. Scand. 2015, 73, 161–168. [Google Scholar] [CrossRef] [PubMed]
  3. Giraldi, L.; Leoncini, E.; Pastorino, R.; Wunsch-Filho, V.; de Carvalho, M.; Lopez, R.; Cadoni, G.; Arzani, D.; Petrelli, L.; Matsuo, K.; et al. Alcohol and cigarette consumption predict mortality in patients with head and neck cancer: A pooled analysis within the International Head and Neck Cancer Epidemiology (INHANCE) Consortium. Ann. Oncol. 2017, 28, 2843–2851. [Google Scholar] [CrossRef] [PubMed]
  4. Irimie, A.I.; Zimta, A.A.; Ciocan, C.; Mehterov, N.; Dudea, D.; Braicu, C.; Berindan-Neagoe, I. The Unforeseen Non-Coding RNAs in Head and Neck Cancer. Genes 2018, 9, 134. [Google Scholar] [CrossRef] [PubMed]
  5. Ferlay, J.; Shin, H.R.; Bray, F.; Forman, D.; Mathers, C.; Parkin, D.M. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int. J. Cancer 2010, 127, 2893–2917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Salyakina, D.; Tsinoremas, N.F. Non-coding RNAs profiling in head and neck cancers. NPJ Genom. Med. 2016, 1, 15004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Irimie, A.I.; Braicu, C.; Sonea, L.; Zimta, A.A.; Cojocneanu-Petric, R.; Tonchev, K.; Mehterov, N.; Diudea, D.; Buduru, S.; Berindan-Neagoe, I. A looking-glass of non-coding RNAs in oral cancer. Int. J. Mol. Sci. 2017, 18, 2620. [Google Scholar] [CrossRef] [PubMed]
  8. Irimie, A.I.; Sonea, L.; Jurj, A.; Mehterov, N.; Zimta, A.A.; Budisan, L.; Braicu, C.; Berindan-Neagoe, I. Future trends and emerging issues for nanodelivery systems in oral and oropharyngeal cancer. Int. J. Nanomed. 2017, 12, 4593–4606. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Feng, H.-M.; Zhao, Y.; Zhang, J.-P.; Zhang, J.-H.; Jiang, P.; Li, B.; Wang, C. Expression and potential mechanism of metabolism-related genes and CRLS1 in non-small cell lung cancer. Oncol. Lett. 2018, 15, 2661–2668. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, T.H.; Hsia, S.M.; Shih, Y.H.; Shieh, T.M. Association of smoking, alcohol use, and betel quid chewing with epigenetic aberrations in cancers. Int. J. Mol. Sci. 2017, 18, 1210. [Google Scholar] [CrossRef] [PubMed]
  11. Leemans, C.R.; Snijders, P.J.F.; Brakenhoff, R.H. The molecular landscape of head and neck cancer. Nat. Rev. Cancer 2018, 18, 269–282. [Google Scholar] [CrossRef] [PubMed]
  12. Vineis, P.; Chadeau-Hyam, M.; Gmuender, H.; Gulliver, J.; Herceg, Z.; Kleinjans, J.; Kogevinas, M.; Kyrtopoulos, S.; Nieuwenhuijsen, M.; Phillips, D.H.; et al. The exposome in practice: Design of the EXPOsOMICS project. Int. J. Hyg. Environ. Health 2017, 220, 142–151. [Google Scholar] [CrossRef] [PubMed]
  13. Van Breda, S.G.J.; Wilms, L.C.; Gaj, S.; Jennen, D.G.J.; Briedé, J.J.; Kleinjans, J.C.S.; de Kok, T.M.C.M. The exposome concept in a human nutrigenomics study: Evaluating the impact of exposure to a complex mixture of phytochemicals using transcriptomics signatures. Mutagenesis 2015, 30, 723–731. [Google Scholar] [CrossRef] [PubMed]
  14. Tang, Z.; Li, C.; Kang, B.; Gao, G.; Li, C.; Zhang, Z. GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017, 45, W98–W102. [Google Scholar] [CrossRef] [PubMed]
  15. Singh, R.K.; Sivabalakrishnan, M. Feature selection of gene expression data for cancer classification: A review. Procedia Comput. Sci. 2015, 50, 52–57. [Google Scholar] [CrossRef]
  16. Zhang, X.; Cha, I.-H.; Kim, K.-Y. Highly preserved consensus gene modules in human papilloma virus 16 positive cervical cancer and head and neck cancers. Oncotarget 2017, 8, 114031–114040. [Google Scholar] [CrossRef] [PubMed]
  17. Braicu, C.; Catana, C.; Calin, G.A.; Berindan-Neagoe, I. NCRNA combined therapy as future treatment option for cancer. Curr. Pharm. Des. 2014, 20, 6565–6574. [Google Scholar] [CrossRef] [PubMed]
  18. UCSC Genome Browser. Available online: https://genome.ucsc.edu (accessed on 12 March 2018).
  19. STRING: Functional Protein Association Network, Version 10.5. Available online: https://string-db.org (accessed on 16 May 2018).
  20. KEGG Pathway Maps. Available online: http://www.genome.jp/kegg/pathway.html (accessed on 10 June 2018).
  21. Panther Clasiffication System. Available online: http://www.pantherdb.org (accessed on 14 June 2018).
  22. miRNet-Network-Based Visual Analysis of miRNAs, Targets and Functions. Available online: http://www.mirnet.ca/faces/home.xhtml (accessed on 24 June 2018).
  23. Campbell, J.D.; Yau, C.; Bowlby, R.; Liu, Y.; Brennan, K.; Fan, H.; Taylor, A.M.; Wang, C.; Walter, V.; Akbani, R.; et al. Genomic, pathway network, and immunologic features distinguishing squamous carcinomas. Cell Rep. 2018, 23, 194–212.e6. [Google Scholar] [CrossRef] [PubMed]
  24. Ferketich, A.K.; Niland, J.C.; Mamet, R.; Zornosa, C.; D’Amico, T.A.; Ettinger, D.S.; Kalemkerian, G.P.; Pisters, K.M.; Reid, M.E.; Otterson, G.A. Smoking status and survival in the national comprehensive cancer network non-small cell lung cancer cohort. Cancer 2013, 119, 847–853. [Google Scholar] [CrossRef] [PubMed]
  25. Bryant, A.; Cerfolio, R.J. Differences in epidemiology, histology, and survival between cigarette smokers and never-smokers who develop non-small cell lung cancer. Chest 2007, 132, 185–192. [Google Scholar] [CrossRef] [PubMed]
  26. Koshiaris, C.; Aveyard, P.; Oke, J.; Ryan, R.; Szatkowski, L.; Stevens, R.; Farley, A. Smoking cessation and survival in lung, upper aero-digestive tract and bladder cancer: Cohort study. Br. J. Cancer 2017, 117, 1224–1232. [Google Scholar] [CrossRef] [PubMed]
  27. Pierce, J.P.; Patterson, R.E.; Senger, C.M.; Flatt, S.W.; Caan, B.J.; Natarajan, L.; Nechuta, S.J.; Poole, E.M.; Shu, X.-O.; Chen, W.Y. Lifetime cigarette smoking and breast cancer prognosis in the after breast cancer pooling project. JNCI J. Nat. Cancer Inst. 2014, 106, djt359. [Google Scholar] [CrossRef] [PubMed]
  28. Champion, M.; Brennan, K.; Croonenborghs, T.; Gentles, A.J.; Pochet, N.; Gevaert, O. Module analysis captures pancancer genetically and epigenetically deregulated cancer driver genes for smoking and antiviral response. EBioMedicine 2018, 27, 156–166. [Google Scholar] [CrossRef] [PubMed]
  29. Yavorski, J.M.; Blanck, G. Smoking correlates with increased cytoskeletal protein-related coding region mutations in the lung and head and neck datasets of the cancer genome atlas. Physiol. Rep. 2016, 4, e13045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Hoffmann, M.; Quabius, E.S.; Tribius, S.; Hebebrand, L.; Gorogh, T.; Halec, G.; Kahn, T.; Hedderich, J.; Rocken, C.; Haag, J.; et al. Human papillomavirus infection in head and neck cancer: The role of the secretory leukocyte protease inhibitor. Oncol. Rep. 2013, 29, 1962–1968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Quabius, E.S.; Moller, P.; Haag, J.; Pfannenschmidt, S.; Hedderich, J.; Gorogh, T.; Rocken, C.; Hoffmann, M. The role of the antileukoprotease SLPI in smoking-induced human papillomavirus-independent head and neck squamous cell carcinomas. Int. J. Cancer 2014, 134, 1323–1334. [Google Scholar] [CrossRef] [PubMed]
  32. Osazuwa-Peters, N.; Adjei Boakye, E.; Chen, B.Y.; Tobo, B.B.; Varvares, M.A. Association Between Head and Neck Squamous Cell Carcinoma Survival, Smoking at Diagnosis, and Marital Status. JAMA Otolaryngol. 2018, 144, 42–50. [Google Scholar] [CrossRef] [PubMed]
  33. Rudin, C.M.; Avila-Tang, E.; Harris, C.C.; Herman, J.G.; Hirsch, F.R.; Pao, W.; Schwartz, A.G.; Vahakangas, K.H.; Samet, J.M. Lung cancer in never smokers: Molecular profiles and therapeutic implications. Clin. Cancer Res. 2009, 15, 5646–5661. [Google Scholar] [CrossRef] [PubMed]
  34. Delagranda, A.; Leterme, G.; Chirpaz, E.; Ferdynus, C.; Fernandez, C.; Rubin, F. Epidemiological features of cancers of the oral cavity, oropharynx, hypopharynx and larynx cancer in Reunion Island. Eur. Ann. Otorhinolaryngol. Head Neck Dis. 2018, 135, 175–181. [Google Scholar] [CrossRef] [PubMed]
  35. Cufari, M.E.; Proli, C.; De Sousa, P.; Raubenheimer, H.; Al Sahaf, M.; Chavan, H.; Shedden, L.; Niwaz, Z.; Leung, M.; Nicholson, A.G.; et al. Increasing frequency of non-smoking lung cancer: Presentation of patients with early disease to a tertiary institution in the UK. Eur. J. Cancer 2017, 84, 55–59. [Google Scholar] [CrossRef] [PubMed]
  36. Farshadpour, F.; Roepman, P.; Hordijk, G.J.; Koole, R.; Slootweg, P.J. A gene expression profile for non-smoking and non-drinking patients with head and neck cancer. Oral Dis. 2012, 18, 178–183. [Google Scholar] [CrossRef] [PubMed]
  37. Cancer Genome Atlas Network. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [Green Version]
  38. Rouissi, K.; Ouerhani, S.; Hamrita, B.; Bougatef, K.; Marrakchi, R.; Cherif, M.; Ben Slama, M.R.; Bouzouita, M.; Chebil, M.; Ben Ammar Elgaaied, A. Smoking and polymorphisms in xenobiotic metabolism and DNA repair genes are additive risk factors affecting bladder cancer in Northern Tunisia. Pathol. Oncol. Res. 2011, 17, 879–886. [Google Scholar] [CrossRef] [PubMed]
  39. Ihsan, R.; Chauhan, P.S.; Mishra, A.K.; Yadav, D.S.; Kaushal, M.; Sharma, J.D.; Zomawia, E.; Verma, Y.; Kapur, S.; Saxena, S. Multiple analytical approaches reveal distinct gene-environment interactions in smokers and non smokers in lung cancer. PLoS ONE 2011, 6, e29431. [Google Scholar] [CrossRef] [PubMed]
  40. Chauhan, P.S.; Ihsan, R.; Mishra, A.K.; Yadav, D.S.; Saluja, S.; Mittal, V.; Saxena, S.; Kapur, S. High order interactions of xenobiotic metabolizing genes and P53 codon 72 polymorphisms in acute leukemia. Environ. Mol. Mutagen. 2012, 53, 619–630. [Google Scholar] [CrossRef] [PubMed]
  41. Petros, W.P.; Younis, I.R.; Ford, J.N.; Weed, S.A. Effects of tobacco smoking & nicotine on cancer treatment. Pharmacotherapy 2012, 32, 920–931. [Google Scholar] [PubMed]
  42. Gümüş, Z.H.; Du, B.; Kacker, A.; Boyle, J.O.; Bocker, J.M.; Mukherjee, P.; Subbaramaiah, K.; Dannenberg, A.J.; Weinstein, H. Effects of tobacco smoke on gene expression and cellular pathways in a cellular model of oral leukoplakia. Cancer Prev. Res. 2008, 1, 100–111. [Google Scholar] [CrossRef] [PubMed]
  43. Irimie, A.I.; Braicu, C.; Zanoaga, O.; Pileczki, V.; Gherman, C.; Berindan-Neagoe, I.; Campian, R.S. Epigallocatechin-3-gallate suppresses cell proliferation and promotes apoptosis and autophagy in oral cancer SSC-4 cells. Onco Targets Ther. 2015, 8, 461–470. [Google Scholar] [PubMed]
  44. Braicu, C.; Mehterov, N.; Vladimirov, B.; Sarafian, V.; Nabavi, S.M.; Atanasov, A.G.; Berindan-Neagoe, I. Nutrigenomics in cancer: Revisiting the effects of natural compounds. Semin. Cancer Biol. 2017, 46, 84–106. [Google Scholar] [CrossRef] [PubMed]
  45. Budisan, L.; Gulei, D.; Zanoaga, O.M.; Irimie, A.I.; Sergiu, C.; Braicu, C.; Gherman, C.D.; Berindan-Neagoe, I. Dietary intervention by phytochemicals and their role in modulating coding and non-coding genes in cancer. Int. J. Mol. Sci. 2017, 18, 1178. [Google Scholar] [CrossRef] [PubMed]
  46. Kim, C.W.; Lee, H.M.; Lee, K.; Kim, B.; Lee, M.Y.; Choi, K.C. Effects of cigarette smoke extracts on cell cycle, cell migration and endocrine activity in human placental cells. Reprod. Toxicol. 2017, 73, 8–19. [Google Scholar] [CrossRef] [PubMed]
  47. Mirghani, H.; Ugolin, N.; Ory, C.; Lefevre, M.; Baulande, S.; Hofman, P.; St Guily, J.L.; Chevillard, S.; Lacave, R. A predictive transcriptomic signature of oropharyngeal cancer according to HPV16 status exclusively. Oral Oncol. 2014, 50, 1025–1034. [Google Scholar] [CrossRef] [PubMed]
  48. Scott, D.A.; Palmer, R.M. The influence of tobacco smoking on adhesion molecule profiles. Tob. Induc. Dis. 2003, 1, 7–25. [Google Scholar] [CrossRef] [PubMed]
  49. Zahra, A.; Rubab, I.; Malik, S.; Khan, A.; Khan, M.J.; Fatmi, M.Q. Meta-Analysis of miRNAs and their involvement as biomarkers in oral cancers. BioMed Res. Int. 2018, 2018, 8439820. [Google Scholar] [CrossRef] [PubMed]
  50. Palmer, R.M.; Stapleton, J.A.; Sutherland, G.; Coward, P.Y.; Wilson, R.F.; Scott, D.A. Effect of nicotine replacement and quitting smoking on circulating adhesion molecule profiles (sICAM-1, sCD44v5, sCD44v6). Eur. J. Clin. Investig. 2002, 32, 852–857. [Google Scholar] [CrossRef]
  51. Scott, D.A.; Stapleton, J.A.; Palmer, R.M.; Wilson, R.F.; Sutherland, G.; Coward, P.Y.; Gustavsson, G.; Odell, E.W.; Poston, R.N. Plasma concentrations of reputed tumor-associated soluble CD44 isoforms (v5 and v6) in smokers are dose related and decline on smoking cessation. Cancer Epidemiol. Biomark. Prev. 2000, 9, 1211–1214. [Google Scholar]
  52. Ionescu, C.; Braicu, C.; Chiorean, R.; Cojocneanu Petric, R.; Neagoe, E.; Pop, L.; Chira, S.; Berindan-Neagoe, I. TIMP-1 expression in human colorectal cancer is associated with SMAD3 gene expression levels: A pilot study. JGLD 2014, 23, 413–418. [Google Scholar] [PubMed]
  53. Braicu, C.; Tudoran, O.; Balacescu, L.; Catana, C.; Neagoe, E.; Berindan-Neagoe, I.; Ionescu, C. The significance of PDGF expression in serum of colorectal carcinoma patients—Correlation with Duke’s classification. Can PDGF become a potential biomarker? Chirurgia 2013, 108, 849–854. [Google Scholar] [PubMed]
  54. Ionescu, S.; Barbu, E.; Ionescu, C.; Costache, A.; Balasoiu, M. Giant gastrointestinal stromal tumor of the stomach. Rom. J. Morphol. Embryol. 2015, 56, 239–245. [Google Scholar] [PubMed]
  55. Scott, D.A.; Todd, D.H.; Coward, P.Y.; Wilson, R.F.; Odell, E.W.; Poston, R.N.; Matthews, J.P.; Palmer, R.M. The acute influence of tobacco smoking on adhesion molecule expression on monocytes and neutrophils and on circulating adhesion molecule levels in vivo. Addict. Biol. 2000, 5, 195–205. [Google Scholar] [CrossRef] [PubMed]
  56. Shiels, M.S.; Katki, H.A.; Freedman, N.D.; Purdue, M.P.; Wentzensen, N.; Trabert, B.; Kitahara, C.M.; Furr, M.; Li, Y.; Kemp, T.J.; et al. Cigarette smoking and variations in systemic immune and inflammation markers. J. Natl. Cancer Inst. 2014, 106, dju294. [Google Scholar] [CrossRef] [PubMed]
  57. Mandal, R.; Şenbabaoğlu, Y.; Desrichard, A.; Havel, J.J.; Dalin, M.G.; Riaz, N.; Lee, K.-W.; Ganly, I.; Hakimi, A.A.; Chan, T.A.; et al. The head and neck cancer immune landscape and its immunotherapeutic implications. JCI Insight 2016, 1, e89829. [Google Scholar] [CrossRef] [PubMed]
  58. Desrichard, A.; Kuo, F.; Chowell, D.; Lee, K.W.; Riaz, N.; Wong, R.J.; Chan, T.A.; Morris, L.G.T. Tobacco Smoking-Associated Alterations in the Immune Microenvironment of Squamous Cell Carcinomas. J. Natl. Cancer Inst. 2018. [Google Scholar] [CrossRef] [PubMed]
  59. Lee, J.; Taneja, V.; Vassallo, R. Cigarette smoking and inflammation: Cellular and molecular mechanisms. J. Dent. Res. 2012, 91, 142–149. [Google Scholar] [CrossRef] [PubMed]
  60. Solomon, B.; Young, R.J.; Rischin, D. Head and neck squamous cell carcinoma: Genomics and emerging biomarkers for immunomodulatory cancer treatments. Semin. Cancer Biol. 2018, in press. [Google Scholar] [CrossRef] [PubMed]
  61. Dok, R.; Nuyts, S. HPV positive head and neck cancers: Molecular pathogenesis and evolving treatment strategies. Cancers 2016, 8, 41. [Google Scholar] [CrossRef] [PubMed]
  62. Thibodeau, B.J.; Geddes, T.J.; Fortier, L.E.; Ahmed, S.; Pruetz, B.L.; Wobb, J.; Chen, P.; Wilson, G.D.; Akervall, J.A. Gene expression Characterization of HPV positive head and neck cancer to predict response to chemoradiation. Head Neck Pathol. 2015, 9, 345–353. [Google Scholar] [CrossRef] [PubMed]
  63. Bishop, J.A.; Lewis, J.S., Jr.; Rocco, J.W.; Faquin, W.C. HPV-related squamous cell carcinoma of the head and neck: An update on testing in routine pathology practice. Semin. Diagn. Pathol. 2015, 32, 344–351. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The interconnected genes with specific miRNAs using miRnet [22] involved in focal adhesion, extracellular matrix (ECM)–receptor interaction, or pathways in cancer interconnected with targeting microRNAs (miRNAs).
Figure 1. The interconnected genes with specific miRNAs using miRnet [22] involved in focal adhesion, extracellular matrix (ECM)–receptor interaction, or pathways in cancer interconnected with targeting microRNAs (miRNAs).
Ijerph 15 01558 g001
Figure 2. Venn diagram used for overlapping the altered gene expression pattern in the case of the three studied groups. Common and specific gene expression signatures for the three groups of HNSCC (head and neck squamous cell carcinomas) patients: non-smokers, ex-smokers, and smokers. (A) For the case of overexpressed genes, 20 common genes from non-smokers vs. smokers and ex-smokers vs. smokers, integrated as network using String, version 10.5) [19]; (B) The case of downregulated genes; (C) Heat maps for the expression level for the three HNSCC patient groups (current smokers, quitters, non-smokers), in dark blue being presented the downregulated genes and in red those overexpressed genes, generated using Gene Spring version 13.0.
Figure 2. Venn diagram used for overlapping the altered gene expression pattern in the case of the three studied groups. Common and specific gene expression signatures for the three groups of HNSCC (head and neck squamous cell carcinomas) patients: non-smokers, ex-smokers, and smokers. (A) For the case of overexpressed genes, 20 common genes from non-smokers vs. smokers and ex-smokers vs. smokers, integrated as network using String, version 10.5) [19]; (B) The case of downregulated genes; (C) Heat maps for the expression level for the three HNSCC patient groups (current smokers, quitters, non-smokers), in dark blue being presented the downregulated genes and in red those overexpressed genes, generated using Gene Spring version 13.0.
Ijerph 15 01558 g002
Figure 3. Gene network generated using Sting program; (A) Network of the interconnected genes; (B) KEGG (Kyoto Encyclopedia of Genes and Genomes) classification based on the altered signature in smoking versus never-smoking HNSCC patients.
Figure 3. Gene network generated using Sting program; (A) Network of the interconnected genes; (B) KEGG (Kyoto Encyclopedia of Genes and Genomes) classification based on the altered signature in smoking versus never-smoking HNSCC patients.
Ijerph 15 01558 g003
Figure 4. Venn diagram overlapping the (A) overexpressed and (B) downregulated genes based on a specific analysis of current smokers versus never smokers with HNSCC stage 1, 2, 3, and 4; (C) String Network for the case of downregulated genes in current smokers versus never smokers in HNSCC stage 1, where the genes involved in cell cycle regulation are highlighted in red.
Figure 4. Venn diagram overlapping the (A) overexpressed and (B) downregulated genes based on a specific analysis of current smokers versus never smokers with HNSCC stage 1, 2, 3, and 4; (C) String Network for the case of downregulated genes in current smokers versus never smokers in HNSCC stage 1, where the genes involved in cell cycle regulation are highlighted in red.
Ijerph 15 01558 g004
Figure 5. The interconnected genes with specific miRNAs using miRnet [22] involved in focal adhesion, ECM–receptor interaction or gap junction.
Figure 5. The interconnected genes with specific miRNAs using miRnet [22] involved in focal adhesion, ECM–receptor interaction or gap junction.
Ijerph 15 01558 g005
Figure 6. Gene expression signature based on HPV16 status with/without correlation with smoking in HNSCC patients. (A) Heat maps representing the expression level in the HNSCC patient group based on smoking and HPV status. For nonsmoking and HPV16-negative patients (Smoking− HPV−, n = 32), we had Smoking+ HPV− (n = 11), Smoking+ HPV+ (n = 11), Smoking− HPV+ (n = 11), in dark blue being presented the downregulated genes and in red those overexpressed genes, generated using Gene Spring version 13.0. (B) Venn diagram showing the differential signature in the case of the overexpressed genes highlighting the main altered pathways as displayed by KEGG classification. (C) Venn diagram to emphasize that the differential signature in the case of the overexpressed genes underlines the main altered pathways as obtained from String Network and KEGG (Kyoto Encyclopedia of Genes and Genomes) classification, with red dots showing the genes involved in cytokine–cytokine receptor interaction and blue dots the cell adhesion molecules.
Figure 6. Gene expression signature based on HPV16 status with/without correlation with smoking in HNSCC patients. (A) Heat maps representing the expression level in the HNSCC patient group based on smoking and HPV status. For nonsmoking and HPV16-negative patients (Smoking− HPV−, n = 32), we had Smoking+ HPV− (n = 11), Smoking+ HPV+ (n = 11), Smoking− HPV+ (n = 11), in dark blue being presented the downregulated genes and in red those overexpressed genes, generated using Gene Spring version 13.0. (B) Venn diagram showing the differential signature in the case of the overexpressed genes highlighting the main altered pathways as displayed by KEGG classification. (C) Venn diagram to emphasize that the differential signature in the case of the overexpressed genes underlines the main altered pathways as obtained from String Network and KEGG (Kyoto Encyclopedia of Genes and Genomes) classification, with red dots showing the genes involved in cytokine–cytokine receptor interaction and blue dots the cell adhesion molecules.
Ijerph 15 01558 g006
Table 1. Succinct presentation of the TCGA characteristics for patients diagnosed with HNSCCs used for gene expression analysis. Head and neck squamous cell carcinomas: HNSCCs, T: tumor, N: node, M0 no metastasis, Mx: metastases presence; HPV: human papilloma virus; TCGA: The Cancer Genome Atlas.
Table 1. Succinct presentation of the TCGA characteristics for patients diagnosed with HNSCCs used for gene expression analysis. Head and neck squamous cell carcinomas: HNSCCs, T: tumor, N: node, M0 no metastasis, Mx: metastases presence; HPV: human papilloma virus; TCGA: The Cancer Genome Atlas.
Clinical ParametersPatients (n = 519)Females/Males
SexFemales134
Males374
Undeclared11
AgeMedian, range61, 19–90
Median, range males59, 19–88
Median, range females64.5, 24–90
Undeclared11
Clinical stage1209/11
29833/65
310126/75
427564/211
Unknown25
Clinical TNMT1N0M0187/11
T1N1M052/3
T1N2M071/6
T1N2M11-/1
T1NxM022/-
T1N0M09633/63
T1N0Mx1-/1
T2N1M0123/9
T2N1Mx11/-
T2N2M0336/27
T2N2M11-/1
T2N2Mx3-/3
T1N3M01-/1
T2NxM01-/1
T3N0M05815/43
T3N0Mx11/-
T2N1M0204/16
T3N1M111/-
T3N2M0468/38
T3N2M111/-
T3N2Mx1-/1
T3N3M01-/1
T3NxM02-/2
T4N0M06617/49
T4N0Mx11/-
T4N1M03913/26
T4N2M06214/48
T4N2Mx1-/1
T4N3M072/5
T4NxM03-/3
TxN1M01-/1
TxN2M01-/1
TxNxMx9-/9
Unknown16
Smoking historySmoker17435/139
Reformed smoker <15 years13426/108
Reformed smoker >15 years7219/53
Reformed, unknown years20/2
Lifelong non-smoker11449/65
Unknown23
Anatomic neoplasm subdivisionAlveolar ridge18
Base of tongue27
Buccal mucosa20
Floor of mouth60
Hard palate7
Hypopharynx9
Larynx114
Lip3
Oral cavity71
Oral tongue127
Oropharynx10
Tonsil42
Unknown11
HPVPositive72
Negative37
Unknown410
Table 2. Genes with altered expression levels, based on a fold change (FC) ± 2, p-value ≤ 0.05 for smoking versus never smoking HNSCC patients.
Table 2. Genes with altered expression levels, based on a fold change (FC) ± 2, p-value ≤ 0.05 for smoking versus never smoking HNSCC patients.
GeneFC (abs)p-ValueRegulationGeneFC (abs)p-ValueRegulation
CDKN2A−2.597370.000603DownEPHA72.394610.000111Up
RYR3−2.359241.45 × 10−6DownpsiTPTE222.3943926.22 × 10−6Up
KRT2−2.317610.002097DownPOU6F22.3932423.2 × 10−6Up
KCNS1−2.314950.000321DownSOHLH12.3906840.000256Up
MYBPC1−2.100910.046549DownLTF2.3808630.012372Up
ARL14−2.082830.000138DownMLXIPL2.3738078 × 10−6Up
IL13RA2−2.059068.82 × 10−5DownGLI22.3722546.52 × 10−7Up
FLRT3−2.023150.002341DownNLGN4Y2.3588720.00092Up
PLA2G2F−2.013070.002746DownPAK72.3556396.22 × 10−6Up
NTS5.2246548 × 10−6UpFIBCD12.3450150.000181Up
RPS4Y14.9240497.27 × 10−5UpGATA42.3433890.000117Up
UGT1A64.5057073.3 × 10−7UpPANX22.3195081.31 × 10−6Up
UPK1B4.3368694.15 × 10−5UpPCYT1B2.3146193.29 × 10−5Up
MGST14.2966364.68 × 10−7UpFGF192.311253.18 × 10−5Up
CYP1A14.2699585.3 × 10−13UpSLC44A42.292020.000356Up
C20orf1143.9856580.001123UpSCN2A2.289850.00011Up
SCGB3A13.9854513.93 × 10−5UpPROM12.28420.001556Up
GPR153.9136922.69 × 10−15UpCYorf15A2.2577710.002582Up
DDX3Y3.7484380.000214UpTFPI22.251930.000222Up
MUC5B3.5736620.001207UpMSI12.2493052.06 × 10−5Up
CNNM13.5089773.02 × 10−8UpADD22.2377416.45 × 10−5Up
CES13.415441.21 × 10−5UpALDH1A12.2348980.000377Up
PRAME3.2890140.000277UpERN22.2109540.000992Up
GPX23.2521813.95 × 10−6UpLGI32.2061280.000183Up
CYP26A13.2412532.07 × 10−5UpPRKY2.2061140.002395Up
NR5A13.1862911.24 × 10−5UpSALL12.2004588.32 × 10−6Up
PPP1R1B3.180471.38 × 10−5UpTBX52.1948918.95 × 10−6Up
FGFBP23.1562942.95 × 10−6UpHOXA42.1816794.15 × 10−6Up
SLC13A53.1425126.04 × 10−6UpTNNI32.181054.06 × 10−5Up
DMBT13.1028930.000721UpPLUNC2.1797970.009085Up
KRTCAP33.0772787.92 × 10−12UpNGB2.1772740.000323Up
BPIL13.0118470.002582UpFOLR12.1751270.000402Up
UGT1A83.008191.43 × 10−5UpGDA2.1632710.00378Up
PTH2R2.8495312.55 × 10−5UpAKR1C32.1509828.66 × 10−5Up
KRT192.8125640.001704UpAZGP12.1423510.004176Up
GAL2.7929546.04 × 10−6UpCCNA12.1385280.001642Up
EIF1AY2.7618050.001147UpPCDH192.1310330.000887Up
WNK22.7563873.74 × 10−5UpGJB72.1284820.000382Up
B4GALNT42.7313525.73 × 10−7UpWDR722.1199210.003053Up
RAB3B2.7216292.09 × 10−7UpCLDN82.1187580.001164Up
FAM132A2.671932.95 × 10−6UpCBS2.1179952.03 × 10−5Up
HOXA72.6664793.29 × 10−5UpMSMB2.1172110.003146Up
PIGR2.6523870.004096UpCFTR2.1128520.00034Up
BCHE2.6390823.2 × 10−6UpNTRK22.1123790.001737Up
UGT82.6353065.71 × 10−5UpFGF132.1081372.2 × 10−5Up
USP9Y2.6295080.001642UpRPL39L2.0851145.8 × 10−6Up
PDIA22.6137512.12 × 10−6UpSLC29A42.0794584.68 × 10−7Up
ZFY2.5581550.0013UpADH72.074140.012696Up
CALB12.5452230.002413UpPIWIL22.0689650.000399Up
AKR1C12.5407837.2 × 10−6UpCYP1B12.0577248.81 × 10−5Up
UTY2.5339130.00193UpCPNE72.0553571.57 × 10−6Up
ATP13A52.5306330.00011UpBRDT2.0412220.001179Up
SLC5A122.5040152.95 × 10−6UpCHGA2.0333922.12 × 10−5Up
FOXJ12.4912340.000516UpABO2.0326630.001283Up
PRDM132.4796014.03 × 10−7UpSTATH2.0234240.020511Up
HORMAD12.4373060.001219UpSCN9A2.0181960.00048Up
UCHL12.4344167.96 × 10−6UpADAMTS202.0111340.000245Up
NPW2.414236.04 × 10−6UpRBM112.0106224.12 × 10−5Up
PNCK2.3962838.83 × 10−5UpZNF5562.0096571.24 × 10−5Up
TCF152.0074994.31 × 10−5Up
Table 3. Gene ontology (GO) classification based on the gene expression signature in smoking versus never-smoking HNSCC patients using the PantherDB online tool [21].
Table 3. Gene ontology (GO) classification based on the gene expression signature in smoking versus never-smoking HNSCC patients using the PantherDB online tool [21].
Ontology FunctionTypeNo. MoleculesPercent (%)
Molecular functionbinding (GO:0005488)3336.7%
catalytic activity (GO:0003824)3033.3%
transporter activity (GO:0005215)1415.6%
receptor activity (GO:0004872)44.4%
signal transducer activity (GO:0004871)44.4%
structural molecule activity (GO:0005198)33.3%
translation regulator activity (GO:0045182)11.1%
antioxidant activity (GO:0016209)11.1%
Biological processcellular process (GO:0009987)4828.1%
metabolic process (GO:0008152)3218.7%
biological regulation (GO:0065007)2313.5%
developmental process (GO:0032502)169.4%
response to stimulus (GO:0050896)158.8%
multicellular organismal process (GO:0032501)137.6%
localization (GO:0051179)84.7%
cellular component organization or biogenesis (GO:0071840)74.1%
biological adhesion (GO:0022610)42.3%
locomotion (GO:0040011)31.8%
immune system process (GO:0002376)10.6%
reproduction (GO:0000003)10.6%
Protein classtransporter (PC00227)1013.5%
hydrolase (PC00121)912.2%
oxidoreductase (PC00176)810.8%
transcription factor (PC00218)810.8%
nucleic acid binding (PC00171)79.5%
signaling molecule (PC00207)68.1%
transferase (PC00220)56.8%
enzyme modulator (PC00095)56.8%
receptor (PC00197)34.1%
extracellular matrix protein (PC00102)22.7%
cytoskeletal protein (PC00085)22.7%
transfer/carrier protein (PC00219)22.7%
cell junction protein (PC00070)22.7%
lyase (PC00144)11.4%
calcium-binding protein (PC00060)11.4%
defense/immunity protein (PC00090)11.4%
membrane traffic protein (PC00150)11.4%
isomerase (PC00135)11.4%
Table 4. Gene ontology (GO) classification based on the gene expression signature in smoking versus never smoking for stage 1 HNSCC patients using the PantherDB online tool, showing the maximum effect of smoking and the minimum effect of tumor environment.
Table 4. Gene ontology (GO) classification based on the gene expression signature in smoking versus never smoking for stage 1 HNSCC patients using the PantherDB online tool, showing the maximum effect of smoking and the minimum effect of tumor environment.
Ontology FunctionTypeNo. MoleculesPercent (%)
Molecular functionbinding (GO:0005488)143337.5%
catalytic activity (GO:0003824)126233.0%
transporter activity (GO:0005215)42911.2%
receptor activity (GO:0004872)2556.7%
signal transducer activity (GO:0004871)2316.0%
structural molecule activity (GO:0005198)1574.1%
antioxidant activity (GO:0016209)230.6%
translation regulator activity (GO:0045182)170.4%
channel regulator activity (GO:0016247)130.3%
Biological processcellular process (GO:0009987)247929.5%
metabolic process (GO:0008152)158918.9%
biological regulation (GO:0065007)91110.9%
response to stimulus (GO:0050896)6968.3%
localization (GO:0051179)5836.9%
cellular component organization or biogenesis (GO:0071840)5816.9%
developmental process (GO:0032502)5346.4%
multicellular organismal process (GO:0032501)5236.2%
immune system process (GO:0002376)1461.7%
locomotion (GO:0040011)1201.4%
biological adhesion (GO:0022610)1101.3%
reproduction (GO:0000003)831.0%
rhythmic process (GO:0048511)320.4%
cell killing (GO:0001906)60.1%
Protein classnucleic acid binding (PC00171)56515.2%
transcription factor (PC00218)50413.5%
hydrolase (PC00121)41011.0%
receptor (PC00197)3559.5%
transporter (PC00227)2697.2%
signaling molecule (PC00207)2516.7%
transferase (PC00220)2406.4%
cytoskeletal protein (PC00085)1925.2%
enzyme modulator (PC00095)1794.8%
oxidoreductase (PC00176)1494.0%
extracellular matrix protein (PC00102)842.3%
membrane traffic protein (PC00150)832.2%
ligase (PC00142)792.1%
calcium-binding protein (PC00060)721.9%
structural protein (PC00211)631.7%
isomerase (PC00135)421.1%
lyase (PC00144)361.0%
defense/immunity protein (PC00090)350.9%
cell adhesion molecule (PC00069)310.8%
cell junction protein (PC00070)310.8%
chaperone (PC00072)240.6%
transfer/carrier protein (PC00219)190.5%
transmembrane receptor regulatory/adaptor protein (PC00226)110.3%
Table 5. Gene ontology (GO) classification based on the gene expression signature in HPV16-positive versus HPV16-negative patients using the PantherDB online tool [21].
Table 5. Gene ontology (GO) classification based on the gene expression signature in HPV16-positive versus HPV16-negative patients using the PantherDB online tool [21].
Ontology FunctionTypeNo. MoleculesPercent (%)
Molecular functionbinding (GO:0005488)690041.1%
catalytic activity (GO:0003824)515630.7%
transporter activity (GO:0005215)14748.8%
receptor activity (GO:0004872)14508.6%
signal transducer activity (GO:0004871)8955.3%
structural molecule activity (GO:0005198)8264.9%
translation regulator activity (GO:0045182)530.3%
antioxidant activity (GO:0016209)230.1%
binding (GO:0005488)690041.1%
Biological processcellular process (GO:0009987)1096828.1%
metabolic process (GO:0008152)658916.9%
biological regulation (GO:0065007)406510.4%
response to stimulus (GO:0050896)34678.9%
developmental process (GO:0032502)33198.5%
multicellular organismal process (GO:0032501)29547.6%
cellular component organization or biogenesis (GO:0071840)26266.7%
localization (GO:0051179)21255.4%
immune system process (GO:0002376)9752.5%
biological adhesion (GO:0022610)8982.3%
locomotion (GO:0040011)6681.7%
reproduction (GO:0000003)2890.7%
rhythmic process (GO:0048511)280.1%
growth (GO:0040007)270.1%
Protein classhydrolase (PC00121)223513.2%
nucleic acid binding (PC00171)16769.9%
signaling molecule (PC00207)15989.5%
transcription factor (PC00218)15519.2%
enzyme modulator (PC00095)14628.7%
receptor (PC00197)13247.8%
cytoskeletal protein (PC00085)9875.9%
transferase (PC00220)9875.9%
transporter (PC00227)9775.8%
oxidoreductase (PC00176)6674.0%
extracellular matrix protein (PC00102)6323.7%
cell adhesion molecule (PC00069)5603.3%
calcium-binding protein (PC00060)3622.1%
membrane traffic protein (PC00150)3312.0%
cell junction protein (PC00070)3071.8%
defense/immunity protein (PC00090)2781.6%
ligase (PC00142)1851.1%
structural protein (PC00211)1731.0%
chaperone (PC00072)1570.9%
transmembrane receptor regulatory/adaptor protein (PC00226)1470.9%
lyase (PC00144)1340.8%
isomerase (PC00135)660.4%
transfer/carrier protein (PC00219)630.4%

Share and Cite

MDPI and ACS Style

Irimie, A.I.; Braicu, C.; Cojocneanu, R.; Magdo, L.; Onaciu, A.; Ciocan, C.; Mehterov, N.; Dudea, D.; Buduru, S.; Berindan-Neagoe, I. Differential Effect of Smoking on Gene Expression in Head and Neck Cancer Patients. Int. J. Environ. Res. Public Health 2018, 15, 1558. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15071558

AMA Style

Irimie AI, Braicu C, Cojocneanu R, Magdo L, Onaciu A, Ciocan C, Mehterov N, Dudea D, Buduru S, Berindan-Neagoe I. Differential Effect of Smoking on Gene Expression in Head and Neck Cancer Patients. International Journal of Environmental Research and Public Health. 2018; 15(7):1558. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15071558

Chicago/Turabian Style

Irimie, Alexandra Iulia, Cornelia Braicu, Roxana Cojocneanu, Lorand Magdo, Anca Onaciu, Cristina Ciocan, Nikolay Mehterov, Diana Dudea, Smaranda Buduru, and Ioana Berindan-Neagoe. 2018. "Differential Effect of Smoking on Gene Expression in Head and Neck Cancer Patients" International Journal of Environmental Research and Public Health 15, no. 7: 1558. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15071558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop