- Research article
- Open Access
The association between HPV gene expression, inflammatory agents and cellular genes involved in EMT in lung cancer tissue
BMC Cancer volume 20, Article number: 916 (2020)
Lung cancer is a leading cause of cancer morbidity and mortality worldwide. Several studies have suggested that Human papillomavirus (HPV) infection is an important risk factor in the development of lung cancer. In this study, we aim to address the role of HPV in the development of lung cancer mechanistically by examining the induction of inflammation and epithelial-mesenchymal transition (EMT) by this virus.
In this case-control study, tissue samples were collected from 102 cases with lung cancer and 48 controls. We examined the presence of HPV DNA and also the viral genotype in positive samples. We also examined the expression of viral genes (E2, E6 and E7), anti-carcinogenic genes (p53, retinoblastoma (RB)), and inflammatory cytokines in HPV positive cases.
HPV DNA was detected in 52.9% (54/102) of the case samples and in 25% (12/48) of controls. A significant association was observed between a HPV positive status and lung cancer (OR = 3.37, 95% C.I = 1.58–7.22, P = 0.001). The most prevalent virus genotype in the patients was type 16 (38.8%). The expression of p53 and RB were decreased while and inflammatory cytokines were increased in HPV-positive lung cancer and HPV-positive control tissues compared to HPV-negative lung cancer and HPV-negative control tissues. Also, the expression level of E-cad and PTPN-13 genes were decreased in HPV- positive samples while the expression level of SLUG, TWIST and N-cad was increased in HPV-positive samples compared to negative samples.
Our study suggests that HPV infection drives the induction of inflammation and EMT which may promote in the development of lung cancer.
Lung cancer is one of the leading causes of cancer morbidity and mortality worldwide . There are several types of primary lung cancer which, are divided into two main groups; small cell lung cancer (SCLC) and non-small cell lung cancer (NSCLC). NSCLC are divided into three common types; squamous cell carcinoma, large cell carcinoma and adenocarcinoma . The pathogenesis of lung cancer is a complex multifactor process with both genetic and environmental factors playing a major role . Infectious agents are emerging as key drivers in the development of cancer [4,5,6,7]. Previously, numerous infectious agents have been shown to be involved in a myriad of lung diseases including cancer, Idiopathic Pulmonary Fibrosis (IPF) and Chronic Obstructive Pulmonary Disease (COPD) [8,9,10].
Human papilloma virus (HPV) is one of the most important human oncogenic viruses , which has previously been shown to be associated with numerous cancers including lung, breast and prostate [1, 6, 11,12,13]. The HPV genome is divided into three main sections; long control region (LCR), early region (E) encoding E1, E2, E4–E7, and late region (L) consisting of L1 and L2 . E6 and E7 are the oncoproteins that act as stimulating factors for host cell proliferation . E6 interacts with p53 and BCL2, while E7 interacts with retinoblastoma (RB); both of which lead to enhanced cell proliferation, resistance to apoptosis and chromosomal instability [16, 17]. These viral proteins enhance tumour development by promoting inflammation and epithelial-mesenchymal transition (EMT) [18, 19].
In response to harmful stimuli and invading pathogens, the innate immune system becomes activated through a variety of receptors, leading to the generation of an acute inflammatory response. This inflammation aids in the removal and clearance of the stimulus. However, should the stimulus fail to be removed the development of chronic inflammation occurs which is strongly associated with cancer .
Chronic inflammation as a result of viral infection is responsible for an estimated 25% of all human cancers [21, 22]. In response to viral infection the generation of a pro-inflammatory response involves activation of numerous transcription factors including NF-κB and the secretion of numerous pro-inflammatory cytokines and metabolites including transforming growth factors like beta (TGF-β), interleukin 1 (IL-1), IL-6, IL-11, Tumour necrosis factor α (TNF-α) and reactive oxygen-nitrogen species (RONS) - all of which play a pro-tumorigenic role in the context of chronic inflammation. This pro-inflammatory tissue microenvironment results in the suppression of anti-humoral immunity and also the promotion of tumour development and metastasis [7, 23, 24].
The second facet of high-risk HPV (hr-HPV) related tumour development is EMT, which plays an important role in solid cancer progression through multiple biochemical changes. EMT is well known to enhance cell migration, invasion and cancer development .
There are several genes involved in EMT, including SLUG, PTPN13, E-cad, N-cad and TWIST. SLUG protein is involved in important cellular events including EMT and also has anti-apoptotic activity . PTPN13 interacts with Fas receptor which is indirectly involved in inhibition of programmed cell death . E-cad and N-cad expression levels have also been connected with survival mechanisms and metastasis of lung cancer cells [28, 29].
In this study we investigated, for the first time, the role of hr-HPV in EMT and lung tumour development. We also assessed the prevalence of HPV in lung tumour samples; examining the expression level of viral and cellular genes and the associations between these expressed genes in EMT and lung tumour development.
Study design and samples
This case-control study was conducted between November 2017 and September 2018. One hundred and two lung cancer samples and forty-eight normal lung tissue samples. Control samples were age and sex matched, with the tissue samples collected from a peripheral region of the surgically removed lung cancers and non-cancer patients with fibrosis. All samples, cases and controls, were fresh tissue with a Tumor Proportion Score (TPS) > 50%. Control samples were age and sex matched. The TNM system was used to denote the stage of cancer as decided by a consultant oncologist and oncological surgeon. Gender, age, smoking status, tumour type and tumour stage were clinical parameters of patients that are shown in (supplementary materials). We had no medical records of HPV infection before cancer diagnosis.
Extraction of nucleic acids
Total DNA extraction from tissue samples was performed by QIAamp® DNA Mini Kit (Qiagen, Hilden, Germany). Quality of extracted DNA was assessed by conducting PCR for β-globin as described before . All samples were deemed suitable for molecular analysis due to β-globin gene amplification.
Total RNA extraction was conducted by RNeasy Mini Kit (Qiagen, Hilden, Germany).
HPV detection and genotyping
HPV genome detection was conducted using PCR for L1 and E7 genes . HPV genotyping was performed by INNO-LiPA HPV Genotyping v2 test (Innogenetics, Ghent, Belgium).
Determination of HPV genome physical status
To determine the physical status of the HPV genome, the E2/E6 ratio was used. An E2/E6 ratio > 0 and < 1 indicates that the virus is in a mixed physical state, with both episomal and integrated forms of the virus .
Quantitative real-time PCR
mRNA level detection of viral genes
Total RNA was extracted and purified from the tissue by using RNEasy Mini kit (QIAGEN, Hilden, Germany). For cDNA synthesis, 1 μg of total RNA was reverse transcribed using the QuantiNova Reverse Transcription Kit (QIAGEN, Germany). CDNA synthesis was performed in a thermal cycler in the following order: 27 °C for 10 min, 38 °C for 15 min, 44 °C for 40 min, 72 °C for 15 min. All the primers which were used to detect viral genes (E2, E6 and E7) are listed in a table in the (supplementary materials). To detect viral genes E2, E6 and E7, Quantitative SYBR green TaqMan Universal PCR Master Mix® (QIAGEN, Germany), one step RT-PCR® kits (QIAGEN, Hilden, Germany) and QuantiNova Reverse Transcription® Kit were used, respectively.
For viral genes we used serial dilutions of E2, E6 and E7 genes cloned in PUC57 vector (GenScript, Jiangsu, China). Serial dilution was containing equivalent amounts of these genes from 72 to 865 million copies per reaction, served as a standard control.
mRNA level detection of cellular genes
cDNA was synthesized using the PrimeScript First Strand cDNA synthesis kit (TaKaRa Bio, Kusatsu, Japan). Quantitative RT-PCR analyses were performed using the Power SYBR Green PCR Master Mix (TaKaRa Bio, Kusatsu, Japan). The relative expression level of each mRNA was normalized using GAPDH. The primers are listed in supplementary materials.
Enzyme linked immunosorbent assay (ELISA)
For tissue homogenization, all fresh tissue samples were weighed and the tissue lysate was prepared according to the manufactures protocol (Invitrogen, CA, USA). Approximately 50 μg of each tissue was excised and washed with ice-cold PBS.
The level of p53, RB, IL-1, IL-6, IL-11, NF-kB, NF-κ PTPN13, E-cadherin, N-cadherin and TWIST was assessed using Abcam’s p53 Simple Step ELISA® Kit (Abcam, Cambridge, MA, USA), Human Retinoblastoma ELISA® kit (Sigma-Aldrich, Saint Louis, USA), Human Retinoblastoma ELISA® kit (Sigma-Aldrich, Saint Louis, USA), Human IL-6 ELISA® Kit, Human IL-1 beta ELISA® Kit, Human IL-11 ELISA® Kit, NF-kB p65 Transcription Factor Assay® Kit (Abcam, Cambridge, MA, USA), Human Tyrosine-Protein Phosphatase Non-Receptor Type 13 (PTPN13) ELISA Kit (MyBiosource, USA), Human E-Cadherin, N-Cadherin ELISA Kit (Abcam, Cambridge, MA, USA).and TWIST ELISA Kit (Aviva Systems Biology, CA, USA).
Quantification of RONS
The RONS level was assessed by OxiSelect™ Intracellular ROS/RNS Assay kit (Cell Biolabs, Inc., San Diego, CA). For this purpose, cell lysate was used and preparation of this based on Kit instructions.
Continuous variables are presented as mean ± standard deviation and categorical variables are presented as N (%). Normality test was checked using Kolmogorov–Smirnov test for the continuous variables. For comparing the central tendency (e.g. mean for normal and median for non-normal variables) between two groups, two-independent samples t-test or Mann-Whitney non-parametric test and between more than two groups, one-way ANOVA or kruskal-wallis test were used. Chi-square/ or Fisher exact test was performed for assessing the associations of the categorical variables. The unit of all expression RT-PCR is (2^-DCt)*1000. Internal normalization was performed using an internal housekeeping or reference gene (GAPDH) and external normalization was applied by standardized approach. In addition, correlation analysis was done by Spearman’s correlation coefficient between viral and cellular factors. All of statistical analyses were analysed using GraphPad Prism 6 and STATA software versions 11.2. False discovery rate was corrected by Benjamini-Hochberg approach for multiple comparisons. A two-sided P-value of less than 0.05 was considered as statistical significance.
In this case-control study, we examined 102 lung cancer cases and 48 controls, with the mean ± SD age; 56.36 ± 12.49 and 57.0 ± 12.24, respectively. Seventy-four (72.5%) of the cases and 31 (64.5%) of the controls were male, respectively. The cases and control groups were matched based on age (p = 0.77). There were three types of lung tumour tissues; squamous-cell carcinoma (51.9%), adenocarcinoma (32.3%) and SCLC (15.7%). The highest and lowest stages of cancer in this study were IIIB (30.4%) and IA and IIB (1.9%) respectively. HPV DNA was detected in 52.9% of the lung cancer specimens and in 25% of control samples. There was a significant association between the presence of HPV and lung tumour (OR = 3.37, 95% C.I = 1.58–7.22, P = 0.001). Genotype 16 was the most frequently isolated genotype in both cases (38.8%) and controls (50%). No significant association was observed between all genotypes and the occurrence of lung tumour (p = 0.651) (supplementary materials). HPV DNA was detected in 55.6% (30 of 53) of squamous-cell carcinoma samples, 54.5% (18 of 33) of adenocarcinoma samples and 37.5% (6 of 16) of SCLC samples. The association between HPV infection and histopathological types of tumour was not statistically significant (p = 0.434). There were no significant differences in the frequency distributions of lung tumour stages between HPV+ and HPV- groups (p = 0.163). More information is presented in supplementary materials.
In the HPV+ lung carcinoma patients, the virus was present in its integrated form in 27.8% of cases. The incidences of episomal and mixed forms of HPV genome were 5.5 and 66.7% respectively. In the control HPV+ group, the incidence of HPV genome status was 25, 0 and 75% integrated, episomal and mixed forms of HPV respectively (Table 1). The gene expression level of viral genes in both types and stages of lung tumour are shown in Table 2. The highest level of viral gene expression was that of E7 which was most highly observed in stage IV samples (mean ± SD:13.56 ± 5.13). The lowest level of viral gene expression examined was E6 in stage IB samples (mean ± SD: 3.0 ± 1.75). The gene expression level of viral factors E2 and E6 were highest in stage IIB and stage IV respectively. Stratification of the samples based on the tumour type reveals the expression level of E7 in adenocarcinoma samples (mean ± SD: 11.94 ± 4.93) and E2 in SCLC (mean ± SD: 3.67 ± 1.15) were the highest and lowest respectively (Table 2). The expression level of viral genes in control samples and tumour samples are illustrated in Fig. 1.
In Table 3, the level of cellular factors such as tumour-suppressors (Rb and p53), inflammatory factors (ILs, IFNs, TGF-β, TNF-α, and NF-κB), EMT factors (PTPN13, SLUG, E-cad, N-cad and TWIST) and RONS are presented. The protein levels of Rb and p53 were significantly downregulated in HPV+ cases and HPV+ controls compared with HPV- cases and controls (p < 0.001). The level of inflammatory factors, were considerably higher in HPV+ cases and controls compared to the HPV- cases and controls groups. The levels of EMT involved factors found to be significantly higher in HPV infected group compare to HPV non-infected group (p < 0.001 for all). Among the EMT involved genes, PTPN13 and E-cad were significantly downregulated in HPV+ cases and controls compared with HPV- cases and controls (p < 0.001). SLUG, N-cad and TWIST were significantly upregulated in HPV+ cases and controls compared with HPV- cases and controls (p < 0.05). The highest expression levels were related with SLUG, N-cad and TWIST in HPV+ compared with HPV- groups (fold change > 15; p < 0.001 for all). More details are presented in Fig. 2. Significant negative correlations were observed between the expression level of viral genes and the protein expression levels of regulatory host proteins, Rb and p53. Among the inflammatory factors examined, the correlations between E2 expression level with IL-1 and TNF-α were statistically significant, and the correlations between IL-6 with E6 and E7 were statistically significant (p < 0.01). The correlation between E2 expression level and PTPN13 was positive but with SLUG, E-cad, N-cad and TWIST was negative. The expression level of E6 significantly correlated with the protein level of PTPN13. The expression level of E7 has the negative correlation with E-cad and N-cad (p < 0.05). Conversely, there were positive correlations between E6 gene expression and IL-1, IL-6, IFN-α and IFN-β protein levels and RONS production (p < 0.05) (Table 3).
Lung cancer is the primary cause of cancer death globally . As such, there is a major unmet clinical need for the development and discovery of prognostic biomarkers for the diagnosis of lung cancer. This need is underlined by the increased mortality rates which are currently being observed in lung cancer worldwide [1, 15]. A plethora of carcinogens are responsible for the initiation and development of various cancers. Of these, viral infections are implicated in approximately 18–20% of cancers [6, 11, 34]. While the prevalence of HPV in lung carcinoma has shown in numerous studies, to date, the role of hr-HPV in the promotion of EMT has not yet been clearly identified. Here, we report for the first time the association between HPV gene expression, inflammatory agents and cellular genes involved in EMT in lung cancer tissue.
In the current study, 52.9% of lung tumour samples were positive for HPV. Moreover, we demonstrate that increased expression of HPV genes is associated with decreased expression of regulatory cellular genes, RB and p53, and as a result increased risk of lung cancer. In an investigation Nadji et al. (2007, Iran) studied 141 lung carcinoma samples and 92 non-cancersamples as controls. Results demonstrated that 25.6% of cases and 9.0% of controls were positive for HPV infection respectively. The reported odds ratio for HPV infection was 3.48 (95% CI 1.522–7.958; P = 0.002) . A similar study by Argyria et al. (2017, Greece) investigated 67 lung tissue samples from a Greek cohort, 12 SCLC and 55 NSCLC, and detected HPV in 3.0% of the samples; with no association between HPV infection and lung cancer. Other studies conducted in this region showed the prevalence of HPV DNA rate was 0–61% . Robinson et al. (2016, USA) examined 70 NSCLC samples for detection of viral DNA. They detected 69% of HPV DNA in their samples . A meta-analysis study which performed by Syrjanen and his research team (2012) showed the prevalence of HPV in Europe (16.9%) in Australia (18.5%), in North America (12.5%) and in China and Taiwan (37.7%). This study shows the role of geographical distribution in HPV prevalence . Another point in the current study, was the highest SCC tumour type (51.9%) that is confirmed by a meta-analysis study (Almasi et al. 2016, Iran) . The most common HPV genotypes detected in cancer patients are HPV 16 and 18 [6, 38,39,40]. Previous investigations have been confirmed this issue and also, a higher prevalence of HPV-16 and 18 in Asian populations compared with European populations (lung samples). It should be noted that HPV-16 is the most prevalent genotype across all geographical areas . In another interesting study, the presence of HPV DNA was examined in the exhaled breath condensate (EBC) of lung cancer patients by Carpagnano et al. (2011, Italy). Their results showed the presence of HPV infection in 16.4% of samples. The authors state that analysis of EBC for HPV infection represents a valid tool for the diagnosis of airway colonisation .
Previous investigations have noted the physical status of HPV DNA as an important marker for tumour progression in other cancers, such as breast cancer . In this study, the highest integrated form was seen in stages III and IV in SCC samples (Table 1). Khodabandehlou et al. previously reported on the physical status of HPV genome in breast cancer samples, with 86% integrated and 14% mixed forms respectively. The largest number of integrated forms was in stage III and IV . Detection of HPV in its integrated form has also been reported in several other cancers [30, 42, 43]. The integration of HPV genome leads to changes in the expression of viral oncogenes (E6 and E7), dysregulating of critical cell cycle checkpoints, increased genetic instability in the host and finally tumour development .
We examined the potential role of HPV in lung cancer pathogenesis in two ways: i) the impact of HPV on the expression of genes involved in EMT, ii) the impact of HPV in the development of chronic inflammation and microenvironment alteration. EMT promotes cancer development through enhancing cellular migration and invasion [25, 45]. Oncoviruses are said to promote EMT in particular cancer cells, enabling the spread of metastatic cells from one location to another . Hr-HPV interacts with EMT factors to promote tumour development. Our results demonstrate that the levels of genes which promote EMT were substantially higher in HPV positive groups compare to HPV negative groups (p < 0.001 for all) (Table 2). This situation could indicate that the viral genes products/proteins may be involved in stimulating of transcription of these genes. We have hypothesized that hr-HPV promotes lung cancer development indirectly through a variety of different mechanisms. For example, HPV induce the production of ROS that leads to cell survival and resistance to programmed cell death . Previous studies have shown that in lung cancers with impaired E-cadherin expression, the frequency of lymph node metastases was significantly higher than tumours with high expression of the E-cadherin [28, 47, 48]. In our study, expression of E-cadherin in HPV+ samples were lower than HPV-negative samples (Fig. 2), with viral E7 detection having a negative correlation with E-cadherin levels (Table 2). Unlike E-cadherin, protein levels of N-cadherin in HPV+ samples were higher than HPV-negative samples; which has previously been shown to be associated with tumour development [49, 50]. On the other hand, TGF-β lead to an increase of N-cadherin and the expression of TGF-β in HPV+ samples was higher than HPV-negative samples (Fig. 2, Table 3) . Another important cellular factor is SLUG. This protein is overexpressed in numerous cancers . High expression levels of SLUG has also been shown to be associated with reduced E-cadherin expression, high histologic grade, lymph node metastasis, postoperative relapse and shorter survival in patients with cancer [51,52,53]. SLUG also has a role in inflammation-dependent tumour development . Our results demonstrate that the gene expression level of this factor was higher in HPV+ samples than HPV-negative samples (Fig. 2, Table 3). Furthermore, expression levels of SLUG have been shown to correlate with lung tumour development and drug resistance . Our results show the over expression of SLUG in HPV+ samples and direct correlation with E6 and E7 (Tables 2 and 4).
The tumour microenvironment is a key factor in tumour development and several epidemiologic and clinical researches have proposed a strong association between inflammation related to chronic infection and lung cancer [20, 56,57,58]. This inflammation affects different aspects of tumor development such as angiogenesis, survival of malignant cells and even tumor response to therapy [59, 60].
Our results demonstrate that the expression of numerous inflammatory factors was higher in HPV+ samples than HPV-negative samples (Table 3). Previously, Stone et al. (2014, Brazil) have shown HPV dependant changes in the tumour microenvironment. Their results showed differences in local inflammation between HPV+ and HPV-negative tumour tissues . Liu et al. (2015, China) have also studied the association between HPV and chronic inflammation, demonstrating that chronic inflammation was higher in oropharyngeal tumour tissue compared to normal tissues (P < 0.001). They propose that HPV infection could be considered as a biomarker/risk in some cancers in individuals with chronic inflammation . Previous investigations have shown microenvironmental alterations, caused by microorganisms, such as cytokine secretion promote epithelial proliferation. This issue was demonstrated in HPV infection and its persistence, which increases the risk of HPV transmission and oropharyngeal carcinogenesis [62,63,64].
In the current study the highest expression level of viral genes was in stage IV and the lowest level was in Stage IA and IB. In the other words, viral genes can be related to chronic inflammation and EMT (Table 2). This issue indicates the important role of these gene products in tumour development and metastasis. Although HPV is an oncovirus, the presence of the virus alone is insufficient for tumorigenesis. In order to promote cancer development, it is necessary to have a proinflammatory tumour microenvironment which occurs due to exposure to environmental factors or altered immunological mechanisms . Also, should be noted that the possibility to get HPV after premalignant lesions appear and how this concomitant infection may promote cancer progression but not lung cancer origin.
A key risk factor associated with HPV infection is smoking status. Previous studies have demonstrated the relationship between smoking and HPV infection in some cancers such as lung and cervical [1, 65]. In an investigation, relationship between HPV infection and cigarette smoking was studied (by Xi et al.). They demonstrated that HPV DNA load (type 16, 18) was associated with status of smoking, and current smokers had a higher HPV DNA load compared to former smokers .
Limitation of the current study were including:
i) Limitations on the number of cases considered for the study and the lack of statistical representation for some tumor stages; ii) protein and RNA samples are pooled representation of the different cell types from the original tumor/control tissue; iii) the absence of medical information regarding a HPV infection before cancer diagnosis for the patients analyzed in this research.
In summary, the presence of HPV was detected in 52.9% of lung cancer samples among which most were at stage III and IV (73.5%). Infection of HPV directly promotes local inflammation which in turn promotes tumorigenesis and cancer development. We demonstrate that HPV is associated with lung cancer development, although the role of hr-HPV in lung cancers requires further study. To the best of our knowledge, this is the first study reporting the role of HPV genes expression in EMT and the association between this virus and chronic inflammation in lung cancer patients.
Availability of data and materials
The datasets used and/or analyzed during the current study could become available through the corresponding author on reasonable request.
Squamous cell carcinoma
Small cell lung cancer
Non-small cell lung cancer
Idiopathic Pulmonary Fibrosis
Chronic Obstructive Pulmonary Disease
Long control region
Tumour necrosis factor α
Reactive oxygen-nitrogen species
Transforming growth factors like beta
Nuclear factor kappa beta
Hasegawa Y, Ando M, Kubo A, Isa S-i, Yamamoto S, Tsujino K, Kurata T, Ou S-H, Takada M, Kawaguchi T. Human papilloma virus in non-small cell lung cancer in never smokers: a systematic review of the literature. Lung Cancer. 2014;83(1):8–13.
Lu C, Onn A, Vaporciyan A. 78: Cancer of the lung, Holland-Frei Cancer medicine. People’s Medical Publishing House; 2010.
Nadji SA, Mokhtari-Azad T, Mahmoodi M, Yahyapour Y, Naghshvar F, Torabizadeh J, Ziaee AA, Nategh R. Relationship between lung cancer and human papillomavirus in north of Iran, Mazandaran province. Cancer Lett. 2007;248(1):41–6.
Etemadi A, Mostafaei S, Yari K, Ghasemi A, Minaei Chenar H, Moghoofei M. Detection and a possible link between parvovirus B19 and thyroid cancer. Tumor Biol. 2017;39(6):1010428317703634.
Goldszmid RS, Dzutsev A, Trinchieri G. Host immune response to infection and cancer: unexpected commonalities. Cell Host Microbe. 2014;15(3):295–305.
Khodabandehlou N, Mostafaei S, Etemadi A, Ghasemi A, Payandeh M, Hadifar S, Norooznezhad AH, Kazemnejad A, Moghoofei M. Human papilloma virus and breast cancer: the role of inflammation and viral expressed proteins. BMC Cancer. 2019;19(1):61.
Moghoofei M, Mostafaei S, Nesaei A, Etemadi A, Sadri Nahand J, Mirzaei H, Rashidi B, Babaei F, Khodabandehlou N. Epstein–Barr virus and thyroid cancer: the role of viral expressed proteins. J Cell Physiol. 2019;234(4):3790–9.
Jafarinejad H, Moghoofei M, Mostafaei S, Salimian J, Jamalkandi SA, Ahmadi A. Worldwide prevalence of viral infection in AECOPD patients: a meta-analysis. Microb Pathog. 2017;113:190–6.
Keyvani H, Moghoofei M, Bokharaei-Salim F, Mostafaei S, Mousavi S-AJ, Monavari SH, Esghaei M. Prevalence of respiratory viruses in Iranian patients with idiopathic pulmonary fibrosis. J Med Microbiol. 2017;66(11):1602–6.
Argyri E, Tsimplaki E, Marketos C, Politis G, Panotopoulou E. Investigating the role of human papillomavirus in lung cancer. Papillomavirus Res. 2017;3:7–10.
Krump NA, You J. Molecular mechanisms of viral oncogenesis in humans. Nat Rev Microbiol. 2018;1.
van Boerdonk RA, Daniels JM, Bloemena E, Krijgsman O, Steenbergen RD, Brakenhoff RH, Grünberg K, Ylstra B, Meijer CJ, Smit EF. High-risk human papillomavirus–positive lung cancer: molecular evidence for a pattern of pulmonary metastasis. J Thorac Oncol. 2013;8(6):711–8.
Moghoofei M, Keshavarz M, Ghorbani S, Babaei F, Nahand JS, Tavakoli A, Mortazavi HS, Marjani A, Mostafaei S, Monavari SH. Association between human papillomavirus infection and prostate cancer: A global systematic review and meta-analysis. Asia Pac J Clin Oncol. 0(0).
Morshed K, Polz-Gruszka D, Szymański M, Polz-Dacewicz M. Human papillomavirus (HPV)–structure, epidemiology and pathogenesis. Otolaryngol Pol. 2014;68(5):213–9.
Zur Hausen H. Papillomaviruses and cancer: from basic studies to clinical application. Nat Rev Cancer. 2002;2(5):342–50.
Jackson S, Harwood C, Thomas M, Banks L, Storey A. Role of Bak in UV-induced apoptosis in skin cancer and abrogation by HPV E6 proteins. Genes Dev. 2000;14(23):3065–73.
Yim E-K, Park J-S. The role of HPV E6 and E7 oncoproteins in HPV-associated cervical carcinogenesis. Cancer Res Treat. 2005;37(6):319.
Georgescu SR, Mitran CI, Mitran MI, Caruntu C, Sarbu MI, Matei C, Nicolae I, Tocut SM, Popa MI, Tampa M. New insights in the pathogenesis of HPV infection and the associated carcinogenic processes: the role of chronic inflammation and oxidative stress. J Immunol Res. 2018;2018:5315816.
DA COSTA RMG, BASTOS MMSM, MEDEIROS R, OLIVEIRA PA: The NFκB signaling pathway in papillomavirus-induced lesions: friend or foe? Anticancer Res 2016, 36(5):2073–2083.
Fernandes JV, TAAdM F, JCV DA, RNO C, MGF DC, Andrade VS, JMG DA. Link between chronic inflammation and human papillomavirus-induced carcinogenesis. Oncol Lett. 2015;9(3):1015–26.
Kakavandi E, Shahbahrami R, Goudarzi H, Eslami G, Faghihloo E. Anoikis resistance and oncoviruses. J Cell Biochem. 2017.
De Marzo AM, Platz EA, Sutcliffe S, Xu J, Grönberg H, Drake CG, Nakai Y, Isaacs WB, Nelson WG. Inflammation in prostate carcinogenesis. Nat Rev Cancer. 2007;7(4):256.
Xing M. Molecular pathogenesis and mechanisms of thyroid cancer. Nat Rev Cancer. 2013;13(3):184.
Wiseman H, Halliwell B. Damage to DNA by reactive oxygen and nitrogen species: role in inflammatory disease and progression to cancer. Biochem J. 1996;313(Pt 1):17.
Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial–mesenchymal transition. Nat Rev Mol Cell Biol. 2014;15(3):178–96.
Cohen ME, Yin M, Paznekas WA, Schertzer M, Wood S, Jabs EW. HumanSLUGGene organization, expression, and chromosome map location on 8q. Genomics. 1998;51(3):468–71.
Gross C, Heumann R, Erdmann KS. The protein kinase C-related kinase PRK2 interacts with the protein tyrosine phosphatase PTP-BL via a novel PDZ domain binding motif. FEBS Lett. 2001;496(2–3):101–4.
Lim SC, Jang IG, Kim YC, Park KO. The role of E-cadherin expression in non-small cell lung cancer. J Korean Med Sci. 2000;15(5):501–6.
Yamauchi M, Yoshino I, Yamaguchi R, Shimamura T, Nagasaki M, Imoto S, Niida A, Koizumi F, Kohno T, Yokota J. N-cadherin expression is a potential survival mechanism of gefitinib-resistant lung cancer cells. Am J Cancer Res. 2011;1(7):823.
Khan N, Castillo A, Koriyama C, Kijima Y, Umekita Y, Ohi Y, Higashi M, Sagara Y, Yoshinaka H, Tsuji T. Human papillomavirus detected in female breast carcinomas in Japan. Br J Cancer. 2008;99(3):408.
Kroupis C, Markou A, Vourlidis N, Dionyssiou-Asteriou A, Lianidou ES. Presence of high-risk human papillomavirus sequences in breast cancer tissues and association with histopathological characteristics. Clin Biochem. 2006;39(7):727–31.
Peitsaro P, Johansson B, Syrjänen S. Integrated human papillomavirus type 16 is frequently found in cervical cancer precursors as demonstrated by a novel quantitative real-time PCR technique. J Clin Microbiol. 2002;40(3):886–91.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Mostafaei S, Keshavarz M, Nahand JS, Hassankiadeh RF, Moradinazar M, Nouri M, Babaei F, Ahadi M, Payandeh M, Esker AS. Viral infections and risk of thyroid Cancer: a systematic review and empirical Bayesian meta-analysis. Pathology-Research and Practice. 2020;152855.
Robinson LA, Jaing CJ, Campbell CP, Magliocco A, Xiong Y, Magliocco G, Thissen JB, Antonia S. Molecular evidence of viral DNA in non-small cell lung cancer and non-neoplastic lung. Br J Cancer. 2016;115(4):497.
SYRJÄNEN K. Detection of human papillomavirus in lung cancer: systematic review and meta-analysis. Anticancer Res. 2012;32(8):3235–50.
Almasi Z, Salehiniya H, Amoori N, Enayatrad M. Epidemiology characteristics and trends of lung cancer incidence in Iran. Asian Pac J Cancer Prev. 2016;17(2):557–62.
Ragin C, Obikoya-Malomo M, Kim S, Chen Z, Flores-Obando R, Gibbs D, Koriyama C, Aguayo F, Koshiol J, Caporaso NE. HPV-associated lung cancers: an international pooled analysis. Carcinogenesis. 2014;35(6):1267–75.
Javanmard D, Moein M, Esghaei M, Naseripour M, Monavari SH, Bokharaei-Salim F, Sadeghipour A. Molecular evidence of human papillomaviruses in the retinoblastoma tumor. Virusdisease:1–7.
Medel-Flores O, Valenzuela-Rodríguez VA, Ocadiz-Delgado R, Castro-Muñoz LJ, Hernández-Leyva S, Lara-Hernández G, Silva-Escobedo J-G, Vidal PG, Sánchez-Monroy V. Association between HPV infection and prostate cancer in a Mexican population. Genet Mol Biol. 2018;41(4):781–9.
Carpagnano GE, Koutelou A, Me N, Martinelli D, Ruggieri C, Di Taranto A, Antonetti R, Carpagnano F, Foschino-Barbaro M. HPV in exhaled breath condensate of lung cancer patients. Br J Cancer. 2011;105(8):1183–90.
Pinatti L, Walline H, Carey T. Human papillomavirus genome integration and head and neck cancer. J Dent Res. 2018;97(6):691–700.
Gao W, Chen Y, Zhang Y, Zhang Q, Zhang L. Nanoparticle-based local antimicrobial drug delivery. Adv Drug Deliv Rev. 2018;127:46–57.
McBride AA, Warburton A. The role of integration in oncogenic progression of HPV-associated cancers. PLoS Pathog. 2017;13(4):e1006211.
Gilmore A: Anoikis. In: Nature Publishing Group; 2005.
Paoli P, Giannoni E, Chiarugi P. Anoikis molecular pathways and its role in cancer progression. Biochim Biophys Acta (BBA) Mol Cell Res. 2013;1833(12):3481–98.
Liu D, Huang C-l, Kameyama K, Hayashi E, Yamauchi A, Kobayashi S, Yokomise H. E-cadherin expression associated with differentiation and prognosis in patients with non–small cell lung cancer. Ann Thorac Surg. 2001;71(3):949–54.
Ohira T, Gemmill RM, Ferguson K, Kusy S, Roche J, Brambilla E, Zeng C, Baron A, Bemis L, Erickson P. WNT7a induces E-cadherin in lung cancer cells. Proc Natl Acad Sci. 2003;100(18):10429–34.
Zhang X, Liu G, Kang Y, Dong Z, Qian Q, Ma X. N-cadherin expression is associated with acquisition of EMT phenotype and with enhanced invasion in erlotinib-resistant lung cancer cell lines. PLoS One. 2013;8(3):e57692.
Yang H, Wang L, Zhao J, Chen Y, Lei Z, Liu X, Xia W, Guo L, Zhang H-T. TGF-β-activated SMAD3/4 complex transcriptionally upregulates N-cadherin expression in non-small cell lung cancer. Lung Cancer. 2015;87(3):249–57.
Alves CC, Carneiro F, Hoefler H, Becker K-F. Role of the epithelial-mesenchymal transition regulator slug in primary human cancers. Front Biosci. 2009;14(1):3035–50.
Shih J-Y, Tsai M-F, Chang T-H, Chang Y-L, Yuan A, Yu C-J, Lin S-B, Liou G-Y, Lee M-L, Chen JJ. Transcription repressor slug promotes carcinoma invasion and predicts outcome of patients with lung adenocarcinoma. Clin Cancer Res. 2005;11(22):8070–8.
Shioiri M, Shida T, Koda K, Oda K, Seike K, Nishimura M, Takano S, Miyazaki M. Slug expression is an independent prognostic parameter for poor survival in colorectal carcinoma patients. Br J Cancer. 2006;94(12):1816.
Walser T, Li R, Jing Z, Tran L, Dubinett S, Grimes B. Overexpression of slug drives malignant phenotypes in models of lung Premalignancy and Cancer. Am J Respir Crit Care Med. 2016;193:A3127.
Shih J-Y, Yang P-C. The EMT regulator slug and lung carcinogenesis. Carcinogenesis. 2011;32(9):1299–304.
Stone SC, Rossetti RAM, Lima AM, Lepique AP. HPV associated tumor cells control tumor microenvironment and leukocytosis in experimental models. Immun inflamm Dis. 2014;2(2):63–75.
Boccardo E, Lepique AP, Villa LL. The role of inflammation in HPV carcinogenesis. Carcinogenesis. 2010;31(11):1905–12.
Sadri Nahand J, Moghoofei M, Salmaninejad A, Bahmanpour Z, Karimzadeh M, Nasiri M, Mirzaei HR, Pourhanifeh MH, Bokharaei-Salim F, Mirzaei H. Pathogenic role of exosomes and microRNAs in HPV-mediated inflammation and cervical cancer: a review. Int J Cancer. 2020;146(2):305–20.
Gomes M, Teixeira AL, Coelho A, Araújo A, Medeiros R. The role of inflammation in lung cancer. In: Inflammation and cancer: Springer; 2014. p. 1–23.
Banat G-A, Tretyn A, Pullamsetti SS, Wilhelm J, Weigert A, Olesch C, Ebel K, Stiewe T, Grimminger F, Seeger W. Immune and inflammatory cell composition of human lung cancer stroma. PLoS One. 2015;10(9):e0139073.
Liu X, Ma X, Lei Z, Feng H, Wang S, Cen X, Gao S, Jiang Y, Jiang J, Chen Q. Chronic inflammation-related HPV: a driving force speeds oropharyngeal carcinogenesis. PLoS One. 2015;10(7):e0133681.
Huh WK. Human papillomavirus infection: a concise review of natural history. Obstet Gynecol. 2009;114(1):139–43.
Williams DA. Inflammatory cytokines and mucosal injury. JNCI Monographs. 2001;2001(29):26–30.
Hormia M, Willberg J, Ruokonen H, Syrjänen S. Marginal periodontium as a potential reservoir of human papillomavirus in oral mucosa. J Periodontol. 2005;76(3):358–63.
Xi LF, Koutsky LA, Castle PE, Edelstein ZR, Meyers C, Ho J, Schiffman M. Relationship between cigarette smoking and human papilloma virus types 16 and 18 DNA load. Cancer Epidemiol Prev Biomarkers. 2009;18(12):3490–6.
The authors are very grateful to the Director and staff of all involved Hospitals, all participant who agreed to participate in this study.
Ethics approval and consent to participate
Ethical approval for this study was obtained from Kermanshah University of Medical Sciences (KUMS) (IR.KUMS.REC.1399.095). Written informed consent form at the time of enrolment was obtained from each patient.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Rezaei, M., Mostafaei, S., Aghaei, A. et al. The association between HPV gene expression, inflammatory agents and cellular genes involved in EMT in lung cancer tissue. BMC Cancer 20, 916 (2020). https://doi.org/10.1186/s12885-020-07428-6
- Human papilloma virus
- Lung Cancer
- Tumour development
- Inflammatory cytokines
- Epithelial-mesenchymal transition (EMT)