PPAR Research
Volume 2020 (2020), Article ID 2510951, 9 pages
https://doi.org/10.1155/2020/2510951

Research Article

PPARG Could Work as a Valid Therapeutic Strategy for the Treatment of Lung Squamous Cell Carcinoma

Shunbin Shi,¹ Guiping Yu,² Bin Huang,² Yedong Mi,² Yan Kang,³ and Julia Pia Simon⁴

¹Department of Thoracic Surgery, Suzhou Ninth People’s Hospital, Suzhou, Jiangsu Province 215200, China
²Department of Cardiothoracic Surgery, Jiangyin Hospital of Southeast University Medical College, Jiangyin, Jiangsu Province 214400, China
³Department of General Surgery, Children’s National Medical Center, Washington, DC 20010, USA
⁴Neuroimaging and Informatics Institute, University of Southern California, California, Los Angeles 90007, USA

Correspondence should be addressed to Guiping Yu; xiaoyuer97103@163.com

Received 28 February 2020; Accepted 6 April 2020; Published 1 June 2020

Guest Editor: Fuquan Zhang

Copyright © 2020 Shunbin Shi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Previous studies showed that PPAR-gamma (PPARG) ligands might serve as potential therapeutic agents for nonsmall cell lung cancer (NSCLC). However, a few studies reported the specific relationship between PPARG and lung squamous cell carcinoma (LSCC). Here, we made an effort to explore the relationship between PPARG and LSCC. First, we used mega-analysis and partial mega-analysis to analyze the effects of PPARG on LSCC by using 12 independent LSCC expression datasets (285 healthy controls and 375 LSCC cases). Then, literature-based molecular pathways between PPARG and LSCC were established. After that, a gene set enrichment analysis (GSEA) was conducted to study the functionalities of PPARG and PPARG-driven triggers within the molecular pathways. Finally, another mega-analysis was constructed to test the expression changes of PPARG and its driven targets. The partial mega-analysis showed a significant downregulated expression of PPARG in LSCC (, ). Twelve diagnostic markers and four prognostic markers were identified within multiple PPARG-LSCC regulatory pathways. Our results suggested that the activation of PPARG expression may inhibit the development and progression of LSCC through the regulation of LSCC upstream regulators and downstream marker genes, which were involved in tumor cell proliferation and protein polyubiquitination/ubiquitination.

1. Introduction

PPARG is a ligand-activated transcription factor belonging to the family of peroxisome proliferator-activated receptors (PPARs) [1], which is widely expressed in many cells and tissues in the human body [2]. PPARG has recently attracted interest as the potential therapeutic target for a variety of malignancies [3]. A number of animal models [4], cell lines [5, 6], and clinical trials [NCT00923949, NCT01199068, NCT01199055] demonstrated that activation of PPAR-gamma impedes lung tumor progression and suggest that PPARG ligands may serve as potential therapeutic agents for nonsmall cell lung cancer (NSCLC), with the emphasis on lung adenocarcinoma [7, 8]. For instance, Ni et al.’s study showed that the activation of PPARG could inhibit the proliferation of EGFR-TKI-resistant lung adenocarcinoma cells and lead to a better survival rate [8].

So far, only a few studies explored the relationship between PPARG and lung squamous cell carcinomas (LSCC) [9, 10]. Kim et al. pointed out that the truncated splice variant of human PPAR gamma1 (hPPARG1(tr)) was strongly expressed in primary LSCC tumorous tissues, and the overexpression of hPPARG1(tr)) could increase the resistance of transfected cells to chemotherapeutic drug- and chemical-induced cell death [10]. However, to our knowledge, no study has systematically studied the role of PPARG in the pathology of LSCC.

To address this issue, we first conducted a meta-analysis to study the gene expression change of PPARG in the case of LSCC. Then, we integrated literature-based pathways and gene set enrichment analysis (GSEA) to study the potential pathways where PPARG could exert influence on the pathologic development of LSCC. Our results may help to understand the potential roles of PPARG in the case of LSCC.

2. Materials and Methods

2.1. Selection of LSCC Expression Data in Mega-Analysis

For the initial selection, we searched the LSCC expression datasets on gene expression omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) by using the keywords “lung squamous cell carcinoma.” Then, we applied the following criteria to filter further: (1) The entry type used in the study was series; (2) The dataset was array expression data; (3) The studies were designed as a comparison between LSCC and healthy controls: and (4) The organism of the dataset was Homo sapiens. Finally, we considered the datasets for which the original data and the corresponding format files were downloadable. Since we calculated the expression using the original data as extracted above, we used the term “mega-analysis” instead of “meta-analysis.”.

2.2. Mega-Analysis and Partial Mega-Analysis Models

In order to identify the relation between PPARG and LSCC, we used a mega-analysis to analyze the expression levels of PPARG in the case of LSCC. In our study, results from using both the random-effects model and the fixed-effects model were compared. To determine the heterogeneity of the datasets, between- and within-study variance was calculated and compared. When the total variance Q was no bigger than the expected value of the between-study variances (df), the model sets the ISq (percentage of the within- over between-study variance) to zero. In this case, the fixed-effects model, instead of the random-effects model, will be selected for the mega-analysis. All analyses were performed using Matlab (version R2017a; https://www.mathworks.com/products/matlab.html).

Moreover, we performed a partial mega-analysis to discover the significance of a gene presented in part of the studies/datasets (e.g., 50% of total studies) but not in all datasets, where 50% top studies/datasets were employed for the mega-analysis of a gene. Here, we define the “top datasets” for a gene as these datasets that demonstrate the bigger absolute value of effect size than the rest of the datasets. It should be mentioned that the top datasets for different genes could be different.

2.3. Analysis of Influential Factors

To estimate the possible influence of several factors (e.g., study date, country of origin, and sample size) on the gene expression in the case of MI, we conducted a multiple linear regression (MLR) analysis and reported the values for each of these factors.

2.4. Construct PPARG-Drive Network and Gene Set Enrichment Analysis

Based on large scale literature data mining, we constructed a diagnostic and a prognostic functional network connecting PPARG and LSCC. In the diagnostic network, we identified the genes that were contra-directionally regulated by PPARG and LSCC. To achieve this goal, we used Pathway Studio (http://www.pathwaystudio.com/) to identify PPARG➔gene relationships and LSCC➔gene relationships with polarity. Each of these relationships has supported by one or two scientific references. Then we identified the overlapped genes within these relationships to construct the PPARG—LSCC diagnostic network. To increase the reliability of the identified network, we limited the genes to these that also demonstrated consistency in terms of their gene expression alteration in the case of LSCC in the mega-analysis. For the prognostic pathway, we followed the same pressure but to identify genes that were downstream targets of PPARG and upstream regulators of LSCC. The reference information supporting the relations identified in these networks was provided in the Supplementary Materials (available here), including the type of the relationship, supporting references, and related sentences from the references where the relationship has been identified.

For these genes within the diagnostic and prognostic networks built above, a gene set enrichment analysis (GSEA) was conducted using Pathway Studio (version 12.1.0.9; http://www.pathwaystudio.com/) against Gene Ontology (GO; http://geneontology.org) and Pathway Studio pathways. The purpose of GSEA was to test the functional profile of the genes involved in the PPARG-driven networks.

3. Results

3.1. Mega-Analysis Based on the Selected LSCC Expression Datasets

There were 4,643 results shown in the GEO datasets identified by the keywords “lung squamous cell carcinoma”. Then, further filters were set as our criteria. A total of 12 datasets satisfied the inclusion criteria for the mega-analysis, which are listed in Table 1. The studies were distributed in 8 different countries, and the study dates ranged from 2 to 15 years ago, including 285 healthy controls and 375 LSCC cases.

Table 1: The 12 qualified LSCC expression datasets for mega-analysis.

GEO ID
Control ()
Case ()
Country
Study age
Sample organism

GSE84784
9
9
Luxembourg
2
Homo sapiens

GSE67061
8
69
China
3
Homo sapiens

GSE30219
14
61
France
5
Homo sapiens

GSE33479
27
14
USA
5
Homo sapiens

GSE32036
59
12
USA
7
Homo sapiens

GSE19188
65
27
Netherlands
9
Homo sapiens

GSE11969
5
35
Japan
10
Homo sapiens

GSE2088
30
48
Japan
10
Homo sapiens

GSE12428
28
34
Netherlands
11
Homo sapiens

GSE6044
5
14
Germany
13
Homo sapiens

GSE1987
7
17
Israel
15
Homo sapiens

GSE12472
28
35
Netherland
10
Homo sapiens

The mega-analysis and partial mega-analysis results for gene PPARG are presented in Table 2. As shown in Figure 1(a), the total variance (Q) was larger than the expected between-study variance (), the within-study variance percentage () was 45.10, the between-study variances were significant, and thus a random-effects model was selected for PPARG in the mega-analysis. However, there were no significant between-study variances (, Q test ), see in Figure 1(b). Thus, the fixed-effects model was selected for PPARG in partial mega-analysis. The LFC of the gene was estimated from about half (top 50%) of the selected datasets. PPARG demonstrated significantly lower expression in the case of LSCC (, ).

Table 2: Analysis of PPARG expression levels in LSCC datasets.

PPARG
Mega-analysis
Partial mega-analysis

Models
Random effects
Fixed effects

# study
11
5

LFC (effect size)
-0.22
-1.08

value
0.094
0.00073

ISq (%)
45.10
0.00

value Q
0.051
0.64

# sample
0.85
0.10

Country
0.0045
2.79e-6

Study age
0.30
0.00034

(a)

(b)

(c)

Figure 1: Mega-analysis results of PPARG using 11 LSCC RNA expression datasets. (a) Mega-analysis results from the random-effects model. (b) Partial mega-analysis results from the fixed-effects model. (c) The influence of population region (country) on the PPARG expression levels. The bar plot on the right of each figure represents the normalized weights for each dataset/study, range (0, 1); the brighter (green) the color, the bigger the weight (labeled right next to the bar). The star (in red) and lines (in blue) on the left are the mean of effect size (log fold change) and 95% confidence interval (CI) of each dataset/study, respectively.

MLR analysis showed that sample size and study age were not significant influential factors for the expression levels of PPARG among the 11 LSCC datasets (). However, the population region (country) was identified as a significant factor that influences the LFC of PPARG in the case of LSCC (, Figure 1(c)). This may partially explain the differential results between the partial mega-analysis and the mega-analysis.

3.2. The LSCC Diagnostic Network Interfered by PPARG

Multiple molecules (12 genes) have been identified through large-scale literature data mining that was contra-directionally influenced by PPARG and LSCC, as shown in Figure 2. According to previous literature reports, a total of 8 molecules (XIAP, UBE2D1, SKP2, ACKR3, MI21, HOXA10, STAT1, and PDPN) were upregulated in LSCC (the genes at the bottom of Figure 2; the arrows with ⊕, Figure 2) but negatively affected by PPARG (the arrows with ┥, Figure 2). These eight genes also presented increased expression levels in the 12 LSCC RNA expression datasets. These results support the literature data mining results and suggest these eight genes as positive markers for LSCC. The mega-analysis results for these genes were provided in the Supplementary Materials→Mega-analysis (available here). The inhibition of these genes by PPARG could exert an anti-LSCC effect during its pathological development (Figure 2).

Figure 2: LSCC diagnostic network interfered with PPARG. Genes highlighted in blue (genes at the top of the figure) were literature-implicated with a downregulation in the case of LSCC, and those highlighted in red (genes at the bottom of the figure) were upregulated according to literature reports. Genes in blue represent a decreased expression level from the mega-analysis using 12 LSCC datasets, while those in red represent an increased expression level.

On the other hand, four genes (MIR223, ANGPT1, CYP2A6, and FOXA2) have been suggested to get suppressed in LSCC (the genes at the top of Figure 2, the arrows with ┥) but were stimulated by PPARG (Figure 2). These four genes also presented decreased expression levels in the mega-analysis, supporting the literature data mining results. Activation of these molecules could be due to other pathways where PPARG inhibits the progress of LA. Detailed information regarding the network presented in Figure 2 can be found in the Supplementary Materials→LSCC (available here) diagnostic network, including the type of the relationships and the supporting references.

3.3. GSEA for the Genes Involved in LSCC Diagnostic Network

The GSEA was performed using Pathway Studio with the purpose of investigating the biological functions of the 12 genes within the LSCC diagnostic network. The GSEA was also confirmed by the mega-analysis, including eight upregulated and five suppressed genes. A total of 10 out of these 12 genes were shared among the top 10 most significantly enriched pathways (, for FDR), which are presented in Table 3. The full 21 pathways/gene sets enriched with were presented in the Supplementary Materials→GSEA1 (available here). Notably, enriched pathways highlighted by the GSEA approach are mainly related to the regulation of protein ubiquitination, the regulation of cell proliferation, and cytokine stimulus.

Table 3: The top 10 genetic pathways enriched by the 12 genes within LSCC_diagnostic network.

Name
GO ID
Overlap
value
Jaccard similarity

GO: positive regulation of protein polyubiquitination
1902916;
3
0.00085
0.14

GO: regulation of smooth muscle cell proliferation
0048660;
5
0.0012
0.022

GO: regulation of protein polyubiquitination
1902914;
3
0.0015
0.097

GO: positive regulation of protein ubiquitination
0031398;
4
0.0036
0.029

GO: positive regulation of protein modification by small protein conjugation or removal
1903322;
4
0.0050
0.025

GO: regulation of response to cytokine stimulus
0060759;
4
0.0052
0.022

GO: regulation of cytokine-mediated signaling pathway
0001959;
4
0.0052
0.023

GO: regulation of endothelial cell proliferation
0001936;
4
0.0080
0.019

GO: regulation of protein ubiquitination
0031396;
4
0.010
0.017

GO: regulation of response to external stimulus
0032101;
6
0.011
0.0060

3.4. LSCC Prognostic Network Interfered by PPARG

As shown in Figure 3, a regulatory pathway connecting PPARG and LSCC was identified, heavily involved in the pathological development of LSCC. Based on the literature reports, three genes, TNF, NOS2, and ACE, could promote the pathological development of LSCC (highlighted in red, the arrows with ⊕, Figure 3). These three genes were deactivated by PPARG. However, according to the mega-analysis results, these genes were downregulated together with PPARG in the case of LSCC. Therefore, PPARG may not necessarily be needed for the deactivation of these genes to inhibit the progress of LSCC. The supporting references for each relation presented in Figure 3 were provided in the Supplementary Materials→LSCC_ (available here)prognostic network, which also include the type of the relationships.

Figure 3: LSCC prognostic network interfered with PPARG. Genes in blue represent a decreased expression level from the mega-analysis using 12 LSCC datasets; entities in red represent an increased expression level. Entities highlighted in blue were literature implicated as the LSCC-inhibitors, while those highlighted in red were LSCC-promoters.

In addition, an LSCC-inhibitor, STK11, has been shown to be activated by PPARG (highlighted in blue, the arrows with ┥, see Figure 3(a)). Mega-analysis showed that this gene showed slightly increased expression in the case of LSCC. Therefore, activation of PPARG may further promote the activation of STK11, which could be a blocker for the pathologic development of LSCC. The mega-analysis results of these four genes in Figure 3 were provided in the Supplementary Materials→Mega-analysis (available here).

3.5. GSEA for the Genes Involved in LSCC Prognostic Network

GSEA results showed that the five genes (PPARG, STK11, NOS2, and TNF) were significantly enriched within 59 pathways/gene sets (; for FDR; see Supplementary Materials→GSEA2 (available here)). We presented the top 10 pathways (4 genes were enriched; ) in Table 4. The pathways were mainly related to cell metabolism and hormone level regulation, which was largely different from that of the LSCC diagnostic network. Also, notably, these five genes enrich more pathways than that of the 12 genes within the LSCC diagnostic network, indicating that these five genes were more functionally linked to each other.

Table 4: The top 10 genetic pathways enriched by the 12 genes within LSCC_diagnostic network.

Name
GO ID
Overlap
value
Jaccard similarity

GO: positive regulation of small molecule metabolic process
0062013
4
0.00054
0.023

GO: regulation of small molecule metabolic process
0062012
4
0.0036
0.0095

GO: fatty acid transport
0015908
3
0.0036
0.037

GO: receptor biosynthetic process
0032800
2
0.0036
0.22

GO: response to oxygen levels
0070482
4
0.0049
0.0073

GO: regulation of inflammatory response
0050727
4
0.0049
0.0079

GO: cellular response to organic cyclic compound
0071407
4
0.0076
0.0063

GO: regulation of hormone levels
0010817
4
0.0077
0.0060

GO: monocarboxylic acid transport
0015718
3
0.0077
0.019

CSF2 -> STAT expression targets
NONE
4
0.0077
0.055

4. Discussion

Previous studies suggested that the activation of PPARG might be associated with the inhibition of NSCLC [4–6]. However, most studies were focused on the cases of lung adenocarcinoma [4–6]. In this study, we aim to explore the possible linkage between PPARG and LSCC. First, we utilized a mega-analysis and a partial mega-analysis to analyze the potential relationship between PPARG and LSCC. Subsequently, we integrated knowledge from large-scale literature data mining and existing LSCC expression data to construct molecular networks connecting PPARG and LSCC, followed by a GAEA analysis to study the functional profile of the molecules involved in the PPARG-drive network. Our results showed that PPARG was significantly downregulated in about half of the LSCC cases, with multiple pathways suggesting an inhibition role of PPARG in the pathologic development and progress of LSCC.

Notably, PPARG did not show a significant decrease in the 11 studies overall (; ), while it demonstrated significant decreased expression in 5 out of the 11 studies (; ). These results indicate that there are influential factors that lead to different expression levels of PPARG among different studies. MLR analysis showed that the population region was a significant factor that influenced the PPARG levels (Figure 1(c)). Specifically, PPARG demonstrated low expression of in the dataset from France (GSE30219), while relatively high expression in the dataset from Germany (; GSE6044). Notably, one of the highest (GSE12428; ) and the lowest (GSE19188; ) expression levels were from the Netherlands, indicating that, besides population region, there could be other factors that influence the PPARG expression levels in LSCC patients. Further investigation showed that the GSE12428 dataset was especially studying LSCC patients who smoke, while the dataset of GSE19188 contained all patients in general. It has been shown that smoking is one of the major risk factors for the development of LSCC [11], and smokers tend to have low PPARG expression levels [12]. These studies may explain the different PPARG levels in these two Netherlands datasets. Moreover, MLR results also suggest that sample size and study date were not significant factors for the PPARG levels. Due to lack of data, we only studied the influence of three factors on the expression levels of PPARG. Further study is needed to test the expression level of PPARG in LSCC patients and the possible influential factors such as age, gender, and complication.

Functional network analysis showed that PPARG could play roles both in the development and progression of LSCC. Specifically, PPARG counter-regulated 12 molecules that were upregulated or downregulated by LSCC (Figure 2). The expression levels of eight LSCC markers (XIAP, UBE2D1, SKP2, ACKR3, MI21, HOXA10, STAT1, and PDPN) were significantly upregulated in LSCC patients [13–20] and were downregulated by PPARG [21–28]. On the other hand, PPARG could activate multiple genes [29–33] inhibited by LSCC [34–38], including MIR 223, PTEN, ANGPT1, CYP2A6, and FOXA2. It should also be noted that the relations between LSCC and the molecules in the network are supported by both the literature data mining and mega-analysis using 12 LSCC expression datasets, which strengthens the reliability of the PPARG-driven network.

GSEA analysis suggested that PPARG may influence the development of LSCC through multiple pathways (Table 3 and Table 4). Besides the regulation of cell proliferation associated with the procession of LSCC, PPARG could also regulate protein polyubiquitination and ubiquitination, which has become increasingly recognized as a controller to regulate the function and signaling of a profusion of proteins. Ubiquitination affects proteins in various cellular processes, including signal transduction, DNA repair, chromosome maintenance, transcriptional activation, cell cycle progression, cell survival, and certain immune cell functions [39]. Thus, it is not surprising that ubiquitin metabolism enzymes prominently feature either oncogenes or tumor suppressors in a variety of cancers and many pathways relevant to cancer. A previous study suggested that targeting those physiological processes may effectively abate the proliferation and facilitate the treatment of lung cancer cells [40]. In the LSCC diagnostic network (Figure 2), E2 ubiquitin-conjugating enzymes (UBE2D1), E3 ubiquitin ligases (XIAP), and SKP2 were involved in the regulation of protein ubiquitination. The expression of XIAP, UBE2D1, and SKP2 downregulated by PPARG at the transcriptional level [21, 22, 41] were overexpressed in LSCC tissues with robust proliferation ability [13–15].

PPARG may also play a role in the progression of LSCC by interfering upstream regulators of LSCC, as shown in Figure 3. For instance, the knock-down of PPARG has been shown to deactivate STK11 [42], while the loss of STK11 could lead to the formation of LSCC [43]. Furthermore, the activation of PPARG could inhibit three promoters of LSCC, including NOS2 [44], ACE [45], and TNF [46]. Thus, increased expression of PPARG may inhibit the formation of LSCC.

The most significant contribution of this study was the identification of the two PPARG-driven networks (Figures 2 and 3) that partially explain the mechanism of the roles of PPARG in the etiology and development of LSCC. However, the integration of literature data-mine and mega-analysis may exclude potential genes/molecules connection PPARG and LSCC. Further study is needed to validate and consummate the networks identified here.

5. Conclusion

Results from this study indicated that the expression of PPARG might be suppressed in LSCC patients. Activation of PPARG expression may inhibit the development and progress of LSCC through the regulation of LSCC upstream regulators and downstream marker genes. Our results indicate the need for further study of the relationship between PPARG and LSCC.

Data Availability

The data of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

All the authors declare no conflict of interest.

Authors’ Contributions

Shunbin Shi and Guiping Yu have contributed equally to this work.

Supplementary Materials

(1) GSEA1 and GAEA2 present GSEA results for the genes within Figure 2 and Figure 3, respectively. (2) Mega-analysis presents the mega-analysis results for PPARG gene and its driven genes within the Figure 2 and Figure 3. (3) Partial_Mega-analysis presents the partial mega-analysis results for PPARG; (4) LSCC_diagnostic network presents the reference information for the network presented in Figure 2. (4) LSCC_prognostic network presents the reference information for the network presented in Figure 3. (Supplementary Materials)

References

V. Bocher, G. Chinetti, J. C. Fruchart, and B. Staels, “Role of the peroxisome proliferator-activated receptors (PPARS) in the regulation of lipids and inflammation control,” Journal de la Societe de biologie, vol. 196, no. 1, pp. 47–52, 2002.

J. Berger and D. E. Moller, “The mechanisms of action of PPARs,” Annual Review of Medicine, vol. 53, no. 1, pp. 409–435, 2002.

J. M. Peters, Y. M. Shah, and F. J. Gonzalez, “The role of peroxisome proliferator-activated receptors in carcinogenesis and chemoprevention,” Nature Reviews Cancer, vol. 12, no. 3, pp. 181–195, 2012.

N. Skrypnyk, X. Chen, W. Hu et al., “PPARα activation can help prevent and treat non-small cell lung cancer,” Cancer Research, vol. 74, no. 2, pp. 621–631, 2014.

V. G. Keshamouni, R. C. Reddy, D. A. Arenberg et al., “Peroxisome proliferator-activated receptor- γ activation inhibits tumor progression in non-small-cell lung cancer,” Oncogene, vol. 23, no. 1, pp. 100–108, 2004.

S. Kaur, A. Nag, G. Gangenahalli, and K. Sharma, “Peroxisome Proliferator Activated Receptor Gamma Sensitizes Non-small Cell Lung Carcinoma to Gamma Irradiation Induced Apoptosis,” Front Genet, vol. 10, p. 554, 2019.

Y. Susaki, M. Inoue, M. Minami et al., “Inhibitory effect of PPARγ on NR0B1 in tumorigenesis of lung adenocarcinoma,” International Journal of Oncology, vol. 41, no. 4, pp. 1278–1284, 2012.

J. Ni, L. L. Zhou, L. Ding et al., “PPARγ agonist efatutazone and gefitinib synergistically inhibit the proliferation of EGFR-TKI-resistant lung adenocarcinoma cells via the PPARγ/PTEN/Akt pathway,” Experimental Cell Research, vol. 361, no. 2, pp. 246–256, 2017.

H. Li, M. C. Weiser-Evans, and R. Nemenoff, “Anti- and Protumorigenic Effects of PPARγ in Lung Cancer Progression: A Double-Edged Sword,” PPAR Research, vol. 2012, Article ID 362085, 12 pages, 2012.

H. J. Kim, J. Y. Hwang, H. J. Kim et al., “Expression of a peroxisome proliferator-activated receptor gamma 1 splice variant that was identified in human lung cancers suppresses cell death induced by cisplatin and oxidative stress,” Clinical Cancer Research, vol. 13, no. 9, pp. 2577–2583, 2007.

M. C. Boelens, A. van den Berg, R. S. Fehrmann et al., “Current smoking-specific gene expression signature in normal bronchial epithelium is enhanced in squamous cell lung cancer,” The Journal of Pathology, vol. 218, no. 2, pp. 182–191, 2009.

B. C. Bergman, L. Perreault, D. M. Hunerdosse, M. C. Koehler, A. M. Samek, and R. H. Eckel, “Intramuscular lipid metabolism in the insulin resistance of smoking,” Diabetes, vol. 58, no. 10, pp. 2220–2227, 2009.

Y. B. Chen, J. Shu, W. T. Yang et al., “XAF1 as a prognostic biomarker and therapeutic target in squamous cell lung cancer,” Chinese Medical Journal, vol. 124, no. 20, pp. 3238–3243, 2011.

L. Hou, Y. Li, Y. Wang et al., “UBE2D1 RNA Expression Was an Independent Unfavorable Prognostic Indicator in Lung Adenocarcinoma, but Not in Lung Squamous Cell Carcinoma,” Disease Markers, vol. 2018, Article ID 4108919, 8 pages, 2018.

K. Zhong, F. Yang, Q. Han, J. Chen, and J. Wang, “Skp2 expression has different clinicopathological and prognostic implications in lung adenocarcinoma and squamous cell carcinoma,” Oncology Letters, vol. 16, no. 3, pp. 2873–2880, 2018.

B. Behnam Azad, A. Lisok, S. Chatterjee et al., “Targeted Imaging of the Atypical Chemokine Receptor 3 (ACKR3/CXCR7) in Human Cancer Xenografts,” Journal of Nuclear Medicine, vol. 57, no. 6, pp. 981–988, 2016.

Y. Song, H. Dou, P. Wang et al., “A novel small-molecule compound diaporine A inhibits non-small cell lung cancer growth by regulating miR-99a/mTOR signaling,” Cancer Biology Therapy, vol. 15, no. 10, pp. 1423–1430, 2014.

A. Clemenceau, O. Boucherat, K. Landry-Truchon et al., “Lung cancer susceptibility genetic variants modulate HOXB2 expression in the lung,” The International Journal of Developmental Biology, vol. 62, no. 11-12, pp. 857–864, 2018.

M. Yang, H. Chen, L. Zhou, K. Chen, and F. Su, “Expression profile and prognostic values of STAT family members in non-small cell lung cancer,” American journal of translational research, vol. 11, no. 8, pp. 4866–4880, 2019.

H. Suzuki, M. Onimaru, T. Koga et al., “High podoplanin expression in cancer cells predicts lower incidence of nodal metastasis in patients with lung squamous cell carcinoma,” Pathology, Research and Practice, vol. 207, no. 2, pp. 111–115, 2011.

K. Bräutigam, J. Biernath-Wüpping, D. O. Bauerschlag et al., “Combined treatment with TRAIL and PPARγ ligands overcomes chemoresistance of ovarian cancer cell lines,” Journal of Cancer Research and Clinical Oncology, vol. 137, no. 5, pp. 875–886, 2011.

P. E. Almeida, A. B. Carneiro, A. R. Silva, and P. T. Bozza, “PPARγ Expression and Function in Mycobacterial Infection: Roles in Lipid Metabolism, Immunity, and Bacterial Killing,” PPAR Research, vol. 2012, 383827 pages, 2012.

J. Meng, Y. Ding, A. Shen et al., “Overexpression of PPARγ can down-regulate Skp2 expression in MDA-MB-231 breast tumor cells,” Molecular and Cellular Biochemistry, vol. 345, no. 1-2, pp. 171–180, 2010.

D. Zhao, Z. Zhu, D. Li, R. Xu, T. Wang, and K. Liu, “Pioglitazone Suppresses CXCR7 Expression To Inhibit Human Macrophage Chemotaxis through Peroxisome Proliferator-Activated Receptor γ,” Biochemistry, vol. 54, no. 45, pp. 6806–6814, 2015.

Y. F. Zhang, H. M. Xu, F. Yu et al., “Crosstalk between MicroRNAs and Peroxisome Proliferator-Activated Receptors and Their Emerging Regulatory Roles in Cardiovascular Pathophysiology,” PPAR Research, vol. 2018, Article ID 8530371, 11 pages, 2018.

R. Yasmeen, J. M. Meyers, C. E. Alvarez et al., “Aldehyde dehydrogenase-1a1 induces oncogene suppressor genes in B cell populations,” Biochimica et Biophysica Acta (BBA) - Molecular Cell Research, vol. 1833, no. 12, pp. 3218–3227, 2013.

V. Costa, M. A. Gallo, F. Letizia, M. Aprile, A. Casamassimi, and A. Ciccodicola, “PPARG: Gene Expression Regulation and Next-Generation Sequencing for Unsolved Issues,” PPAR Research, vol. 2010, Article ID 409168, 17 pages, 2010.

X. Zhu, Y. Chen, W. Zhu et al., “Oroxylin A inhibits Kaposi’s sarcoma-associated herpes virus (KSHV) vIL-6-mediated lymphatic reprogramming of vascular endothelial cells through modulating PPARγ/Prox1 axis,,” Journal of Medical Virology, vol. 91, no. 3, pp. 463–472, 2019.

W. Ying, A. Tseng, R. C. A. Chang et al., “MicroRNA-223 is a crucial mediator of PPARγ-regulated alternative macrophage activation,” Journal of Clinical Investigation, vol. 125, no. 11, pp. 4149–4159, 2015.

B. Farrow and B. M. Evers, “Activation of PPARγ increases PTEN expression in pancreatic cancer cells,” Biochemical and Biophysical Research Communications, vol. 301, no. 1, pp. 50–53, 2003.

T. Fink, L. Abildtrup, K. Fogd et al., “Induction of adipocyte-like phenotype in human mesenchymal stem cells by hypoxia,” Stem Cells, vol. 22, no. 7, pp. 1346–1355, 2004.

J. Deng, L. Guo, and B. Wu, “Circadian Regulation of Hepatic Cytochrome P450 2a5 by Peroxisome Proliferator-Activated Receptorγ,” Drug Metabolism and Disposition, vol. 46, no. 11, pp. 1538–1545, 2018.

H. S. Kim, Y. C. Hwang, S. H. Koo et al., “PPAR-γ activation increases insulin secretion through the up-regulation of the free fatty acid receptor GPR40 in pancreatic β-cells,” PLoS One, vol. 8, no. 1, article e50128, 2013.

P. Luo, Q. Wang, Y. Ye et al., “MiR-223-3p functions as a tumor suppressor in lung squamous cell carcinoma by miR-223-3p-mutant p53 regulatory feedback loop,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 74, 2019.

R. D. Jolly, K. G. Thompson, and B. G. Winchester, “Bovine mannosidosis--a model lysosomal storage disease,” Bovine mannosidosis--a model lysosomal storage disease, vol. 11, no. 6, pp. 273–278, 1975.

S. Yao, S. S. Dong, J. M. Ding et al., “Sex-specific SNP-SNP interaction analyses within topologically associated domains reveals ANGPT1 as a novel tumor suppressor gene for lung cancer,” Genes Chromosomes Cancer, vol. 59, no. 1, pp. 13–22, 2019.

H. Hosono, M. Kumondai, T. Arai et al., “CYP2A6 genetic polymorphism is associated with decreased susceptibility to squamous cell lung cancer in Japanese smokers,” Drug Metabolism and Pharmacokinetics, vol. 30, no. 4, pp. 263–268, 2015.

L. Tang, J. Liu, L. Zhu et al., “Curcumin Inhibits Growth of Human NCI-H292 Lung Squamous Cell Carcinoma Cells by Increasing FOXA2 Expression,” Frontiers in Pharmacology, vol. 9, p. 60, 2018.

W. H. McBride, K. S. Iwamoto, R. Syljuasen, M. Pervan, and F. Pajonk, “The role of the ubiquitin/proteasome system in cellular responses to radiation,” Oncogene, vol. 22, no. 37, pp. 5755–5773, 2003.

Y. Tang, Y. Geng, J. Luo et al., “Downregulation of ubiquitin inhibits the proliferation and radioresistance of non-small cell lung cancer cells _in vitro_ and _in vivo_,” Scientific Reports, vol. 5, no. 1, p. 9476, 2015.

S. T. Wang, H. J. Ho, J. T. Lin, J. J. Shieh, and C. Y. Wu, “Simvastatin-induced cell cycle arrest through inhibition of STAT3/SKP2 axis and activation of AMPK to promote p27 and p21 accumulation in hepatocellular carcinoma cells,” Cell Death Disease, vol. 8, no. 2, p. e2626, 2017.

J. Ji, T. F. Xue, X. D. Guo et al., “Antagonizing peroxisome proliferator-activated receptor γ facilitates M1-to-M2 shift of microglia by enhancing autophagy via the LKB1-AMPK signaling pathway,” Aging Cell, vol. 17, no. 4, p. e12774, 2018.

C. Xu, C. M. Fillmore, S. Koyama et al., “Loss of Lkb1 and Pten leads to lung squamous cell carcinoma with elevated PD-L1 expression,” Cancer Cell, vol. 25, no. 5, pp. 590–604, 2014.

T. E. Cullingford, “The ketogenic diet; fatty acids, fatty acid-activated receptors and neurological disorders,” Prostaglandins, Leukotrienes, and Essential Fatty Acids, vol. 70, no. 3, pp. 253–264, 2004.

A. Vallée, B. L. Lévy, and J. Blacher, “Interplay between the renin-angiotensin system, the canonical WNT/β-catenin pathway and PPARγ in hypertension,” Current Hypertension Reports, vol. 20, no. 7, p. 62, 2018.

W. A. Hsueh and R. Law, “The central role of fat and effect of peroxisome proliferator-activated receptor–γ on progression of insulin resistance and cardiovascular disease,” The American Journal of Cardiology, vol. 92, no. 4A, pp. 3J–9J, 2003.