Biological Network Analysis of PPAR and Related Signaling PathwaysView this Special Issue
Research Article | Open Access
Min Zhao, Xiaoyang Li, Yunxiang Zhang, Hongming Zhu, Zhaoqing Han, Yan Kang, "PPARG Drives Molecular Networks as an Inhibitor for the Pathologic Development and Progression of Lung Adenocarcinoma", PPAR Research, vol. 2020, Article ID 6287468, 7 pages, 2020. https://doi.org/10.1155/2020/6287468
PPARG Drives Molecular Networks as an Inhibitor for the Pathologic Development and Progression of Lung Adenocarcinoma
Previous studies showed that low PPARG expression was associated with poor prognosis of lung adenocarcinoma (LA) with limited mechanisms identified. We first conducted a large-scale literature-based data mining to identify potential molecular pathways where PPARG could exert influence on the pathological development of LA. Then a mega-analysis using 13 independent LA expression datasets and a Pathway Enrichment Analysis (PEA) was conducted to study the gene expression levels and the functionalities of PPARG and the PPARG-driven triggers within the molecular pathways. Finally, a protein-protein interaction (PPI) network was established to reveal the functional connection between PPARG and its driven molecules. We identified 25 PPARG-driven molecule triggers forming multiple LA-regulatory pathways. Mega-analysis using 13 LA datasets supported these pathways and confirmed the downregulation of PPARG in the case of LA (). Results from the PEA and PPI analysis suggested that PPARG might inhibit the development of LA through the regulation of tumor cell proliferation and transmission-related molecules, including an LA tumor cell suppressor MIR145. Our results suggested that increased expression of PPARG could drive multiple molecular triggers against the pathologic development and prognosis of LA, indicating PPARG as a valuable therapeutic target for LA treatment.
Lung adenocarcinoma is one of the most common histological subtypes of Nonsmall-cell lung carcinoma that accounts for about 85% of all cases of lung cancer worldwide . The overall 5-year survival rate of lung cancer is low even after surgical treatment (about 69.6%) . Therefore, more effective strategies of therapy are necessary.
Locating on chromosome 3 (base pairs 12,287,485 to 12,434,356), PPARG encodes a member of the peroxisome proliferator-activated receptor (PPAR) subfamily of nuclear receptors—Peroxisome proliferator-activated receptor γ (PPARγ), which have been shown to possesses an antagonistic function against LA (PMID: 22843091). However, so far, limited knowledge of this PPARG-inhibiting-LA mechanism is known . On the one hand, it has been shown that the expression of PPARG was reduced in LA progression cells , and the low PPARG expression was strictly correlated with poor prognosis of stage IA LA . On the other hand, increased expression of PPARG has been positively associated with a better survival rate of LA patients . Ni et al.’s work showed that overexpression of PPARγ could inhibit the drug resistance effect of gefitinib in the treatment of LA by reducing the proliferation of gefitinib-resistant cells . Susaki et al. showed that PPARG could inhibit the tumorigenic potential of NR0B1 in LA . However, more studies are needed to identify the underlying mechanism of the role that PPARG plays in the pathological development of LA.
To dissect the role of PPARG in LA at the genetic level, we employed Pathway Studio (http://www.pathwaystudio.com) knowledge database to undertake large-scale literature mining effort and integrated its results with an analysis of multiple LA expression datasets. We identified a set of PPARG-driven molecular triggers, possibly contributing to inhibition of the development of LA through a quantitative regulation. Our results might add new insights into the understanding of the LA-inhibition role of PPARG.
2. Materials and Method
This study is organized as follows. First, a large-scale literature-based data mining was performed to identify genes as the disease markers and the regulators of LA. Subsequently, regulations of PPARG on these LA genes were identified under the Pathway Studio environment. After that, a mega-analysis was performed using 13 independent LA gene expression datasets to test the expression changes of PPARG and the LA genes that were regulated by PPARG. Finally, a Pathway Enrichment Analysis (PEA) has been conducted to explore the functionality of the PPARG-driven molecular triggers, with protein-protein interaction (PPI) network built. All data and analysis results were organized in an excel file named as PPARG_LA, which is downloadable at http://www.gousinfo.com/database/Data_Genetic/PPARG_LA.xlsx.
2.1. Literature-Based Pathway Analysis
Assisted by using Pathway Studio (http://www.pathwaystudio.com), we conducted a large-scale literature-based functional pathway analysis to investigate the potential biological associations between PPARG and LA. Specifically, we identified the genes influenced by PPARG and also regulating LA to build the connections between PPARG and LA. Only relationships with polarity were selected within the Pathway Studio database. Each and all of the relations identified were supported by one or more references (1255 references in total; please refer to the worksheet “PPARG-LA Regulation Pathway” within the file PPARG_LA). In the PPARG_LA→PPARG-LA Regulation Pathway, the reference information supporting the relations identified in the PPARG-LA regulatory pathways was provided, including the types of associations, the number of underlying supporting references, and the sentences where these associations had been identified and described. The expression changes of PPARG and its driven genes involved in the pathways were tested using a mega-analysis approach described as follows.
2.2. Gene Expression Data Selected for Mega-Analysis
Following the initial search with “Lung adenocarcinoma”, 634 microarray expression datasets were identified on gene expression omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) . Subsequently, the following criteria were applied: (1) the organism used in the study was Homo sapiens; (2) the data type was microarray expression profiling; (3) the studies were limited to a comparison between LA and healthy controls; and (4) the original data and the corresponding format file were downloadable. A total of 13 datasets satisfied the inclusion criteria for the mega-analysis, which are listed in Table 1.
2.3. Mega-Analysis Models
For PPARG and the 25 genes involved in the PPARG-LA regulatory pathway, the log2 fold-change (LFC) of the gene expression level was used to indicate the effect size. Both fixed-effects and random-effects models were employed to investigate and compare the effect size (doi:10.1002/jrsm.12). The heterogeneity of the mega-analysis was analyzed to study the variance within and between different studies. In the case that the total variance (Q) was equal to or smaller than the expected between-study variance (df), the was set at 0, and a fixed-effects model was selected for the mega-analysis. Otherwise, a random-effects model was selected. Q-p represents the probability that the total variance was only due to within-study variance. The current study presented all the mega-analysis results in the worksheet “mega-analysis” of the excel file PPARG_LA (http://www.gousinfo.com/database/Data_Genetic/PPARG_LA). The full name and description of related stats were given within “mega-analysis.” All analyses were performed using Matlab (version R2017a; https://www.mathworks.com/products/matlab.html). Here, we used the term “mega-analysis” instead of “meta-analysis” to reflect that the LFC of each gene was calculated from the original data instead of results within publications.
2.4. Multiple Linear Regression Analysis
A multiple linear regression (MLR) model was employed to investigate the possible influence of sample size, country of origin, and study date on the gene expression in the case of LA. values were reported for each of these factors.
2.5. Pathway Enrichment Analysis and Protein-Protein Interaction Analysis
To test the functional profile of the genes involved in the PPARG-LA regulation pathway, a Fisher’s Exact Test based pathway enrichment analysis (PEA; Pathway Studio: Find Pathways/Groups Enriched with Selected Entities) was conducted using Pathway Studio (version 220.127.116.11; http://www.pathwaystudio.com) against Gene Ontology (GO; http://geneontology.org) and Pathway Studio pathways. Statistics for the enriched pathways were provided, including false discovery rate (FDR) corrected -value and Jaccard similarity.
Moreover, based on the PEA results, a protein-protein interaction (PPI) network was constructed. Two genes were recognized as connected if they were identified to play roles within at least one common pathway (from Pathway Studio Pathway Collection) or functional group (from GO groups).
3.1. PPARG-LA Contradirectional Common Targets
Pathway analysis has identified multiple molecules that were contradirectionally influenced by PPARG and LA, as shown in Figure 1. The expression levels of these genes, including PPARG, were tested in the mega-analysis and color-coded with the literature-based pathway (see the color bar in Figure 1).
According to previous literature reports, a total of 13 molecules upregulated in LA were negatively affected by PPARG (genes highlighted in red), and a total of six molecules suppressed in LA were stimulated by PPARG (genes highlighted in blue). The detailed information regarding the network presented in Figure 1 can be found in PPARG_LA⟶PPARG-LA Regulation Pathway with each network-related entry, including the type of the relationship, supporting references, and related sentences from the references where the relationship has been identified.
To note, five LA-unregulated genes presented increased expression levels in the mega-analysis using 13 LA datasets, including COL1A1, SPP1, CXCL14, MMP9, and CCNB1. The depression of these genes by PPARG could exert an anti-LA effect during its pathological development. On the other hand, four out of six LA-suppressed genes presented decreased expression levels in the mega-analysis, including three genes (CAV1, PTEN, and FAS) and one microRNA (MIR145). The activation of these molecules could be other pathways where PPARG inhibits the progress of LA.
Notably, PPARG presented a decreased expression level ( value = 1.07e-05, ). This was consistent with previous studies . Heterogeneity analysis showed that there was no significant between-study variance of PPARG expression levels (). Therefore, the fixed-effects model was used for the mega-analysis of PPARG. Moreover, MLR analysis suggested that none of the three parameters (sample size, country of origin, and study date) was a significant influence factor for the PPARG expression changes ( value > 0.070). For detailed info of the mega-analysis results of PPARG and other molecules, please refer to PPARG_LA⟶Mega-analysis.
3.2. PPARG-LA Regulation Pathway
We also identified a regulatory pathway through which PPARG may bar the pathological development of LA, as shown in Figure 2. According to literature reports, there were seven LA promoters (highlighted in red in Figure 2) deactivated by PPARG. Out of these molecular triggers, two genes presented increased expression levels in LA patients according to the mega-analysis results, including CCR7 and TLR2.
Literature data mining also revealed three LA inhibitors that could be activated by PPARG (see Figure 2; highlighted in blue). Only one of them, MIR145, presented a decreased expression level in the case of LA according to the results of mega-analysis, which was consistent with the negative-regulation relationship between LA and MIR145 that has been identified in literature data mining (see Figure 1). For more detailed information regarding the network presented in Figure 2, please refer to PPARG_LA⟶PPARG-LA Regulation Pathway.
3.3. PEA Results and PPI Network
To investigate the biological functions of the 26 genes within the PPARG-LA regulatory pathways (Figures 1 and 2), a pathway enrichment analysis was executed by using Pathway Studio. A total of 24 out of these 26 genes were shared among the top 10 most significantly enriched pathways ( value < 1.40e-14, for FDR), which are presented in Table 2. The full 115 pathways and GO terms enriched with value < 1.00e-6, which encompassed all 26 genes, were presented in PPARG_LA⟶PEA. Notably, a majority of the shared pathways highlighted by the PEA approach were related to cell proliferation and cell migration, which were implicated with LA.
Based on the significantly enriched pathways identified from PEA, a PPI network has been constructed, as shown in Figure 3 The number between two entities is the number of shared pathways out of the 115 pathways from PEA results. Please note that most of the molecules paly roles with all other molecules, suggesting the shared functionality of these molecules. Specifically, PPARG plays a role within 52 pathways and connecting with other molecules through pathways on average. This supports the literature-based relationship presented in the PPARG-LA regulatory pathways (Figures 1 and 2).
In this study, we aimed at exploring the possible genetic mechanisms of the linkage between PPARG and better survival in LA. First, we conducted a large-scale literature-based data mining to construct functional PPARG-LA regulatory pathways (Figures 1 and 2). The 25 genes within these literature-based pathways were identified as connected to both PPARG and LA with polarity. Then, PEA was conducted to study the pathological functions of these 25 genes. After that, an LA expression data-based mega-analysis was performed to explore the expression levels of these 25 genes in the case of LA. Our results showed that PPARG could exert influence on both the development and progression of LA, which may add new insights into the understanding of the PPARG-LA association. We provided detailed descriptions of the PPARG-LA regulatory pathways as follows.
There were 19 molecules that have been identified as common targets but were contradirectionally regulated by PPARG and LA. These molecules have been reported to present altered expression levels in case of LA, and part of them (eight proteins and one microRNA) have been confirmed from the mega-analysis using 13 LA datasets employed in this study (please see Figure 1, the entities with body color matched with the highlighting color). With both literature and experimental data support, the pathways constructed with these molecules have high potential to be the PPARG→LA regulating channel. Specifically, the expression levels of five LA markers (i.e., COL1A1, SPP1, CXCL14, MMP9, and CCNB1) were significantly unregulated in LA patients [6–10], which were confirmed in the mega-analysis in this study. PPARG has been shown to inhibit hepatic stellate cell proliferation and COL1A1 expression in vitro and in vivo . Similarly, the expressions of SPP1, CXCL14, MMP9, and CCNB1 have also been negatively regulated by PPARG [12–15]. On the other hand, PPARG could stimulate multiple molecules that inhibited by LA, including CAV1, PTEN, FAS, and MIR145. More descriptions of these PPARG regulations were provided in PPARG_LA→PPARG-LA Regulation Pathway. These pathways might partially explain the mechanisms of the development-blocking effect of PPARG on LA.
Moreover, PPARG has been shown to play roles in the prognostic pathway of LA (Figure 2). Both the literature data mining and the mega-analysis results support the suppression of four LA promoters (IL1B, PTGS2, END1, and TNF) and stimulation of one LA inhibitors, MIR145. As an LA tumor cell suppressor , MIR145 was found to be involved in multiple other types of cancers, including breast cancer , colon cancer , and acute myeloid leukemia . However, as shown in Figure 1, in the case of LA, the expression levels of MIR145 will be downregulated, which also got confirmed from the 13-LA-datasets mega-analysis (see PPARG_LA→Mega-analysis). PPARG could activate the expression of MIR145 by directly binding to a PPAR response element in its promoter at 1207/-1194 bp from the transcription start site , which could be an important mechanism underlying the LA inhibition effect of PPARG.
PEA results indicated that each of the 25 PPARG-driven molecular triggers has at least two functional pathways shared with PPARG (see Figure 3). Most of these pathways were linked to cell proliferation, cell migration, and motility, indicating that PPARG may influence the LA development and progression through the cell metabolism and motivation pathways, which are important in the etiology of LA . More interestingly, we see that PPARG was enriched in all top six pathways (see Table 2). These results support the functional association between PPARG and the 25 molecular triggers identified from the literature data mining (PPARG_LA→PPARG-LA Regulation Pathway).
To note, the expression of PPARG has been downregulated in LA, which is consistent with previous findings . Therefore, our results suggest that the activation of PPARG could be a valid therapeutic strategy for the treatment of LA.
One of the limitations of this study was that the PPARG-driven LA pathways built only explored the genes connecting PPARG and LA. There could be more “bridge items” (e.g., functional class and compounds) to reveal more mechanisms underlying the PPARG_LA relation.
This study confirmed the downregulation of PPARG in the case of LA and revealed multiple pathways through which PPARG could play roles as a LA blocker. Our results shed light on the understanding of the PPARG-LA association, suggesting PPARG as a valuable therapeutic target for the treatment of LA.
The data in our study are available from the corresponding author upon reasonable request. The supplementary data and results were available at http://www.gousinfo.com/database/Data_Genetic/PPARG_LA.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Min Zhao, Xiaoyang Li, and Yunxiang Zhang contributed equally to this work.
- J. Ansari, R. E. Shackelford, and H. El-Osta, “Epigenetics in non-small cell lung cancer: from basics to therapeutics,” Translational Lung Cancer Research, vol. 5, no. 2, pp. 155–171, 2016.
- N. Sawabata, E. Miyaoka, H. Asamura et al., “Japanese lung cancer registry study of 11,663 surgical cases in 2004: demographic and prognosis changes over decade,” Journal of Thoracic Oncology, vol. 6, no. 7, pp. 1229–1235, 2011.
- 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.
- R. Edgar, M. Domrachev, and A. E. Lash, “Gene Expression Omnibus: NCBI gene expression and hybridization array data repository,” Nucleic Acids Research, vol. 30, no. 1, pp. 207–210, 2002.
- Y. Zhang, B. Chen, Y. Wang et al., “NOTCH3 overexpression and posttranscriptional regulation by miR-150 were associated with EGFR-TKI resistance in lung adenocarcinoma,” Oncology Research Featuring Preclinical and Clinical Cancer Therapeutics, vol. 27, no. 7, pp. 751–761, 2019.
- H. M. Fathallah-Shaykh, M. Rigen, L. J. Zhao et al., “Mathematical modeling of noise and discovery of genetic expression classes in gliomas,” Oncogene, vol. 21, no. 47, article BF1205654, pp. 7164–7174, 2002.
- R. Shaykhiev, R. Sackrowitz, T. Fukui et al., “Smoking-induced CXCL14 expression in the human airway epithelium links chronic obstructive pulmonary disease to lung cancer,” American Journal of Respiratory Cell and Molecular Biology, vol. 49, no. 3, pp. 418–425, 2013.
- Y. L. Wen and L. Li, “Correlation between matrix metalloproteinase-9 and vascular endothelial growth factor expression in lung adenocarcinoma,” Genetics and Molecular Research, vol. 14, no. 4, pp. 19342–19348, 2013.
- T. Wang, X. B. Yan, J. J. Zhao et al., “Gene associated with retinoid-interferon-induced mortality-19 suppresses growth of lung adenocarcinoma tumor in vitro and in vivo,” Lung Cancer, vol. 72, no. 3, pp. 287–293, 2011.
- S. Zheng and A. Chen, “Disruption of transforming growth factor-beta signaling by curcumin induces gene expression of peroxisome proliferator-activated receptor-gamma in rat hepatic stellate cells,” American Journal of Physiology-Gastrointestinal and Liver Physiology, vol. 292, no. 1, pp. G113–G123, 2007.
- M. Shan, R. You, X. Yuan et al., “Agonistic induction of PPARγ reverses cigarette smoke-induced emphysema,” Journal of Clinical Investigation, vol. 124, no. 3, pp. 1371–1381, 2014.
- M. Takahashi, Y. Takahashi, K. Takahashi et al., “CXCL14 enhances insulin-dependent glucose uptake in adipocytes and is related to high-fat diet-induced obesity,” Biophysical Research Communications, vol. 364, no. 4, pp. 1037–1042, 2007.
- E. L. Schiffrin, “Peroxisome proliferator-activated receptors and cardiovascular remodeling,” American Journal of Physiology-Heart and Circulatory Physiology, vol. 288, no. 3, pp. H1037–H1043, 2005.
- J. Wan, Z. Xiao, S. Chao et al., “Pioglitazone modulates the proliferation and apoptosis of vascular smooth muscle cells via peroxisome proliferators-activated receptor-gamma,” Diabetology & Metabolic Syndrome, vol. 6, no. 1, p. 101, 2014.
- Y. F. Chang, K. H. Lim, Y. W. Chiang et al., “STAT3 induces G9a to exacerbate HER3 expression for the survival of epidermal growth factor receptor-tyrosine kinase inhibitors in lung cancers,” BMC Cancer, vol. 19, no. 1, article 6217, p. 959, 2019.
- M. Götte, C. Mohr, C. Y. Koo et al., “miR-145-dependent targeting of junctional adhesion molecule A and modulation of fascin expression are associated with reduced breast cancer cell motility and invasiveness,” Oncogene, vol. 29, no. 50, pp. 6569–6580, 2010.
- J. Zhang, H. Guo, H. Zhang et al., “Putative tumor suppressor miR-145 inhibits colon cancer cell growth by targeting oncogene friend leukemia virus integration 1 gene,” Cancer, vol. 117, no. 1, pp. 86–95, 2011.
- D. T. Starczynowski, R. Morin, A. McPherson et al., “Genome-wide identification of human microRNAs located in leukemia-associated genomic alterations,” Blood, vol. 117, no. 2, pp. 595–607, 2011.
- A. V. Das and R. M. Pillai, “Implications of miR cluster 143/145 as universal anti-oncomiRs and their dysregulation during tumorigenesis,” Cancer Cell International, vol. 15, no. 1, p. 92, 2015.
- Z. Li, J. Xia, M. Fang, and Y. Xu, “Epigenetic regulation of lung cancer cell proliferation and migration by the chromatin remodeling protein BRG1,” Oncogenesis, vol. 8, no. 11, p. 66, 2019.
Copyright © 2020 Min Zhao 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.