Next Article in Journal
Apoptotic MSCs and MSC-Derived Apoptotic Bodies as New Therapeutic Tools
Next Article in Special Issue
Review on Advanced Cancer Modeling for a Cancer Study
Previous Article in Journal
Personalized Selection of a CFTR Modulator for a Patient with a Complex Allele [L467F;F508del]
Previous Article in Special Issue
Chemokine Homeostasis in Healthy Volunteers and during Pancreatic and Colorectal Tumor Growth in Murine Models
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 44
Issue 11
10.3390/cimb44110350
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Immune Pathway and Gene Database (IMPAGT) Revealed the Immune Dysregulation Dynamics and Overactivation of the PI3K/Akt Pathway in Tumor Buddings of Cervical Cancer

by
Yeseul Choi
1,2,
Nora Jee-Young Park
3,4,5,
Tan Minh Le
1,2,
Eunmi Lee
1,2,
Donghyeon Lee
1,2,
Hong Duc Thi Nguyen
1,2,
Junghwan Cho
5,
Ji-Young Park
3,4,
Hyung Soo Han
1,2,5,6,* and
Gun Oh Chong
5,7,8,*
1
Department of Biomedical Science, Graduate School, Kyungpook National University, Daegu 41944, Korea
2
BK21 Four Program, School of Medicine, Kyungpook National University, Daegu 41944, Korea
3
Department of Pathology, School of Medicine, Kyungpook National University, Daegu 41944, Korea
4
Department of Pathology, Kyungpook National University Chilgok Hospital, Daegu 41404, Korea
5
Clinical Omics Institute, Kyungpook National University, Daegu 41405, Korea
6
Department of Physiology, School of Medicine, Kyungpook National University, Daegu 41944, Korea
7
Department of Obstetrics and Gynecology, School of Medicine, Kyungpook National University, Daegu 41944, Korea
8
Department of Obstetrics and Gynecology, Kyungpook National University Chilgok Hospital, Daegu 41404, Korea
*
Authors to whom correspondence should be addressed.
Curr. Issues Mol. Biol. 2022, 44(11), 5139-5152; https://doi.org/10.3390/cimb44110350
Submission received: 5 September 2022 / Revised: 17 October 2022 / Accepted: 21 October 2022 / Published: 23 October 2022
(This article belongs to the Special Issue Advances in Molecular Pathogenesis Regulation in Cancer)
Browse Figures
Review Reports Versions Notes

Abstract

:
Tumor budding (TB) is a small cluster of malignant cells at the invasive front of a tumor. Despite being an adverse prognosis marker, little research has been conducted on the tumor immune microenvironment of tumor buddings, especially in cervical cancer. Therefore, RNA sequencing was performed using 21 formalin-fixed, paraffin-embedded slides of cervical tissues, and differentially expressed genes (DEGs) were analyzed. Immune Pathway and Gene Database (IMPAGT) was generated for immune profiling. “Pathway in Cancer” was identified as the most enriched pathway for both up- and downregulated DEGs. Kyoto Encyclopedia of Genes and Genomes Mapper and Gene Ontology further revealed the activation of the PI3K/Akt signaling pathway. An IMPAGT analysis revealed immune dysregulation even at the tumor budding stage, especially in the PI3K/Akt/mTOR axis, with a high efficiency and integrity. These findings emphasized the clinical significance of tumor buddings and the necessity of blocking the overactivation of the PI3K/Akt/mTOR pathway to improve targeted therapy in cervical cancer.

1. Introduction

Tumor budding (TB) is the formation of a cluster of fewer than five malignant cells at the invasive edge of a tumor [1,2,3]. First established in colorectal cancer, tumor budding is believed to be related to the epithelial–mesenchymal transition (EMT), which involves the loss of cell–cell connectivity along with epithelial markers, such as E-cadherin, and the gain of mesenchymal markers, such as N-cadherin or vimentin [4,5,6,7]. After morphologically transitioning, tumor cells gain a locomotive ability, allowing them to migrate and invade secondary sites and eventually leading to metastasis. The newly developed disseminated tumor cells undergo another transition called the mesenchymal–epithelial transition, which recapitulates the pathology of primary epithelial tumors and promotes tumor formation at the attached site. Briefly, TB, which exhibits EMT traits, is associated with tumor initiation, invasion, metastasis, and overall cancer progression, and it is a clinically adverse prognosis marker [3,8,9,10,11,12,13].
Despite its clinical significance, the tumor immune microenvironment (TIME) of TB has not been extensively studied, especially in cervical cancer. To date, immune gene databases such as ImSig (Immune Cell Gene Signatures) [14] or platforms such as CIBERSORT [15] or TIMER [16] have been used for immune profiling. However, these platforms do not provide a list of genes involved in specific immune pathways. Moreover, CIBERSORT only presents relative cell expression, which makes it difficult to comprehend the degree of differential gene expression.
Therefore, we used a slightly different approach than traditional analysis methods. Instead of using diverse databases to infer biological implications and deduce potential clinically meaningful biomarkers in cancer etiology, we generated a new immune database, namely, Immune Pathway and Gene Database (IMPAGT), based on the gene set of CIBERSORT. Using IMPAGT, the tumor immune microenvironment of tumor buddings in cervical cancer was investigated.

2. Materials and Methods

2.1. Patient Selection

Following approval by the Institutional Review Board of Kyungpook National University Chilgok Hospital (approval no.: 2020-10-003), 21 patients with cervical cancer were selected for the study. Cervical cancer tissue from each patient was collected after radical hysterectomy with pelvic and/or para-aortic lymphadenectomy to treat early-stage and locally advanced cervical cancer. All patients were clinically staged according to the revised 2009 International Federation of Gynecologic Obstetrics staging system for cervical cancer [17].

2.2. Pathological Process and Tumor Budding Counts

All cervical tissue specimens were collected after biopsy, and formalin-fixed, paraffin-embedded (FFPE) slides were generated. The FFPE slides were stained with hematoxylin and eosin (H&E), and multiple sections of H&E-stained slides were examined by pathologists at Kyungpook National University Chilgok Hospital. Based on the criterion used in a previous study (high tumor budding, ≥3 tumor buddings), among 21 cervical cancer samples, 15 had high tumor budding, whereas the remaining samples had low tumor budding [18]. The medical records of all patients were reviewed, and their clinicopathologic parameters were collected. The patient information is described in Table 1.

2.3. RNA Sequencing

RNA was extracted from FFPE slides of 21 patients using the ReliaPrep™ FFPE Total RNA Miniprep kit (Promega, Madison, WI, USA). Extracted RNA was then synthesized into cDNA using the TruSeq RNA Exome Kit (Illumina, San Diego, CA, USA), and an RNA sequencing dataset was generated using the Nova 6000 platform (Illumina). The raw data were quantified using the STAR [19] and RSEM [20] quantification tools, and differentially expressed gene (DEG) analysis was performed using DESeq2 [21]. The threshold for identifying DEGs was adjusted to absolute log2 fold change ≥ 1 and p < 0.05. The DEGs were plotted in a volcano plot graph using the “ggplot” package in R.

2.4. CIBERSORT and LM22

CIBERSORT is a computational platform designed by the Alizadeh Laboratory at Stanford University (Palo Alto, CA, USA) that uses the machine learning technique, support vector regression, to deconvolute and characterize diverse cell types from bulk sequencing data [15]. The platform requires two input data files for analysis: (1) a vector consisting of gene expression profiles (GEPs) and (2) a “signature matrix” containing signature genes for the cell subset of interest.
For immune-related analysis, LM22 was used as a signature matrix. According to Alizadeh et al., LM22 contains 547 genes that distinguish 22 mature human hematopoietic populations isolated from peripheral blood or in vitro culture conditions, including seven T cell types, naïve and memory B cells, plasma cells, NK cells, and myeloid subsets. CIBERSORT analysis using this LM22 matrix provides proportional gene expression measures of the input GEP file corresponding to 22 distinct immune cell types [15].
The 547 gene symbols of LM22 are listed in Table S1.

2.5. Generation of IMPAGT

STRINGdb [22] is a database that provides protein–protein interaction (PPI) networks with functional enrichment analysis using diverse functional classification tools, such as Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) [23]. For our analysis, all 547 gene symbols of LM22 were input into the “Multiple proteins” tab of STRINGdb (version 11.5) using “Homo sapiens” as the corresponding organism. All outputted KEGG pathways were retrieved. Each pathway was individually searched using the KEGG pathway database webpage (https://www.genome.jp/kegg/pathway.html (accessed on 13 June 2022)), and all corresponding genes were collected using the R package “KEGGREST”. An LM22-based IMPAGT was generated, with the KEGG pathways of LM22 presented in columns and the corresponding gene symbols presented in rows.

2.6. Intersection of DEGs in Tumor Budding with IMPAGT

DEGs overlapping in tumor budding and IMPAGT were identified using the basic R function “intersect”. Then, the intersecting DEGs were further searched using IMPAGT, and upregulated (orange) and downregulated (blue) DEGs were highlighted. Subsequently, the number of highlighted DEGs in each pathway was counted.

2.7. Enriched Pathway Visualization

Pathways with more than 10 DEGs were considered significantly enriched in the TB environment of cervical cancer. Among these pathways, the most enriched pathway was visualized using several pathway visualization tools, such as KEGG Mapper Color [24], DAVID [25], and ShinyGo [26]. The final visualization result was obtained from KEGG Mapper.
From the KEGG Mapper Color website (https://www.genome.jp/kegg/mapper/color.html (accessed on 20 June 2022)), “hsa” was selected as the search mode to identify human-specific pathways. The lists of DEGs were uploaded, with upregulated DEGs colored red and downregulated DEGs colored cyan.

2.8. Functional Enrichment Analysis Using DAVID

GO analysis was performed using the DAVID database (https://david.ncifcrf.gov/tools.jsp (accessed on 20 July 2022)) [25,27]. In total, 133 upregulated DEGs that overlapped with IMPAGT genes were entered as a gene list of interest, and “OFFICIAL_GENE_SYMBOL” was selected as the identifier. Using “Homo Sapiens; 9606” as the species, the gene list was analyzed. In the GO analysis, only biological process (BP) terms were retrieved.

3. Results

3.1. DEGs in Samples with High or Low Tumor Budding

The RNA sequencing revealed 1253 DEGs between patients with high tumor budding and those with low tumor budding. Among these DEGs, 903 genes were upregulated and 349 were downregulated. A volcano plot of the top 30 up- and downregulated DEGs is presented in Figure 1.

3.2. Creation of IMPAGT

When the 547 genes of LM22 were analyzed in STRINGdb, the resultant PPI network was closely interrelated, forming an intricate ball cluster (Figure 2). There were 519 nodes and 6049 edges with significant PPI enrichment values of less than 1.0 × 10−16.
STRINGdb also identified the KEGG pathways associated with genes in the LM22 dataset. In total, 51 pathways from the KEGG database were analyzed, all of which were related to immune reactions of the human body. From these pathways, those related to direct viral infection, such as influenza A pathways, or those irrelevant to tumor budding, such as asthma pathways, were excluded, leaving 26 immune pathways. Meanwhile, 4131 genes were collected from the 26 pathways using KEGGREST. When overlapping genes were filtered, 2534 genes were left. This list of 2534 genes was considered the “IMPAGT gene set” which could be used for immune analysis. IMPAGT and the IMPAGT gene set are presented in Table S2.

3.3. Intersection of DEGs with the IMPAGT gene set

In total, 133 upregulated DEGs overlapped with the IMPAGT gene set, and they were engaged in all 26 immune pathways. Considering pathways with more than 10 DEGs as significantly enriched pathways, seven significantly enriched pathways were associated with upregulated DEGs: “Pathway in Cancer (Hsa 05200)”, “Antigen processing and presentation (Hsa 04612)”, “Transcriptional misregulation in Cancer (Hsa 05202)”, “MAPK signaling pathway (Hsa 04010)”, “NOD-like receptor signaling pathway (Hsa 04621)”, “Cell adhesion molecules (Hsa 04514)”, and “Cytokine-cytokine receptor interaction (Hsa 04060)”.
Among 31 downregulated DEGs, there were 31 overlapping DEGs, which were involved in 17 different immune pathways. The only significantly enriched pathway was “Pathway in Cancer”, and the remaining pathways contained fewer than five DEGs.
To summarize, “Pathway in Cancer” was the most strongly enriched pathway, indicating that it was the most strongly activated and dysregulated pathway in our sample. It was considered to be the most significantly enriched pathway in the TB environment of cervical cancer.
The significantly enriched pathways for both upregulated and downregulated DEGs, along with the number of DEGs involved in each pathway, are listed in Table 2.

3.4. Pathway Visualization Revealed Specific Enrichment of the PI3K/Akt Signaling Pathway

“Pathway in Cancer” was selected as the most significantly enriched pathway. Subsequently, it was visualized using KEGG Mapper Color. To clarify the overall involvement of TB DEGs in the pathway, all 1253 DEGs were inserted into KEGG Mapper.
The DEG-mapped “Pathway in Cancer” revealed the points at which the upregulated and downregulated DEGs engaged the pathway (Figure 3). Most DEGs were mainly involved in the “Signaling molecules and interaction” and “Signal transduction” metabolism pathways, such as the PI3K/Akt, JAK/STAT, MAPK, and Notch signaling pathways. Because most DEGs were involved in the PI3K/Akt pathway, this pathway was examined in more detail. All signal transduction pathways were directed to “Sustained Angiogenesis”, “Evading apoptosis”, and “Proliferation”.

3.5. GO Analysis of 133 Overlapping Upregulated DEGs

GO analysis on 133 upregulated DEGs that intersected with IMPAGT identified BP functions closely related to the PI3K/AKT pathway. When all significant BP terms were hierarchically ordered by percentage relative to the total input genes with p-values of <0.05, ”Signal transduction (GO: 0007165)” and “Immune response (GO: 0006955)” were the top two BP terms. Other significantly enriched BP terms included “Positive regulation of cell proliferation (GO: 0008284)”, “Positive regulation of apoptotic process (GO: 0043065)”, “Regulation of apoptotic process (GO: 0042981)”, “Regulation of cell proliferation (GO: 0042127)”, and “Positive regulation of protein phosphorylation (GO: 0001934)” (Table 3).
Meanwhile, the 133 overlapping DEGs were further visualized using STRINGdb, and the PPI network revealed a tendency to cluster into three groups, which were termed “Apoptotic process”, “Protein Phosphorylation”, and ”Cell proliferation” (Figure 4).

4. Discussion

This study aimed to investigate the tumor immune microenvironment (TIME) of TB in cervical cancer using the newly generated immune database IMPAGT. For immune profiling analyses, researchers mainly use immune gene databases, such as ImSig [14] or ImmuCellDB [28]. The other popular approach is to use bioinformatics applications, such as CIBERSORT [15], a deconvolution tool that characterizes immune cell compositions. However, all of these methods have limitations. First, the usage of several different databases for data interpretation is extremely time-consuming, and this is associated with a high risk of bias. Second, although CIBERSORT allows users to estimate the immune cell composition, it only provides proportional values of cell composition that are summed in one, thereby making it difficult to identify the differential expression of the input genes. Moreover, it does not provide lists of genes involved in specific immune pathways.
Considering these limitations, our study holds significance in terms of several aspects. First, IMPAGT allows researchers to interpret their data with high efficiency and integrity and a low risk of bias. IMPAGT, the “two in one” immune database, provides immune pathways and the corresponding genes of each pathway simultaneously. The integrated database enables users to investigate their genes of interest in a more rapid and efficient manner than individually searching and merging separate results from multiple databases. Furthermore, instead of obtaining “one gene–one pathway” associations based on the researcher’s interpretation, IMPAGT provides several immune pathways in which certain genes of interest are involved. In particular, it proposes “one gene–many immune pathways” interconnections. As it offers diverse immune pathways in which specific genes are engaged, IPMAGT encourages researchers to view and interpret the tumor immune microenvironment (TIME) from a more holistic, comprehensive, and integrative perspective.
Second, IMPAGT presents the dysregulation of the immune system even at the tumor budding stage in cervical cancer. “Pathway in Cancer” was the most strongly enriched and activated pathway for both up- and downregulated DEGs. This highlighted the presence of abnormal shifts in gene expression in TB, indicating the activation of carcinogenesis dynamics even from a small cluster of tumor cells. This process supported the initiation of cancer from tumor budding and further revalidated their clinical significance.
In addition, DEG-mapped pathway visualization of “Pathway in Cancer” using KEGG Mapper revealed overactivation of “Signal transduction pathways”, such as the PI3K/Akt, JAK/STAT, and Notch signaling pathways. These signaling pathways eventually pointed toward “Sustained angiogenesis”, “Evading apoptosis”, and “Proliferation”. A detailed investigation of DEG-mapped pathways revealed that up- and downregulated DEGs were mainly involved in the PI3K/Akt signaling pathway. Furthermore, the GO analysis of 133 IMPAGT-overlapped upregulated DEGs revealed that the immune-related genes were directly associated with the characteristics of the PI3K/Akt pathway.
The PI3K/Akt/mTOR signaling pathway is one of the most common aberrantly stimulated pathways in solid tumors [29]. The pathway is critical for processes involved in cancer progression, including cell proliferation, metabolism, apoptosis, and metastasis [30,31,32]. Moreover, this pathway appears to be important for tumor initiation [30,33]. PI3K is a intracellular lipid kinase that is usually activated by receptor tyrosine kinase (RTK) in response to various stimuli, such as growth factors, cytokines, hormones, and insulin [34]. Among three classes of PI3Ks, Class 1A PI3K is the most widely studied in relation to cancer [33]. Class 1A PI3Ks are heterodimers consisting of a catalytic subunit (p110α, p110β, or p110γ) and a regulatory subunit (p85α, p85β, p50α, p55α, or p55γ) [33,34,35,36]. When RTK is activated, regulatory p85 releases p110 catalytic subunits and generates PI3K, which phosphorylates PIP2 to PIP3 [37]. The accumulation of PIP3 on the cell membrane leads to the activation of downstream Akt upon the phosphorylation of Thr308 by PDK1 and Ser473 by mTORC2 [34,38,39]. Subsequently, fully activated Akt phosphorylates multiple downstream target proteins and substrates, such as mTORC1, p21, GSK3, and FOXO, and plays a dominant role in the signal transduction of the entire PI3K signaling pathway [33,37].
The oncogenic stimulation of this pathway is caused by several mechanisms: inappropriate activation of RTK; somatic mutations in PI3KCA, which encodes p110α of PI3K; and the loss of the tumor suppressor PTEN and other molecular signals [33,36,40,41]. PIK3CA is the second most frequently mutated oncogene across cancers, comprising 39% of mutations found in cervical cancer [42,43]. It also has been found that PIK3CA-mutant cervical cancer exhibits poor prognostic features, such as a higher tumor mutation burden and worse survival [31,32,44,45]. PI3K/Akt/mTOR dysregulation in cancer caused by the aforementioned reasons is associated with several cancer hallmarks, such as cell cycle, proliferation, autophagy, apoptosis, angiogenesis and metabolism, metastasis, and the EMT [39,46,47]. As these hallmarks were depicted in our enriched pathway visualization, each point related to the PI3K/Akt/mTOR pathway was examined in more detail.
The PI3K/Akt/mTOR axis is related to antiapoptosis and proliferation in cancer cells. Apoptosis is a programmed cell death pathway that is activated by both the extrinsic death receptor and intrinsic mitochondrial pathways. Peng et al. argued that PI3K/Akt/mTOR signaling leads to the binding of extrinsic death complexes, such as FasL and Fas, which activate caspase 8 and trigger apoptosis with other caspases [46]. Additionally, they mentioned that the balance between proapoptotic (Bad, Bax) and antiapoptotic proteins (Bcl-2, Bcl-x) is lost, thereby promoting the evasion of apoptosis [38,45,46]. Moreover, Akt can negatively regulate transcription factors, such as NF-κB or FOXO, and increase the production of antiapoptotic and prosurvival genes [43,45,48]. In addition to the inhibition of apoptosis, Akt also regulates cell survival through p27, GSK3, or the mTOR complex. To elaborate, Akt regulates the cell cycle via p27 and produces proteins involved in the G1/S phase of the cell cycle, such as cyclin D1, through GSK3 [46,49]. The overactivation of mTORC2 leads to the hyperaccumulation of Akt, which stimulates mTORC1 to phosphorylate translation effectors, such as S6K or 4EBP1, to enhance protein synthesis for cell survival [39,45,46,50].
The PI3K/Akt/mTOR pathway is also critical in sustaining angiogenesis by serving as a master controller of angiogenic signaling in the endothelium [33]. Tumor angiogenesis is caused by an imbalance between proangiogenic factors such as VEGF and antiangiogenic factors such as angiostatin or endostatin [51]. Peng et al. insisted that hypoxia stabilizes HIF-1a and stimulates the PI3K/Akt/mTOR axis in malignant cells [46]. Such activation of the axis boosts VEGF production and enhances vascular permeability, resulting in leaky vessels, sluggish blood flow, and high intestinal pressure, which eventually leads to sustained angiogenesis [46].
Although it was not specifically identified in the DEG-mapped pathway, PI3K/Akt/mTOR signaling is also closely associated with the EMT [52,53]. As mentioned above, epithelial cells change their morphology into a more mesenchymal–fibroblast phenotype. The morphological change is accompanied by decreased intercellular cell–cell adhesion and enhanced motility. Such alteration implies a decreased expression of the epithelial marker E-cadherin and an increased expression of mesenchymal markers, such as vimentin and fibronectin, which are regulated by Snail or Twist signals [10]. Qureshi et al. reported that PI3K signaling is stimulated by EGFR along with hypoxia, and it affects transcription factors, such as Snail and Twist [11]. In addition, NF-κB, mTORC2, and TGF-β are linked to the EMT via the activation of the PI3K/Akt/mTOR signaling pathway [54,55,56,57,58,59,60]. Owing to its importance, several clinical trials have examined the PI3/Akt/mTOR pathway. Several types of inhibitors have been utilized, including pan-PI3K inhibitors targeting all four isoforms of Class 1 PI3K, isoform-selective inhibitors targeting a single isoform of Class 1 PI3K, dual PI3K/mTOR inhibitors, and specific target inhibitors [36,41,42,45]. These small-molecule inhibitors support the inappropriate activation of the PI3K/Akt/mTOR signaling pathway in malignancies in cervical cancer. For example, treatment with rapamycin (mTOR inhibitor) blocks angiogenesis and enhances autophagy in cervical cancer cell lines [46,61]. The PI3K inhibitor LY294002 increases apoptosis and radiosensitization in cervical cancer cell lines [62,63].
Despite the necessity and continuous trials of PI3K inhibitors, only a few treatments have been approved for clinical use for several reasons [64]. First, sufficient doses needed to generate beneficial effects cannot be given to patients because of the risk of “on-target side effects”. Diabetes can be induced by high doses of PI3K inhibitors because the pathway is involved in insulin signaling. In addition, as PI3K/Akt/mTOR is engaged with many effectors, unexpected “off-target side effects” could develop. Furthermore, a recent study proposed a noncanonical PI3K pathway in the nuclei of cells, which is different from the canonical PI3K pathway located in the plasma membrane [65]. Such findings highlight the importance of selective pharmacological inhibition specifically targeting molecules of the PI3K/Akt/mTOR axis without affecting other pharmacodynamics markers or causing unintended potential toxicity [43,66].
On that note, a recently developed technology, termed proteolysis-targeting chimeras (PROTACs), can be used to solve issues related to small-molecule inhibitors. PROTACs are bivalent molecules consisting of a ligand of the protein of interest (POI) and a covalently linked ligand of an E3 ubiquitin ligase. PROTACs recruit E3 ligase in proximity. After binding and forming the complex, the POI is transported to the 26S proteasome for degradation [67,68]. PROTACS repeatedly search for the targeted protein, which facilitates POI elimination at extremely low doses [69,70]. As PROTACs specifically eliminate POIs with better selectivity, low toxicity, and fewer on-/off-target side effects, they may be adapted as clinically useful treatments [71]. Additional studies of PROTACs targeting PI3K/Akt/mTOR signaling pathways are being conducted [70,72].
The limitations of our study include the small number of samples with pathway-level results. Because the current data focused on the most enriched pathway related to the TIME of TB in cervical cancer, specific details, such as gene-level biomarkers or specific aspects of molecular reaction in the PI3K/Akt/mTOR pathway, have not been elucidated. However, these issues can be resolved by planned follow-up research on PI3K PROTACs with more clinical samples and validation studies in cervical cancer cell lines. Moreover, protein-level analyses, such as Immunohistochemistry and Western blotting, are planned in the next follow-up study as well. Despite these limitations, our study validated the presence of immune dysregulation in TB of cervical cancer, especially the dysregulation of the PI3K/Akt/mTOR signaling pathway. In addition, the findings emphasized the clinical significance of tumor budding in relation to tumorigenesis. Furthermore, the study asserted the necessity of in-depth investigations on the PI3K/Akt/mTOR axis to develop more specifically targeted patient-tailored treatments in cervical cancer.

5. Conclusions

The newly generated database IMPAGT revealed the existence of immune dysregulation dynamics and the specific enrichment of the PI3K/Akt/mTOR signaling pathway in the tumor immune microenvironment of cervical cancer with high tumor budding. Further research focused on the PI3K/Akt/mTOR axis and selective therapeutic approaches targeting this pathway could help develop better targeted therapy strategies for cervical cancer.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb44110350/s1, Table S1: List of 547 gene symbols of LM22; Table S2: IMPAGT and IMPAGT Gene Set.

Author Contributions

Conceptualization, Y.C., N.J.-Y.P., H.S.H. and G.O.C.; Formal analysis, Y.C., N.J.-Y.P. and H.S.H.; Investigation, N.J.-Y.P., J.-Y.P., H.S.H. and G.O.C.; Methodology, Y.C., T.M.L., E.L., D.L. and H.D.T.N.; Project administration, N.J.-Y.P., H.S.H. and G.O.C.; Writing—original draft preparation, Y.C.; Writing—review and editing, T.M.L., E.L., D.L., H.D.T.N., J.-Y.P. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (Nos. NRF-2020R1F1A1076816 and NRF-2022R1F1A1061498) and a grant of the Brain Pool Program through the NRF funded by the Ministry of Science and ICT (No. NRF-2020H1D3A2A02102040).

Institutional Review Board Statement

This study was approved by the Institutional Review Board of Kyungpook National University Chilgok Hospital (KNUCH 2020-10-003).

Informed Consent Statement

Patient consent was waived due to the retrospective study.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lugli, A.; Kirsch, R.; Ajioka, Y.; Bosman, F.; Cathomas, G.; Dawson, H.; El Zimaity, H.; Fléjou, J.-F.; Hansen, T.P.; Hartmann, A.; et al. Recommendations for reporting tumor budding in colorectal cancer based on the International Tumor Budding Consensus Conference (ITBCC) 2016. Mod. Pathol. 2017, 30, 1299–1311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Prall, F. Tumour budding in colorectal carcinoma. Histopathology 2007, 50, 151–162. [Google Scholar] [CrossRef]
  3. Li, H.; Xu, F.; Li, S.; Zhong, A.; Meng, X.; Lai, M. The tumor microenvironment: An irreplaceable element of tumor budding and epithelial-mesenchymal transition-mediated cancer metastasis. Cell Adh. Migr. 2016, 10, 434–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Schmalhofer, O.; Brabletz, S.; Brabletz, T. E-cadherin, β-catenin, and ZEB1 in malignant progression of cancer. Cancer Metastasis Rev. 2009, 28, 151–166. [Google Scholar] [CrossRef]
  5. Grigore, A.D.; Jolly, M.K.; Jia, D.; Farach-Carson, M.C.; Levine, H. Tumor budding: The name is EMT. partial EMT. J. Clin. Med. 2016, 5, 51. [Google Scholar] [CrossRef]
  6. Koelzer, V.H.; Zlobec, I.; Lugli, A. Tumor budding in colorectal cancer-Ready for diagnostic practice? Hum. Pathol. 2016, 47, 4–19. [Google Scholar] [CrossRef] [PubMed]
  7. van Wyk, H.; Park, J.; Roxburgh, C.; Horgan, P.; Foulis, A.; McMillan, D.C. The role of tumour budding in predicting survival in patients with primary operable colorectal cancer: A systematic review. Cancer Treat. Rev. 2015, 41, 151–159. [Google Scholar] [CrossRef] [PubMed]
  8. Yao, D.; Dai, C.; Peng, S. Mechanism of the mesenchymal-epithelial transition and its relationship with metastatic tumor formation. Mol. Cancer Res. 2011, 9, 1608–1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Pastushenko, I.; Blanpain, C. EMT transition states during tumor progression and metastasis. Trends Cell Biol. 2019, 29, 212–226. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, M.Y.; Chou, C.Y.; Tang, M.J.; Shen, M.R. Epithelial-mesenchymal transition in cervical cancer: Correlation with tumor progression, epidermal growth factor receptor overexpression, and snail up-regulation. Clin. Cancer Res. 2008, 14, 4743–4750. [Google Scholar] [CrossRef]
  11. Qureshi, R.; Arora, H.; Rizvi, M.A. EMT in cervical cancer: Its role in tumour progression and response to therapy. Cancer Lett. 2015, 356, 321–331. [Google Scholar] [CrossRef]
  12. Zhang, X.; Liu, S.; Zhu, Y. A-kinase-interacting protein 1 promotes EMT and metastasis via PI3K/Akt/IKKβ pathway in cervical cancer. Cell Biochem. Funct. 2020, 38, 782–791. [Google Scholar] [CrossRef]
  13. Dawson, H.; Lugli, A. Molecular and pathogenetic aspects of tumor budding in colorectal cancer. Front. Med. 2015, 2, 11. [Google Scholar] [CrossRef] [Green Version]
  14. Nirmal, A.J.; Regan, T.; Shih, B.B.; Hume, D.A.; Sims, A.H.; Freeman, T.C. Immune cell gene signatures for profiling the microenvironment of solid tumors. Cancer Immunol. Res. 2018, 6, 1388–1400. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, B.; Khodadoust, M.S.; Liu, C.L.; Newman, A.M.; Alizadeh, A.A. Profiling tumor infiltrating immune cells with CIBERSORT. Methods Mol. Biol. 2018, 1711, 243–259. [Google Scholar]
  16. Li, T.; Fan, J.; Wang, B.; Traugh, N.; Chen, Q.; Liu, J.S.; Li, B.; Liu, X.S. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017, 77, e108–e110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Pecorelli, S.; Zigliani, L.; Odicino, F. Revised FIGO staging for carcinoma of the cervix. Int. J. Gynecol. Obstet. 2009, 105, 107–108. [Google Scholar] [CrossRef] [PubMed]
  18. Chong, G.O.; Park, N.J.-Y.; Han, H.S.; Cho, J.; Kim, M.-G.; Choi, Y.; Yeo, J.Y.; Lee, Y.H.; Hong, D.G.; Park, J.Y. Intratumoral budding: A novel prognostic biomarker for tumor recurrence and a potential predictor of nodal metastasis in uterine cervical cancer. Eur. J. Surg. Oncol. 2021, 47, 3182–3187. [Google Scholar] [CrossRef]
  19. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  20. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 1–21. [Google Scholar] [CrossRef] [PubMed]
  22. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar] [CrossRef] [PubMed]
  23. Kanehisa Laboratories. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999, 27, 29–34. Available online: http://www.kegg.jp/ (accessed on 26 July 2022). [CrossRef] [PubMed] [Green Version]
  24. Kanehisa, M.; Sato, Y.; Kawashima, M. KEGG mapping tools for uncovering hidden features in biological data. Protein Sci. 2022, 31, 47–53. [Google Scholar] [CrossRef]
  25. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, 216–221. [Google Scholar] [CrossRef]
  26. Ge, S.X.; Jung, D. ShinyGO: A graphical enrichment tool for ani-mals and plants. bioRxiv 2018, 36, 315150. [Google Scholar]
  27. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  28. Chen, Z.; Na, H.; Wu, A. ImmuCellDB: An indicative database of immune cell composition from different tissues and disease conditions in mouse and human. Front. Immunol. 2021, 12, 3226. [Google Scholar] [CrossRef]
  29. Liu, P.; Cheng, H.; Roberts, T.M.; Zhao, J.J. Targeting the phosphoinositide 3-kinase pathway in cancer. Nat. Rev. Drug Discov. 2009, 8, 627–644. [Google Scholar] [CrossRef] [Green Version]
  30. Courtney, K.D.; Corcoran, R.B.; Engelman, J.A. The PI3K pathway as drug target in human cancer. J. Clin. Oncol. 2010, 28, 1075–1083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. McIntyre, J.B.; Wu, J.S.; Craighead, P.S.; Phan, T.; Köbel, M.; Lees-Miller, S.P.; Ghatage, P.; Magliocco, A.M.; Doll, C.M. PIK3CA mutational status and overall survival in patients with cervical cancer treated with radical chemoradiotherapy. Gynecol. Oncol. 2013, 128, 409–414. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, L.; Wu, J.; Ling, M.T.; Zhao, L.; Zhao, K.N. The role of the PI3K/Akt/mTOR signalling pathway in human cancers induced by infection with human papillomaviruses. Mol. Cancer 2015, 14, 87. [Google Scholar] [CrossRef] [PubMed]
  33. Yuan, T.L.; Cantley, L.C. PI3K pathway alterations in cancer: Variations on a theme. Oncogene 2008, 27, 5497–5510. [Google Scholar] [CrossRef] [Green Version]
  34. Shi, X.; Wang, J.; Lei, Y.; Cong, C.; Tan, D.; Zhou, X. Research progress on the PI3K/AKT signaling pathway in gynecological cancer (Review). Mol. Med. Rep. 2019, 19, 4529–4535. [Google Scholar] [CrossRef] [Green Version]
  35. Engelman, J.A.; Luo, J.; Cantley, L.C. The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat. Rev. Genet. 2006, 7, 606–619. [Google Scholar] [CrossRef] [PubMed]
  36. Engelman, J.A. Targeting PI3K signalling in cancer: Opportunities, challenges and limitations. Nat. Rev. Cancer 2009, 9, 550–562. [Google Scholar] [CrossRef]
  37. Graupera, M.; Potente, M. Regulation of angiogenesis by PI3K signaling networks. Exp. Cell Res. 2013, 319, 1348–1355. [Google Scholar] [CrossRef]
  38. Chang, F.; Lee, J.T.; Navolanic, P.M.; Steelman, L.S.; Shelton, J.G.; Blalock, W.L.; A Franklin, R.; A McCubrey, J. Involvement of PI3K/Akt pathway in cell cycle progression, apoptosis, and neoplastic transformation: A target for cancer chemotherapy. Leukemia 2003, 17, 590–603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Liu, M.; Li, C.-M.; Chen, Z.-F.; Ji, R.; Guo, Q.-H.; Li, Q.; Zhang, H.-L.; Zhou, Y.-N. Celecoxib regulates apoptosis and autophagy via the PI3K/Akt signaling pathway in SGC-7901 gastric cancer cells. Int. J. Mol. Med. 2014, 33, 1451–1458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Wong, K.K.; Engelman, J.A.; Cantley, L.C. Targeting the PI3K signaling pathway in cancer. Curr. Opin. Genet. Dev. 2010, 20, 87–90. [Google Scholar] [CrossRef] [Green Version]
  41. Noorolyai, S.; Shajari, N.; Baghbani, E.; Sadreddini, S.; Baradaran, B. The relation between PI3K/AKT signalling pathway and cancer. Gene 2019, 698, 120–128. [Google Scholar] [CrossRef] [PubMed]
  42. Jiang, N.; Dai, Q.; Su, X.; Fu, J.; Feng, X.; Peng, J. Role of PI3K/AKT pathway in cancer: The framework of malignant behavior. Mol. Biol. Rep. 2020, 47, 4587–4629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Samuels, Y.; Ericson, K. Oncogenic PI3K and its role in cancer. Curr. Opin. Oncol. 2006, 18, 77–82. [Google Scholar] [CrossRef] [PubMed]
  44. Voutsadakis, I.A. PI3KCA mutations in uterine cervix carcinoma. J. Clin. Med. 2021, 10, 220. [Google Scholar] [CrossRef]
  45. Martini, M.; De Santis, M.C.; Braccini, L.; Gulluni, F.; Hirsch, E. PI3K/AKT signaling pathway and cancer: An updated review. Ann. Med. 2014, 46, 372–383. [Google Scholar] [CrossRef]
  46. Peng, Y.; Wang, Y.; Zhou, C.; Mei, W.; Zeng, C. PI3K/Akt/mTOR pathway and its role in cancer therapeutics: Are we making headway? Front. Oncol. 2022, 12, 819128. [Google Scholar] [CrossRef]
  47. Ye, B.; Jiang, L.L.; Xu, H.T.; Zhou, D.W.; Li, Z.S. Expression of pi3k/akt pathway in gastric cancerand its blockade suppresses tumor growth and metastasis. Int. J. Immunopathol. Pharmacol. 2012, 25, 627–636. [Google Scholar] [CrossRef] [Green Version]
  48. Hoxhaj, G.; Manning, B.D. The PI3K–AKT network at the interface of oncogenic signalling and cancer metabolism. Nat. Rev. Cancer 2020, 20, 74–88. [Google Scholar] [CrossRef]
  49. Jeong, S.-J.; Dasgupta, A.; Jung, K.-J.; Um, J.-H.; Burke, A.; Park, H.U.; Brady, J.N. PI3K/AKT inhibition induces caspase-dependent apoptosis in HTLV-1-transformed cells. Virology 2008, 370, 264–272. [Google Scholar] [CrossRef] [Green Version]
  50. Shiojima, I.; Walsh, K. Role of Akt signaling in vascular homeostasis and angiogenesis. Circ. Res. 2002, 90, 1243–1250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Karar, J.; Maity, A. PI3K/AKT/mTOR Pathway in Angiogenesis. Front. Mol. Neurosci. 2011, 4, 51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Herrerias, M.M.; Scott Budinger, G.R. Revisiting mTOR and Epithelial—Mesenchymal transition. Am. J. Respir. Cell Mol. Biol. 2020, 62, 669–670. [Google Scholar] [CrossRef]
  53. Chen, H.-T.; Liu, H.; Mao, M.-J.; Tan, Y.; Mo, X.-Q.; Meng, X.-J.; Cao, M.-T.; Zhong, C.-Y.; Liu, Y.; Shan, H.; et al. Crosstalk between autophagy and epithelial-mesenchymal transition and its application in cancer therapy. Mol. Cancer 2019, 18, 1–19. [Google Scholar] [CrossRef] [Green Version]
  54. Hao, Y.; Baker, D.; Dijke, P.T. TGF-β-mediated epithelial-mesenchymal transition and cancer metastasis. Int. J. Mol. Sci. 2019, 20, 2767. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lamouille, S.; Connolly, E.; Smyth, J.W.; Akhurst, R.J.; Derynck, R. TGf-β-induced activation of mTOR complex 2 drives epithelial-mesenchymal transition and cell invasion. Development 2012, 139, 1259–1273. [Google Scholar] [CrossRef]
  56. Roshan, M.K.; Soltani, A.; Soleimani, A.; Kahkhaie, K.R.; Afshari, A.R.; Soukhtanloo, M. Role of AKT and mTOR signaling pathways in the induction of epithelial-mesenchymal transition (EMT) process. Biochimie 2019, 165, 229–234. [Google Scholar] [CrossRef] [PubMed]
  57. Choi, D.; Kim, C.L.; Kim, J.E.; Mo, J.S.; Jeong, H.S. Hesperetin inhibit EMT in TGF-β treated podocyte by regulation of mTOR pathway. Biochem. Biophys. Res. Commun. 2020, 528, 154–159. [Google Scholar] [CrossRef]
  58. Zheng, S.; Luo, J.; Xie, S.; Lu, S.; Liu, Q.; Xiao, H.; Luo, W.; Huang, W.; Liu, K. Tumor budding of cervical squamous cell carcinoma: Epithelial-mesenchymal transition-like cancer stem cells? PeerJ 2022, 10, e13745. [Google Scholar] [CrossRef]
  59. Xu, W.; Yang, Z.; Lu, N. A new role for the PI3K/Akt signaling pathway in the epithelial-mesenchymal transition. Cell Adh. Migr. 2015, 9, 317–324. [Google Scholar] [CrossRef] [Green Version]
  60. Zuo, J.-H.; Zhu, W.; Li, M.-Y.; Li, X.-H.; Yi, H.; Zeng, G.-Q.; Wan, X.-X.; He, Q.-Y.; Li, J.-H.; Qu, J.-Q.; et al. Activation of EGFR promotes squamous carcinoma SCC10A cell migration and invasion via inducing EMT-like phenotype change and MMP-9-mediated degradation of E-cadherin. J. Cell Biochem. 2011, 112, 2508–2517. [Google Scholar] [CrossRef]
  61. Lu, R.; Yang, Z.; Xu, G.; Yu, S. miR-338 modulates proliferation and autophagy by PI3K/AKT/mTOR signaling pathway in cervical cancer. Biomed. Pharmacother. 2018, 105, 633–644. [Google Scholar] [CrossRef]
  62. Lee, C.M.; Fuhrman, C.B.; Planelles, V.; Peltier, M.R.; Gaffney, D.K.; Soisson, A.P.; Dodson, M.K.; Tolley, H.D.; Green, C.L.; Zempolich, K.A. Phosphatidylinositol 3-kinase inhibition by LY294002 radiosensitizes human cervical cancer cell lines. Clin. Cancer Res. 2006, 12, 250–256. [Google Scholar] [CrossRef] [PubMed]
  63. Liu, Y.; Cui, B.; Qiao, Y.; Zhang, Y.; Tian, Y.; Jiang, J.; Ma, D.; Kong, B. Phosphoinositide-3-Kinase inhibition enhances radiosensitization of cervical cancer in vivo. Int. J. Gynecol. Cancer 2011, 21, 100–105. [Google Scholar] [CrossRef]
  64. Janku, F.; Yap, T.A.; Meric-Bernstam, F. Targeting the PI3K pathway in cancer: Are we making headway? Nat. Rev. Clin. Oncol. 2018, 15, 273–291. [Google Scholar] [CrossRef]
  65. Chen, M.; Choi, S.; Wen, T.; Chen, C.; Thapa, N.; Lee, J.H.; Cryns, V.L.; Anderson, R.A. A p53-phosphoinositide signalosome regulates nuclear Akt activation. Nat. Cell Biol. 2022, 24, 1099–1113. [Google Scholar] [CrossRef]
  66. Dituri, F.; Mazzocca, A.; Giannelli, G.; Antonaci, S. PI3K functions in cancer progression, anticancer immunity and immune evasion by tumors. Clin. Dev. Immunol. 2011, 2011, 1–10. [Google Scholar] [CrossRef] [Green Version]
  67. Li, X.; Song, Y. Proteolysis-targeting chimera (PROTAC) for targeted protein degradation and cancer therapy. J. Hematol. Oncol. 2020, 13, 1–14. [Google Scholar] [CrossRef] [PubMed]
  68. Yu, F.; Cai, M.; Shao, L.; Zhang, J. Targeting protein kinases degradation by PROTACs. Front. Chem 2021, 9, 679120. [Google Scholar] [CrossRef] [PubMed]
  69. Scheepstra, M.; Hekking, K.; van Hijfte, L.; Folmer, R.H. Bivalent ligands for protein degradation in drug discovery. Comput. Struct. Biotechnol. J. 2019, 17, 160–176. [Google Scholar] [CrossRef]
  70. Sun, X.; Gao, H.; Yang, Y.; He, M.; Wu, Y.; Song, Y.; Tong, Y.; Rao, Y. Protacs: Great opportunities for academia and industry. Signal Transduct. Target. Ther. 2019, 4, 64. [Google Scholar] [CrossRef] [Green Version]
  71. Li, W.; Gao, C.; Zhao, L.; Yuan, Z.; Chen, Y.; Jiang, Y. Phthalimide conjugations for the degradation of oncogenic PI3K. Eur. J. Med. Chem. 2018, 151, 237–247. [Google Scholar] [CrossRef] [PubMed]
  72. Wang, H.; Li, C.; Liu, X.; Ma, M. Design, synthesis and activity study of a novel PI3K degradation by hijacking VHL E3 ubiquitin ligase. Bioorg. Med. Chem. 2022, 61, 116707. [Google Scholar] [CrossRef] [PubMed]
Cimb 44 00350 g001 550
Figure 1. A volcano plot of differentially expressed genes (DEGs) of high tumor budding in cervical cancer. Top 30 up-regulated and top 30 downregulated DEGs are plotted. The upper right corner indicates over-expressed DEGs and upper left corner indicates down-expressed DEGs.
Figure 1. A volcano plot of differentially expressed genes (DEGs) of high tumor budding in cervical cancer. Top 30 up-regulated and top 30 downregulated DEGs are plotted. The upper right corner indicates over-expressed DEGs and upper left corner indicates down-expressed DEGs.
Cimb 44 00350 g001
Cimb 44 00350 g002 550
Figure 2. Protein–protein interaction (PPI) network of LM22 on STRINGdb. STRINGdb analysis of LM22 presented a cluster that is very interrelated, an almost ball-shaped aggregate. In total, 51 pathways were identified which were all related to immune reactions of the human body.
Figure 2. Protein–protein interaction (PPI) network of LM22 on STRINGdb. STRINGdb analysis of LM22 presented a cluster that is very interrelated, an almost ball-shaped aggregate. In total, 51 pathways were identified which were all related to immune reactions of the human body.
Cimb 44 00350 g002
Cimb 44 00350 g003 550
Figure 3. Enriched pathway visualization with KEGG Mapper Color. DEG-mapped pathway visualization presented that DEGs were mainly involved in “Signal transduction metabolism”, such as the PI3K/Akt, JAK/STAT, and Notch signaling pathways, and were specifically enriched in “PI3K/Akt signaling pathway”. All pathways were directed towards “Sustained Angiogenesis”, “Evading Apoptosis”, and “Proliferation”.
Figure 3. Enriched pathway visualization with KEGG Mapper Color. DEG-mapped pathway visualization presented that DEGs were mainly involved in “Signal transduction metabolism”, such as the PI3K/Akt, JAK/STAT, and Notch signaling pathways, and were specifically enriched in “PI3K/Akt signaling pathway”. All pathways were directed towards “Sustained Angiogenesis”, “Evading Apoptosis”, and “Proliferation”.
Cimb 44 00350 g003
Cimb 44 00350 g004 550
Figure 4. BP clusters of 133 upregulated DEGs overlapped with IMPAGT. When PPI network was clustered with BP traits, it presented a tendency to cluster into three groups: “Apoptotic process”, “Protein phosphorylation”, and “Cell proliferation”.
Figure 4. BP clusters of 133 upregulated DEGs overlapped with IMPAGT. When PPI network was clustered with BP traits, it presented a tendency to cluster into three groups: “Apoptotic process”, “Protein phosphorylation”, and “Cell proliferation”.
Cimb 44 00350 g004
Table
Table 1. Clinicopathologic characteristics of 21 cervical cancer samples. All 21 patients were staged following revised 2009 FIGO (International Federation of Gynecology and Obstetrics) staging system. (n, %) indicates the numbers and percentage of samples in each parameter out of total 21.
Table 1. Clinicopathologic characteristics of 21 cervical cancer samples. All 21 patients were staged following revised 2009 FIGO (International Federation of Gynecology and Obstetrics) staging system. (n, %) indicates the numbers and percentage of samples in each parameter out of total 21.
VariablesN = 21
Age (years)47.4 ± 11.0
FIGO stage (n, %)
 IB111 (52.4)
 IB22 (9.5)
 IIA24 (19.0)
 IIB4 (19.0)
Histology (n, %)
 Squamous cell carcinoma15 (71.4)
 Adenocarcinoma5 (23.8)
 Adenosquamous carcinoma1 (4.8)
HPV type (n, %)
 HPV 1616 (52.4)
 HPV 181 (4.8)
 HPV 311 (4.8)
 HPV 331 (4.8)
 HPV 351 (4.8)
 Negative infection3 (14.3)
 Unknown5 (23.8)
Primary tumor size (cm)3.3 ± 1.6
Lymphovascular invasion (n, %)15 (71.4)
Deep stromal invasion (n, %)18 (85.7)
Parametrial invasion (n, %)8 (38.1)
Positive vaginal margin (n, %)1 (4.8)
Lymph node metastasis (n, %)11 (52.4)
High tumor budding (≥3, n, %)15 (71.4)
Recurrence (n, %)6 (28.6)
Death (n, %)4 (19.0)
FIGO = International Federation of Gynecology and Obstetrics; HPV = human papilloma virus.
Table
Table 2. A table of significantly enriched pathways. Pathways with more than 10 DEGs were considered as significantly enriched pathways. A list of significantly enriched pathways from up- and downregulated DEGs are presented with KEGG pathway ID and numbers of DEGs involved in each pathway.
Table 2. A table of significantly enriched pathways. Pathways with more than 10 DEGs were considered as significantly enriched pathways. A list of significantly enriched pathways from up- and downregulated DEGs are presented with KEGG pathway ID and numbers of DEGs involved in each pathway.
Up/Down DEGsKEGG Pathway IDPathwayNumber of DEGs Involved
UpHsa 05200Pathway in Cancer20
Hsa 04612Antigen processing and presentation13
Hsa 05202Transcriptional misregulation in Cancer13
Hsa 04010MAPK Signaling Pathway12
Hsa 04621NOD-like receptor signaling Pathway12
Hsa 04514Cell adhesion molecules11
Hsa 04060Cytokine-cytokine receptor interaction 10
Hsa 04151PI3K/Akt signaling pathway9
DownHsa 05200Pathway in Cancer10
Table
Table 3. Significant BP terms of 133 upregulated DEGs overlapped with IMPAGT.
Table 3. Significant BP terms of 133 upregulated DEGs overlapped with IMPAGT.
GO IDBP Termsp-Value
GO: 0007165Signal transduction1.3 × 10−4
GO: 0006955Immune response 9.8 × 10−10
GO: 0008284Positive regulation of cell proliferation1.3 × 10−2
GO: 0043065Positive regulation of apoptotic process7.2 × 10−3
GO: 0042981Regulation of apoptotic process6.0 × 10−3
GO: 0042127Regulation of cell proliferation5.2 × 10−3
GO: 0001934Positive regulation of protein phosphorylation 1.4 × 10−2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Choi, Y.; Park, N.J.-Y.; Le, T.M.; Lee, E.; Lee, D.; Nguyen, H.D.T.; Cho, J.; Park, J.-Y.; Han, H.S.; Chong, G.O. Immune Pathway and Gene Database (IMPAGT) Revealed the Immune Dysregulation Dynamics and Overactivation of the PI3K/Akt Pathway in Tumor Buddings of Cervical Cancer. Curr. Issues Mol. Biol. 2022, 44, 5139-5152. https://doi.org/10.3390/cimb44110350

AMA Style

Choi Y, Park NJ-Y, Le TM, Lee E, Lee D, Nguyen HDT, Cho J, Park J-Y, Han HS, Chong GO. Immune Pathway and Gene Database (IMPAGT) Revealed the Immune Dysregulation Dynamics and Overactivation of the PI3K/Akt Pathway in Tumor Buddings of Cervical Cancer. Current Issues in Molecular Biology. 2022; 44(11):5139-5152. https://doi.org/10.3390/cimb44110350

Chicago/Turabian Style

Choi, Yeseul, Nora Jee-Young Park, Tan Minh Le, Eunmi Lee, Donghyeon Lee, Hong Duc Thi Nguyen, Junghwan Cho, Ji-Young Park, Hyung Soo Han, and Gun Oh Chong. 2022. "Immune Pathway and Gene Database (IMPAGT) Revealed the Immune Dysregulation Dynamics and Overactivation of the PI3K/Akt Pathway in Tumor Buddings of Cervical Cancer" Current Issues in Molecular Biology 44, no. 11: 5139-5152. https://doi.org/10.3390/cimb44110350

APA Style

Choi, Y., Park, N. J. -Y., Le, T. M., Lee, E., Lee, D., Nguyen, H. D. T., Cho, J., Park, J. -Y., Han, H. S., & Chong, G. O. (2022). Immune Pathway and Gene Database (IMPAGT) Revealed the Immune Dysregulation Dynamics and Overactivation of the PI3K/Akt Pathway in Tumor Buddings of Cervical Cancer. Current Issues in Molecular Biology, 44(11), 5139-5152. https://doi.org/10.3390/cimb44110350

Article Metrics

Back to TopTop