Kai Hong1,2
Yingjue Zhang3
Lingli Yao1,2Jiabo Zhang1Xianneng Sheng1Lihua Song1,2
Yu Guo1*
Yangyang Guo1*- 1Department of Thyroid and Breast Surgery, Ningbo City First Hospital, Ningbo, Zhejiang, China
- 2Medicine School, Ningbo University, Ningbo, Zhejiang, China
- 3Department of Molecular Pathology, Division of Health Sciences, Graduate School of Medicine, Osaka University, Suita, Osaka, Japan
Background: Understanding interior molecular mechanisms of tumorigenesis and cancer progression contributes to antitumor treatments. The angiotensin II receptor-associated protein (AGTRAP) has been confirmed to be related with metabolic products in metabolic diseases and can drive the progression of hepatocellular carcinoma and colon carcinoma. However, functions of AGTRAP in other kinds of cancers are unclear, and a pan-cancer analysis of AGTRAP has not been carried out.
Methods and materials: We downloaded data from The Cancer Genome Atlas and Genotype-Tissue Expression dataset and The Human Protein Atlas databases and then used R software (version 4.1.1) and several bioinformatic tools to conduct the analysis.
Results: In our study, we evaluated the expression of AGTRAP in cancers, such as high expression in breast cancer, lung adenocarcinoma, and glioma and low expression in kidney chromophobe. Furthermore, our study revealed that high expression of AGTRAP is significantly related with poor prognosis in glioma, liver cancer, kidney chromophobe, and so on. We also explored the putative functional mechanisms of AGTRAP across pan-cancer, such as endoplasmic reticulum pathway, endocytosis pathway, and JAK-STAT signaling pathway. In addition, the connection between AGTRAP and tumor microenvironment, tumor mutation burden, and immune-related genes was proven.
Conclusion: Our study provided comprehensive evidence of the roles of AGTRAP in different kinds of cancers and supported the relationship of AGTRAP and tumorous immunity.
Introduction
Cancer is one of the major causes of death in humans. Treatments for various types of cancers have become one of the biggest parts of healthcare spending. Publications reported that abnormal expression and mutation of certain genes including methylation-related genes, immune checkpoint genes, immunomodulatory genes, and mismatch repair (MMR) genes could promote the progression of cancers (Chowell et al., 2021; An and Duan, 2022; Bennett et al., 2022; Kraehenbuehl et al., 2022). In addition, in the past several decades, the tumor microenvironment (TME) has been demonstrated to be a vital factor for cancerous development, such as immune cell infiltration or Hashimoto’s thyroiditis for cancer progression (Zeng et al., 2018; Gan et al., 2021; Mao et al., 2021). Relaying on abundant online databases, many gene-related analyses of cancers were conducted, and a great number of cancer-related genes have been found (Kim et al., 2015; Song et al., 2017; Barneh et al., 2019). Pan-cancer analysis is a valuable way to explore functions and mechanisms of specific genes in various cancers. According to pan-cancer analyses, more detailed research can be designed and conducted.
The renin–angiotensin system (RAS) has been demonstrated to regulate cardiovascular and body fluid homeostasis which causes diseases such as type 2 diabetes mellitus, hypertension, and dyslipidemia (Arendse et al., 2019). Ang II type 1 receptor (AT1R) is an important receptor, and it interacts with the AT1R-associated protein which is encoded by AGTRAP gene, playing a vital role in activating RAS (Daviet et al., 1999; Lopez-Ilasaca et al., 2003). Furthermore, some studies supported that AGTRAP is also associated with tumor progression. Liu et al. (2021) concluded that AGTRAP expression in hepatocellular carcinoma samples is higher than that in normal tissues and adjacent tissues, which potentially leads to a poor prognosis through regulating nuclear factor-kappa B subunit 1 (NF-KB) and mitogen-activated protein kinase (MAPK) signaling pathways. Sanz-Pamplona et al. (2016) reported that overexpression of AGTRAP was found in colon carcinoma which was related to poor outcome in patients. In addition, AGTRAP-BRAF fusion was validated through qRT-PCR and fluorescence in situ hybridization (FISH) analysis in gastric cancer (Figueiredo et al., 2017). Still, the underlying functions and mechanisms of AGTRAP are unclear. Further explorations are urgently needed to reveal the precise functions and mechanisms of AGTRAP in cancer patients.
TME is made up of tumor, stromal, endothelial, and immune cells (Downs-Canner et al., 2022). According to previous studies, the component immune cells play an important role in controlling cancer progression (Lam et al., 2021). Hence, more and more researchers tend to apply TME to predict prognosis and drug response. DNA MMR is a pathway of correcting DNA strands which are wrongly synthesized (Mas-Ponte and Supek, 2020). Studies presented that a decrease in MMR activities is related to high microsatellite instability (MSI-H) and high tumor burden (TMB) of tumors (Mandal et al., 2019; Pestana et al., 2022). Simultaneously, in MSI-H tumors, increased tumor-infiltrating lymphocytes (TILs) can be detected (Michael-Robinson et al., 2001). Importantly, researchers have demonstrated how cancers with MMR deficiency benefit from the immune checkpoint blockade can be predicted according to the diversity of MSI (Sahin et al., 2019). In addition, metabolic recombination of tumor cells is a vital method that cancers take to meet their need for growth and their homeostasis in TME (Lunt and Vander Heiden, 2011; Scharping et al., 2021). An and Duan, (2022) suggested that RNA methylation modification plays an important role in metabolic reprogramming which can greatly influence the progression of cancers. These novel markers can be used to predict prognoses of cancer and drug responses and to reveal potential mechanisms of other cancer-related molecules.
In order to explore the functions and mechanisms of AGTRAP in cancers which were highly unknown, we extracted data from online databases to conduct a pan-cancer analysis. Our study analyzed the expression level of AGTRAP and whether it influenced survival, potential mechanisms of AGTRAP, and correlation between AGTRAP and TME, immune cell infiltration, RNA methylation, MMR, immune checkpoint, immunoregulator, TMB, and neoantigen.
Materials and methods
Gene and protein expression analysis
First, we analyzed the expression difference of AGTRAP (ENSG00000177674) with the detection data of tumor samples and adjacent normal tissues in The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/) project from Tumor Immune Estimation Resource (TIMER, http://timer.cistrome.org/). TIMER is a resource of immune infiltration, gene expression profile, clinical prognosis, and somatic mutations and copy number alterations for various kinds of cancer (Li et al., 2017). Inputting “AGTRAP” in the “Gene_DE” module, we found that there are some deficiencies in several cancers [e.g., The Cancer Genome Atlas Adrenocortical Carcinoma (TCGA-ACC) and The Cancer Genome Atlas Uterine Serous Carcinoma (TCGA-USC)]. To further demonstrate the abnormal expression of AGTRAP in diverse cancer types, Gene Expression Profiling Interactive Analysis, version 2 tool (GEPIA2, http://gepia2.cancer-pku.cn/#analysis) was used to conduct the analysis for the combination of TCGA and Genotype-Tissue Expression (GTEx) datasets (Tang et al., 2019). TCGA tumor samples were compared to the normal tissues from TCGA and GTEx [p-value cutoff = 0.05, log2 fold change (FC) cutoff = 1]. In order to explore the proteomic expression profile of AGTRAP, we conducted the proteomic expression analysis by UALCAN tool (http://ualcan.path.uab.edu/analysis-prot.html) of the CPTAC (Clinical Proteomic Tumor Analysis Consortium) dataset (Chandrashekar et al., 2017). The keyword “AGTRAP” was input in the “CPTAC analysis” module. The protein expression of AGTRAP in breast cancer, colon cancer, head and neck squamous carcinoma, pancreatic adenocarcinoma, glioblastoma multiforme, and hepatocellular carcinoma was presented, respectively. To discuss AGTRAP expression in different stages of cancer, the “Pathological Stage Plot” module of GEPIA2 was used to visualize the results with violin plots. The log2 [transcripts per million (TPM) +1] was used for log-scale.
Survival analysis
To explore correlations between AGTRAP expression and different kinds of cancer, the log-rank test and a multivariate Cox proportional hazard regression model were performed. The “overall survival (OS),” “disease-specific survival (DSS),” “disease-free survival (DFS),” and “progression-free survival (PFS)” of various cancers were analyzed using the “survival” R package and visualized by forest plots. The hazard ratios (HRs) and 95% confidence intervals (CIs) calculated through Cox analysis were used to identify the survival difference affected by AGTRAP abnormal expression in various types of cancer. The median expression of AGTRAP was used to divide samples into low- and high-expression groups.
Protein–protein interaction network
A protein–protein interaction (PPI) network for AGTRAP was constructed using the STRING website tool (https://string-db.org/) (Szklarczyk et al., 2019). The medium confidence genes were conserved (interaction score ≥0.4). Then, correlation analysis of AGTRAP and its related genes in pan-cancer was visualized using the GEPIA2 tool (Tang et al., 2019). We used Pearson’s test for the correlation analysis, nonlog scale for calculation, and log-scale axis for visualization. Moreover, the expression of AGTRAP-related genes across pan-caner was visualized by the heatmap using the GEPIA2 tool, and log2 (TPM +1) was used for log-scale (Tang et al., 2019). In addition, OS and DFS of AGTRAP-related genes were also analyzed using the GEPIA2 tool (Tang et al., 2019).
Enrichment analysis
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were conducted to analyze the potential functions of AGTRAP (Kanehisa and Goto, 2000; The Gene Ontology Consortium, 2019). First, we selected the most relevant top 200 genes of AGTRAP through the GEPIA2 tool (Supplementary Material S1). Then, functional annotations and enrichment pathways of AGTRAP and its relevant genes were identified. Moreover, gene set enrichment analysis (GSEA) of single AGTRAP was separately performed in different kinds of cancers (Subramanian et al., 2005). Enrichment analysis was conducted using “limma,” “org.Hs.eg.db,” “clusterprofiler,” and “enrichplot” R packages (Yu et al., 2012).
Immune score and cell infiltration analysis
We downloaded the standardized pan-cancer dataset from the University of California, Santa Cruz (UCSC, https://xenabrowser.net/) database including TCGA, Therapeutically Applicable Research to Generate Effective Treatments (TARGET), and GTEx cohorts and then extracted the expression data of AGTRAP gene in each sample (Kent et al., 2002). The results of AGTRAP expression were transformed to log(x+0.001), and the gene expression profiles of each kind of tumors were extracted. Afterward, we evaluated the correlation between AGTRAP and immune scores (ImmuneScore, StromaScore, and MicroenvironmentScore) and immune cell infiltration in different kinds of cancers by xCell algorithms using the “IOBR” R package (Zeng et al., 2021). Finally, the significance was identified using Pearson’s correlation coefficient by the “psych” R package from the 10,180 samples of 44 types of cancers (Revelle and Revelle, 2015).
Immune checkpoint, RNA modification, and immunoregulator correlation analysis
In order to evaluate the correlation between immune checkpoint, RNA modification, and immunoregulation genes and AGTRAP, Spearman correlation analysis was conducted to analyze the data downloaded from the UCSC database (Myers and Sirois, 2004). Regarding the immune checkpoint, inhibitory and stimulatory genes were analyzed, respectively. In addition, the correlation between “N6-methyladenosine (m(6)A),” “N1-methyladenosine (m(1)A),” and “5-methylcytosine (m(5)C)” modifications and AGTRAP was presented, respectively. Moreover, the relationship between immunoregulation genes and AGTRAP was estimated.
Tumor mutational burden, microsatellite instability, neoantigen, and mismatch repair mutation correlation analysis
The mutation data of each kind of cancer were calculated using the “MAF Tools” R package (Mayakonda et al., 2018). MSI scores, neoantigen data, and MMR data were obtained from previous studies. MMR genes included: epithelial cell adhesion molecule (EPCAM), PMS1 homolog 2 (PMS2), MutS homolog 6 (MSH6), MutS protein homolog 2 (MSH2), and MutL protein homolog 1 (MLH1) (Latham et al., 2019). Finally, cancers with a number less than 3 were removed from the results. The correlation was analyzed between TMB, MSI scores, neoantigen, and MMR and AGTRAP expression by Spearman’s method (Revelle and Revelle, 2015).
Immunohistochemistry
The immunohistochemistry results of AGTRAP pathological expression in diverse cancers were obtained from The Human Protein Atlas database (THPA, https://www.proteinatlas.org/) (Thul and Lindskog, 2018).
Cell culture and transfection
The breast, pancreatic, and gastric cancer cells (BT-549, PANC-1, and HGC-27) and normal cells (MCF-10A, HPNE, and GES-1) were obtained from the Cell Bank of the Shanghai Institute of Biochemistry and Cell Biology (Shanghai, China). The cell lines were cultivated in Dulbecco’s modified Eagle’s medium or RPMI-1640, both added with 10% fetal bovine serum (FBS), 100 μg/ml streptomycin, and 100 U/ml penicillin (Gibco). The cells were infected with si-AGTRAP (GenePharma).
qRT-PCR
All cells were washed with PBS and then homogenized in TRIzol reagent (Invitrogen). Total RNA was extracted and reverse transcribed into the cDNA template, and then SYBR Green Real-Time PCR Master Mix Plus (Toyobo) was used for qRT-PCR. β-Actin was used as an endogenous reference gene. The primer sequences for amplification were as follows: forward 5′-TGGGGCTGCATTGTATTCTCA-3′ and reverse 5′-AGCCACCCAGAAACATGCTTA-3′ for AGTRAP and forward 5′-TGACGTGGACATCCGCAAAG-3′ and reverse 5′-CTGGAAGGTGGACAGCGAGG-3′ for β-actin.
Western blot
To measure the protein expression of AGTRAP and its related signaling pathway, Western blot was conducted on breast cancer, pancreatic cancer, and gastric cancer. Cell proteins were extracted in accordance with manufacturer’s instruction. The concentration of proteins was identified using BCA protein assay kit (Beyotime, Shanghai, China). A volume of 35 µL of proteins from each sample was used for electrophoresis. Membranes were incubated overnight with primary antibodies followed by blocking the membranes with 5% skimmed milk. Subsequently, secondary antibodies (1:5,000) were used for further incubation. Finally, the protein bands were visualized through an enhanced chemiluminescent reagent (Thermo Fisher).
The AGTRAP antibody (F-6) (sc-271367) was purchased from Santa Cruz Biotechnology (Heidelberg, Germany). The mTOR antibody (ab2732) and AKT antibody (ab8805) were purchased from Abcam (Beijing, China). The p-mTOR antibody (CST 5536S) and p-AKT antibody (CST 4060S) were purchased from Cell Signaling Technology (CST, United States). The GAPDH antibody (AP0063) was purchased from Bioworl Technology (United States).
Statistical analysis
The Wilcoxon test was used to compare the data of tumor samples and normal tissues. K–M analysis and Cox regression analysis were used to analyze the survivals. In addition, Spearman’s test or Pearson’s test was used to conduct the correlation analyses. p-values less than 0.05 were considered statistically significant. Annotations were used to show statistical significance: *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001. R software (version 4.1.1), GraphPad Prism (version 8.4.3), and several bioinformatical tools such as Sangerbox (http://vip.sangerbox.com/home.html) and TIMER were used to conduct the statistical analyses.
Results
Angiotensin II receptor-associated protein expression in pan-cancer
First, AGTRAP mRNA expression in tumor samples and normal tissues was analyzed across TCGA dataset (Figure 1A). AGTRAP expression in tumor samples of BLCA (bladder urothelial carcinoma), BRCA (breast invasive carcinoma), CHOL (cholangiocarcinoma), COAD (colon adenocarcinoma), ESCA (esophageal carcinoma), HNSC (head and neck squamous cell carcinoma), LIHC (liver hepatocellular carcinoma), LUAD (lung adenocarcinoma), READ (rectum adenocarcinoma), STAD (stomach adenocarcinoma), THCA (thyroid carcinoma), UCEC (uterine corpus endometrial carcinoma) (p < 0.001), KIRC (kidney renal clear cell carcinoma), and KIRP (kidney renal papillary cell carcinoma) (p < 0.05) is significantly higher than that in normal tissues. Moreover, in KICH (kidney chromophobe), AGTRAP expression in tumor samples is significantly lower than that in normal tissues (p < 0.001). Notably, HPV-negative HNSC has a higher AGTRAP expression than HPV-positive HNSC.
FIGURE 1. Expression and protein level of AGTRAP in pan-cancer. (A) mRNA expression levels of AGTRAP in different cancers from TCGA data in TIMER. *p < 0.05, **p < 0.01, ***p < 0.001. (B) Analysis of TCGA and GTEx datasets showed differential mRNA expression of AGTRAP in various cancers. *p < 0.05, **p < 0.01, ***p < 0.001. (C) AGTRAP protein expression level in pan-cancer from CPTAC samples. Z-values represented standard deviations from the median across samples for the given cancer types.
In view of the deficiency of several kinds of normal tissue data, TCGA and GTEx datasets were pooled together for analyzing. As shown in Figure 1B, significantly higher AGTRAP expression in tumor samples than in normal tissues was presented in GBM (glioblastoma multiforme), LGG (brain lower grade glioma), UCEC, BRCA, CESC (cervical squamous cell carcinoma and endo-cervical adenocarcinoma), ESCA, STES (stomach and esophageal carcinoma), KIRP, COAD, PRAD (prostate adenocarcinoma), STAD, HNSC, LIHC, SKCM (skin cutaneous melanoma), BLCA, THCA, READ, OV (ovarian serous cystadenocarcinoma), PAAD (pancreatic adenocarcinoma), TGCT (testicular germ cell tumors), UCS (uterine carcinosarcoma), CHOL (p < 0.0001), LUAD (p < 0.001), and KIPAN (pan-kidney) (p < 0.01). Significantly lower AGTRAP expression in tumor samples than in normal tissues was demonstrated in LUSC (lung squamous cell carcinoma), ALL (acute lymphoblastic leukemia), LAML (acute myeloid leukemia), KICH (p < 0.0001), and PCPG (pheochromocytoma and paraganglioma) (p < 0.05).
Abnormal levels of cancer-related proteins are strongly related to cancer progression. Thus, in addition to the mRNA expression analysis, AGTRAP data from CPTAC was analyzed by comparing the tumor samples to normal tissues. Finally, significant high expression of AGTRAP was found in breast cancer, colon cancer, head and neck squamous carcinoma, pancreatic adenocarcinoma, glioblastoma multiforme, and hepatocellular carcinoma (Figure 1C).
In order to evaluate the expression of AGTRAP in different pathological stages in pan-cancer, we conducted the “Pathological Stage Plot” module in GEPIA2. A significant difference of AGTRAP expression in different pathological stages was demonstrated across various cancer types, such as ESCA, KICH, BLCA, KIRC, LIHC, PAAD, and UCS (p < 0.05) (Figure 2).
FIGURE 2. Stage-dependent expression level of AGTRAP. Stages (stage I, stage II, stage III, and stage IV) of pan-cancer were analyzed using TCGA dataset. The log2 (TPM+1) was used for log-scale.
Survival analysis of abnormal angiotensin II receptor-associated protein expression in pan-cancer
TCGA cancer cases were divided into low- and high-expression subsets according to AGTRAP expression, and then the correlation between AGTRAP expression and prognosis was analyzed. Highly expressed AGTRAP was associated to significantly poor OS in LGG, LIHC, and UVM (uveal melanoma) (p < 0.05) (Figure 3A), significantly poor DFS in LGG and STAD (p < 0.05) (Figure 3B), significantly poor DSS in KIRC, LGG, LIHC, STAD, and UVM (Figure 4A), and significantly poor PFS in LGG, HNSC, LIHC, and STAD (p < 0.05) (Figure 4B). Contrarily, highly expressed AGTRAP showed better OS in MESO (mesothelioma) and THCA (p < 0.05) (Figure 3A), better DFS in COAD (p < 0.05) (Figure 3B), better DSS in BRCA and THCA (p < 0.05) (Figure 4A), and better PFS in BRCA (p < 0.05) (Figure 4B).
FIGURE 3. Correlation between the expression of AGTRAP and prognosis in pan-caner. (A) Relationship of overall survival (OS) and AGTRAP expression in pan-cancer. (B) Relationship of disease-free survival (DFS) and AGTRAP expression in pan-cancer.
FIGURE 4. Correlation between the expression of AGTRAP and prognosis in pan-caner. (A) Relationship of disease-specific survival (DSS) and AGTRAP expression in pan-cancer. (B) Relationship of progression-free survival (PFS) AGTRAP expression in pan-cancer.
Protein–protein interaction network and correlation and survival analyses of angiotensin II receptor-associated protein-related genes
The top 10 related genes of AGTRAP were identified through STRING, namely, AGTR1, JAK2, RAF1, CAMK2A, CAMK2B, CAMK2D, CAMK2G, CLCN6, ARAF, and PITPNC1 (Figure 5A). As shown in Figure 5B, the expression of AGTRAP was positively correlated with AGTR1 and ARAF. A significantly negative correlation was confirmed between AGTRAP and CAMK2A, CAMK2B, CAMK2D, CAMK2G, CLCN6, JAK2, and PITPNC1. In addition, difference in the expression of AGTRAP-related genes was found in various cancers in TCGA pan-cancer datasets (Figure 5C). A significant difference in OS and RFS of the AGTRAP-related genes was presented in the heatmaps. Notably, the same gene shows contrary affects in different cancers. Survival analysis showed that the high expression of AGTRAP-related genes mostly leads to poor prognoses in ACC, BLCA, KICH, LGG, LIHC, LUSC, MESO, SARC (sarcoma), and UVM and leads to better outcomes in CHOL, KIRC, KIRP, and PAAD (Figures 5D,E). In addition, to better represent the expression of AGTRAP-related genes in pan-cancer, we conducted the analysis in the TIMER database. The results showed that the AGTRAP-related genes have significantly abnormal expression in diverse cancer types. (Figure 6).
FIGURE 5. Analysis of top 10 AGTRAP-related genes. (A) AGTRAP interaction network. (B) Expression correlation between AGTRAP and AGTRAP-related genes (AGTR1, JAK2, RAF1, CAMK2A, CAMK2B, CAMK2D, CAMK2G, CLCN6, ARAF, and PITPNC1). (C) Expression of AGTRAP-related genes in different cancer samples and normal tissues. (D) Overall survival (OS) analysis of AGTRAP-related genes in pan-cancer. (E) Disease-free survival (DFS) analysis of AGTRAP-related genes in pan-cancer.
FIGURE 6. Expression of AGTRAP-related genes in pan-cancer. (A) AGTR1. (B) PITPNC1. (C) ARAF. (D) CLCN6. (E) CAMK2B. (F) CAMK2G. (G) CAMK2D. (H) CAMK2A. (I) RAF1. (J) JAK2.
Angiotensin II receptor-associated protein and related gene functional enrichment analysis
To further explore the molecular mechanisms of AGTRAP in pan-cancer, enrichment analysis was conducted using the data of AGTRAP and its top 200 related genes. As shown in Figure 7A, AGTRAP is mostly involved in “protein processing in the endoplasmic reticulum pathway,” “alcoholism pathway,” “endocytosis pathway,” “prostate cancer pathway,” and “amino sugar nucleotide sugar metabolism pathway.” In addition, AGTRAP is mostly enriched in hallmarks such as “androgen response,” “glycolysis,” “protein secretion,” “peroxisome,” and “DNA repair.” Regarding the GO enrichment analysis, “intracellular transport,” “intracellular protein transport,” “cellular component disassembly,” and “golgi vesicle transport” are the most relative functional annotations of biological process; “endoplasmic reticulum,” “organelle subcompartment,” “golgi apparatus,” and “nuclear outer membrane endoplasmic reticulum membrane network” are mostly related to cellular component; “transporter activity,” “cell adhesion molecule binding,” “cadherin binding,” and “cation transmembrane transport activity” are significantly associated to molecular function (Figure 7B).
FIGURE 7. Enrichment analysis of top 200 AGTRAP-related genes. (A) KEGG enrichment analysis of AGTRAP-related genes. (B) GO enrichment analysis of AGTRAP-related genes.
Afterward, enrichment pathways of single AGTRAP in different kinds of cancers were evaluated to further reveal the potential mechanism (Figure 8). From GSEA, enrichment of “olfactory transduction pathway” was found in most types of cancers including DLBC, ESCA, BLCA, BRCA, CESC, COAD, GBM, and KIRC. In addition, the “neuroactive ligand–receptor interaction pathway” was enriched in BLCA, BRCA, CHOL, GBM, HNSC, and KIRC. Notably, the “autoimmune thyroid disease pathway,” “cytosolic DNA sensing pathway,” “olfactory transduction pathway,” “RIG-I-like receptor signaling pathway,” “neuroactive ligand–receptor interaction pathway,” and “taste transduction pathway” were enriched in DLBC, BRCA, and COAD, which means that AGTRAP might have some similar functions in these three kinds of cancers. Importantly, AGTRAP also enriched in some immune-related pathways including “chemokine signaling pathway,” “cytokine–receptor interaction pathway,” “natural killer cell-mediated cytotoxicity pathway,” “autoimmune thyroid disease pathway,” “cytosolic DNA sensing pathway,” “regulation of autophagy pathway,” “RIG-I-like receptor signaling pathway,” “antigen processing and presentation pathway,” “JAK-STAT signaling pathway,” “toll-like receptor signaling pathway,” and “neuroactive ligand–receptor interaction pathway” and some metabolism-related pathways including “taurine and hypotaurine metabolism pathway,” “drug metabolism cytochrome P450 pathway,” “metabolism of xenobiotics by the cytochrome P450 pathway,” “retinol metabolism pathway,” “glycine, serine, and threonine metabolism pathway,” and “arachidonic acid metabolism pathway,” which confirmed that AGTRAP is a gene whose functions are associated with immunity and metabolism.
FIGURE 8. GSEA of AGTRAP in pan-cancer.
Tumor microenvironment correlation analysis
As shown in Figure 9, a significantly positive association was found between AGTRAP and ImmuneScore in LGG, NB (neuroblastoma), BLCA, SARC, UVM, PCPG, KIPAN, LIHC, KICH, LAML, ALL, and GBM; between AGTRAP and StromalScore in LGG, NB, PCPG, COAD, LUSC, and THYM; and between AGTRAP and MicroenvironmentScore in LGG, NB, BLCA, SARC, UVM, PCPG, LAML, ALL, and GBM. Moreover, a significantly negative correlation was confirmed between AGTRAP and ImmuneScore in THCA, BRCA, THYM, OV, STES, and STAD; between AGTRAP and StromalScore in PRAD, KIPAN, LIHC, TGCT, CESC, MESO, OV, STES, and HNSC; and between AGTRAP and MicroenvironmentScore in PRAD, THCA, THYM, OV, STES, HNSC, and STAD.
FIGURE 9. Analysis of correlation between the expression of AGTRAP and immune cell infiltration, ImmuneScore, StromaScore, and MicroenvironmentScore.
Then, immune cell infiltration analysis was conducted to further expose the association between AGTRAP expression and tumor microenvironment in pan-cancer (Figure 9). The heatmap shows that in almost all kinds of cancers, AGTRAP expression is variously related with immune and stromal cells. Notably, LGG, NB, BLCA, SARC, UVM, and PCPG are cancers which have more positive coexpression of AGTRAP and cell infiltration, while OV, THYM, LAML, UCEC, and STES have more negative coexpression. In addition, M2 macrophages, mesangial cells, monocytes, and epithelial cells show positive coexpression in most kinds of cancers and CD4+ T cells including memory and naïve cells, myocytes, naïve B cells, neutrophils, pericytes, plasma cells, and platelets show negative coexpression in most cancers.
Immune checkpoint, RNA methylation, and immunomodulator correlation analysis
The correlation of AGTRAP expression and 60 immune checkpoint genes including inhibitory and stimulatory genes was analyzed (Figure 10). AGTRAP expression in different kinds of cancer is variously associated with immune checkpoint genes. In NB, SARC, UVM, DLBC, GBM, LGG, TGCT, KIRC, COAD, PCPG, SKCM, KIPAN, BLCA, CHOL, THYM, LIHC, PAAD, ACC, and OV, AGTRAP expression was more strongly related with immune checkpoint genes than others. In addition, CD276, TGFB1, C10orf54, IL10, CD274, BTN3A1, BTN3A2, CXCL10, CCL5, GZMA, PRF1, CD40, TNFRSF18, TNFRSF4, TNFRSF14, TNFRSF9, ITGB2, ICAM1, HMGB1, CX3CL1, CD70, ICOS, CD27, CD40LG, ENTPD1, IL1A, CD80, CXCL9, IFNG, TNF, IL1B, and IL2RA were more significantly associated with AGTRAP. From the result, it can obviously be concluded that AGTRAP is more related with the stimulatory immune checkpoint.
FIGURE 10. Analysis of correlation between the expression of AGTRAP and immune checkpoint genes.
Regarding the RNA modification genes, correlation analysis was performed for these genes and AGTRAP across various cancers. As shown in Figure 11, AGTRAP expression in ALL, NB, SARC, HNSC, LUSC, LGG, OV, ACC, KICH, TGCT, BLCA, GBM, CHOL, LIHC, PAAD, ESCA, KIPAN, KIRC, UVM, COAD, WT (Wilms’ tumor), THCA, STAD, and STES relates with more modification genes. Moreover, higher correlation modification genes include: m(1)A writer: TRMT61A; m(1)A readers: TYHDF1 and TYHDF2; m(1)A eraser: ALKBH3; m(5)C writers: NSUN4 and NSUN5; m(5)C reader: ALYREF; m(5)C eraser: TET2; m(6)A writers: KIAA1429; m(6)A readers HNRNPA2B1, HNRNPC, YTHDF1, YTHDF2, and ELAVL1; and m(6)A eraser: ALKBH5.
FIGURE 11. Analysis of correlation between the expression of AGTRAP and RNA modification genes.
Figure 12 shows that the immunoregulatory genes of CESC, LUSC, SKCM, COAD, THYM, STES, NB, ACC, PCPG, KIPAN, KIRC, DLBC, UVM, GBM, LGG, OV, LIHC, PAAD, BLCA, and SARC were highly related with AGTRAP. Importantly, MHC-related genes were more associated with AGTRAP expression than the chemokine-related and receptor-related genes in pan-cancer.
FIGURE 12. Analysis of correlation between the expression of AGTRAP and immunoregulatory genes.
Tumor mutational burden, microsatellite instability, neoantigen, and mismatch repair mutation correlation analysis
TMB, MSI, and neoantigen are the three novel biomarkers related with immunotherapy responses. Correlation analyses were performed between these biomarkers and AGTRAP expression in pan-cancer. The AGTRAP expression level is significantly associated with TMB in CHOL (R < -0.2), ACC, and KICH (R > 0.2) (Figure 13A). Moreover, a significant correlation between AGTRAP and MSI is found in KIPAN and USC (R < -0.2), and GBM (R > 0.2) (Figure 13B). In addition, the relationship between AGTRAP and neoantigen was estimated, and CHOL and TGCT showed a negative correlation (R < -0.2), while DLBC and COAD showed a positive correlation (R > 0.2) (Figure 13C).
FIGURE 13. Analysis of relationship between the expression of AGTRAP and novel biomarkers in order to predict immunotherapy responses. (A) Correlation between the expression of AGTRAP and tumor mutation burden (TMB). (B) Correlation between the expression of AGTRAP and microsatellite instability (MSI). (C) Correlation between the expression of AGTRAP and neoantigen. (D) Correlation between the expression of AGTRAP and five mismatch repair (MMR) genes (EPCAM, PMS2, MSH6, MSH2, and MLH1).
MMR genes play a vital role in maintaining gene stability, such as correcting base mismatches and deletions and insertion errors (Liu et al., 2017). Therefore, coexpression analysis of MMR genes and AGTRAP in pan-cancer was conducted and presented by the heatmap. Figure 13D shows that all five MMR genes we chose are significantly related with various cancers, in which PMS2, MSH6, and MSH2 are more related with AGTRAP than the other MMR genes. In addition, GBM, HNSC, LUSC, PRAD, THCA, and THYM are cancers in which AGTRAP expression is significantly related with all the five MMR genes.
Immunohistochemistry of the angiotensin II receptor-associated protein in pan-cancer
Immunohistochemistry analysis highlighted the expression of AGTRAP in breast cancer, cervix cancer, colorectal cancer, endometrial cancer, glioma, head and neck cancer, liver cancer, lung cancer, lymphoma, melanoma, ovarian cancer, pancreatic cancer, prostate cancer, renal cancer, skin cancer, stomach cancer, testis cancer, thyroid cancer, and urothelial cancer (Figure 14).
FIGURE 14. Immunohistochemistry of AGTRAP expression (normal tissue is on the left and tumor tissue is on the right) (A) Breast cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/breast+cancer). (B) Cervix cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/cervical+cancer). (C) Colorectal cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/colorectal+cancer). (D) Endometrial cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/endometrial+cancer). (E) Glioma (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/glioma). (F) AGTRAP expression in head and neck cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/head+and+neck+cancer). (G) Renal cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/renal+cancer). (H) Skin cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/skin+cancer). (I) Stomach cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/stomach+cancer). (J) Liver cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/liver+cancer). (K) Lung cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/lung+cancer). (L) Lymphoma (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/lymphoma). (M) Ovarian cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/ovarian+cancer). (N) Pancreatic cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/pancreatic+cancer). (O) Prostate cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/prostate+cancer). (P) Testis cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/testis+cancer). (Q) Thyroid cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/thyroid+cancer). (R) Urothelial bladder cancer (https://www.proteinatlas.org/ENSG00000177674-AGTRAP/pathology/urothelial+cancer).
qRT-PCR of the angiotensin II receptor-associated protein in pan-cancer
To further verify the abnormal expression of AGTRAP in pan-cancer, we conducted qRT-PCR to detect the mRNA expression in breast cancer, pancreatic cancer, and gastric cancer. The results showed that AGTRAP is significantly highly expressed in breast cancer, pancreatic cancer, and gastric cancer, which verifies our bioinformatic analysis (Supplementary Figure S1).
Western blot of the angiotensin II receptor-associated protein and the signaling pathway in pan-cancer
In addition, we performed Western blot in pan-cancer. Subsequently, high expression of AGTRAP was identified in breast cancer, pancreatic cancer, and gastric cancer (Supplementary Figure S2A). To explore the mechanism underlying AGTRAP, the protein expression of AKT/mTOR axis was measured by Western blot. The results showed that p-AKT and p-mTOR expression of si-AGTRAP was significantly lower than the expression of control and si-NC (Supplementary Figure S2B).
Discussion
As a protein that regulates visceral obesity-related metabolism, AGTRAP has been widely researched in metabolic disorders (Maeda et al., 2013). Accumulating evidence confirmed that AT1R can regulate the RAS signaling pathway, thereby affecting the handling of renal sodium (Tamura et al., 2022). The overexpression of AT1R would excessively activate the RAS signaling, leading to the retention of sodium which is related to hypertension (Arendse et al., 2019). As a binding protein of the carboxyl-terminal domain, the AT1R-associated protein can suppress the hyperactivation of AT1R (Tamura et al., 2013). Previous studies confirmed that AT1R-associated proteins can inhibit the angiotensin-dependent hypertension by decreasing the sodium reabsorption (Ohsawa et al., 2014). In addition to metabolic functions, several bioinformatic analyses demonstrated that AGTRAP can be used as a component of prognostic signatures in brain lower grade glioma, tongue squamous cell carcinomas, and melanoma, revealing the potential of prognostic ability of AGTRAP (Loftus et al., 2017; Zeng et al., 2019; Xiong et al., 2020). However, details of AGTRAP about expression, survival, mechanism, and treatment response in diverse cancers are highly unclear.
In our study, the high expression of AGTRAP was found in 14 kinds of cancer and low expression in 1 cancer type from TCGA datasets. Moreover, the analysis of combination of TCGA and GTEx datasets suggested that AGTRAP was highly expressed in 26 kinds of cancer and lowly expressed in six cancer types. In addition, the abnormal protein expression of AGTRAP has been confirmed in the CPTAC samples across breast cancer, colon cancer, head and neck cancer, pancreatic cancer, glioblastoma multiforme, and hepatocellular carcinoma. Survival analysis showed that 12 kinds of cancers have poorer OS significantly related with the overexpression of AGTRAP and 2 types of cancers have better OS associated with overexpressed AGTRAP. According to our results, we found that our study is consistent with the past studies which researched AGTRAP in glioma and hepatocellular carcinoma (Xiong et al., 2020; Liu et al., 2021). However, in our study, the high expression of AGTRAP in colon cancer showed better survival which is contrary to a previous study (Sanz-Pamplona et al., 2016). In view of the deficiency of studies of AGTRAP in cancers, relevant research studies are urgently needed.
The PPI network revealed the most related 10 genes of AGTRAP, namely, AGTR1, JAK2, RAF1, CAMK2A, CAMK2B, CAMK2D, CAMK2G, CLCN6, ARAF, and PITPNC1. As shown in Table 1, most of them have been confirmed to influence tumor progression in various cancers, despite some genes showing contrary functions in different cancers. Based on our KEGG analysis, AGTRAP is related with protein processing in ER also called unfolded protein response (UPR) or endoplasmic reticulum (ER) stress (Hetz and Saxena, 2017; Hetz et al., 2020). ER stress can degrade misfolded proteins, decrease protein translation, and induce apoptosis (Wang et al., 2017; You et al., 2021). Researchers have reported that CHOP, Bcl-2 family, caspase-12, and JNK were components of ER stress related to the apoptotic pathway, whose cancer-related mechanisms have been demonstrated in various cancers (Ashkenazi et al., 2017; Kim et al., 2018; Tan et al., 2006). Therefore, there might be some potential mechanisms between these molecules and AGTRAP in cancers that need to be further explored. Furthermore, alcoholism is another related pathway of AGTRAP in cancer progression. From past studies, associations between alcoholism and various kinds of cancer have been widely explored. Lin et al. (2020) found that longer alcoholism history is related with the higher risk of colorectal cancer. In the study of Sun et al. (2020), an anti-alcoholism drug disulfiram showed the effect of causing immunogenic cell death in radiation-resistant breast cancer stem cells. Interestingly, Zhang et al. (2019) supported that disulfiram can induce autophagy-dependent apoptosis of pancreatic and breast cancer cells by evoking ER stress, which is accomplished through activation of IRE1α. Therefore, there might be some synergetic effects between ER stress and alcoholism, through which AGTRAP plays its role in cancerous development. Moreover, our analysis showed the relationship between AGTRAP and endocytosis. As a method of capturing materials from the extracellular environment, there takes place a special endocytic process named macropinocytosis, which can degrade protein to supply amino acids for cancer cell development and remove damaged material from the plasma membrane to restore membrane integrity (Commisso et al., 2013). Notably, it is this function of catching materials that makes it possible to mediate targeting therapy (Xiao et al., 2021). Nanomedicines (NMs) can effectively increase therapeutic efficacy and decrease toxicity for cancer therapy. Through different endocytosis processes, NMs can enter cells (Patel et al., 2019; Wang et al., 2020). If it is possible that drugs can selectively enter into cancer cells instead of delivering into healthy cells, cancer treatment would be highly efficient with low toxicity. In light of the KEGG analysis, AGTRAP may have potential mechanisms related with endocytosis in cancers, perhaps mediating the endocytosis process, which means that it is likely to use AGTRAP to reveal selectivity of drugs in cancer cells, which is worthy to be explored. In addition, AGTRAP is associated with the prostate cancer pathway according to enrichment analysis, which is consistent with our expression analysis that AGTRAP is highly expressed in prostate cancer. We also found the correlation between AGTRAP and androgen response, and since prostate cancer is strongly related with androgen, we believe that our analysis revealed the possible connection between AGTRAP and prostate cancer. According to the GO analysis, the most related biological process of AGTRAP is intracellular transport, which plays a role in direct movement of materials within cells (Clark, 2020). In addition, as the highest related molecular function of AGTRAP, transporter activity including lipid transporter activity, protein transporter activity, and transmembrane transporter activity have the function of enabling the directed movement of substances (Fuji et al., 2009; Nishimura and Stefan, 2020). Furthermore, combining the cellular component result, in which ER is the most enriched annotation, we inferred that AGTRAP is possible to influence cancer progression through ER-related material transportation. In our GSEA, we demonstrated that AGTRAP enriched in many immune-related pathways: “natural killer cell-mediated cytotoxicity pathway,” “chemokine signaling pathway,” “RIG-I-like receptor signaling pathway,” “toll-like receptor signaling pathway,” and “JAK-STAT signaling pathway”; and some metabolism-related pathways: “taurine and hypotaurine metabolism pathway” and “glycine, serine, and threonine metabolism pathway.” Evidence has confirmed that natural killer (NK) cells play a vital role in immunosurveillance and immunoediting, correlating with other immunity and assisting immunotherapy (Malmberg et al., 2017; Pugh-Toole et al., 2022). Veneziani et al. (2022) found that the function of NK cells can be strengthened by toll-like receptor 8 agonists, meaning that the interaction between NK cells and toll-like receptor signaling pathway, both of which are potential pathways that AGTRAP acts on. In addition, the enrichment of the chemokine signaling pathway was consisted in the abundant correlation between immunomodulators and AGTRAP, meaning that AGTRAP has functions in regulating immunomodulators, promoting tumor progression. In addition, RIG-I-like receptor and toll-like receptor are two nucleic acid sensors whose actions for innate immunity and potential for immunotherapy have been identified (Kanzler et al., 2007; Loo and Gale, 2011). Furthermore, previous studies proved that taurine, serine–glycine, and threonine metabolism could promote tumorigenesis and cancer cell proliferation in various kinds of cancers through diverse pathways (Locasale, 2013; Kim et al., 2014; Lin et al., 2018). In conclusion, our study indicated that AGTRAP might regulate these metabolism-related and immune-related pathways, leading to tumor development.
TABLE 1. AGTRAP-related genes.
A total of three novel markers used to predict immunotherapy-responses including TMB, MSI, and neoantigen were evaluated in pan-cancer. Immune checkpoint inhibitors (ICIs) have been demonstrated to be efficient in treating various kinds of cancers (Darvin et al., 2018). Currently, despite having found a lot of biomarkers for ICI and OS such as different immune checkpoints, PD-L1 is the only one which has been widely confirmed, such as the validated effectiveness of pembrolizumab for nonsmall cell lung cancer selected with PD-L1 > 50% and the combination of atezolizumab and nab-paclitaxel for metastatic triple negative breast cancer with PD-L1 > 1% (Schmid et al., 2018; Reck et al., 2019). Nevertheless, the prediction ability of PD-L1 is imperfect because of its heterogeneity and lability (Kim and Chung, 2019). Therefore, TMB, which is defined as the number of nonsynonymous exonic mutations per megabase, is used to predict the response of ICI (Fumet et al., 2020). Another biomarker, MSI-H, is demonstrated to have better prognoses in colon cancer patients (Latham et al., 2019). Moreover, neoantigens as mutant peptides are produced through antitumor immune response, which is elicited by somatic alterations where mutated peptide fragments are generated and then presented on the class I major histocompatibility complex (MHC-1) (Yarchoan et al., 2017). Since neoantigens can impel anticancer immunity and is tumor-specific, it has been regarded as potential therapeutic targets (Schumacher and Schreiber, 2015). Through our correlation analysis, the overexpression of AGTRAP predicts better immunotherapy responses in ACC, KICH, GBM, DLBC, and COAD. In view of the prediction ability of TMB, MSI, and neoantigen in diverse cancers, potential of AGTRAP for predicting drug responses combined with these novel biomarkers in pan-cancer is promising. Immune cells play important roles in tumorigenesis. In cancerous progression, immune cells can effectively eliminate tumor cells in early stages, but cancer cells gradually generate the tumor immune escape and influence the state of immune cells (Bremnes et al., 2011; Gonzalez et al., 2018). Among immune cells, neutrophils, natural killer cells, macrophages, dendritic cells, B cells, and T cells are those which are mostly related with cancer development. From our results, we found that correlations between AGTRAP expression and diverse immune cells can be found in almost all kinds of cancers. In addition, AGTRAP expression is significantly related with ImmuneScore in 19 types of cancers and with StromalScore in 17 types of cancers, which also reveals the relationship between AGTRAP and cancer immunity. Importantly, we found that AGTRAP expression was positively related with M2 macrophage infiltration among 25 cancer types. M2 macrophage is a kind of tumor-associated macrophages (TAMs) that play a vital role in forming TME (Pan et al., 2020). Previous research confirmed that M2 macrophage has functions of promoting angiogenesis, reconstructing tissue, repairing injury, and facilitating tumorigenesis and progression (Zhu et al., 2017; Annamalai et al., 2018). Chen et al. (2017) demonstrated that increased M2 macrophage was associated with poor prognosis and metastasis in breast cancer. Also, Yin et al. (2016) uncovered that M2 macrophage could promote spheroid formation and tumor growth stimulating metastasis in ovarian cancer. By combining previous studies and our analysis, we could speculate that AGTRAP might increase the M2 macrophage infiltration in TME which leads to tumor development. All of the factors we discussed previously have been used to evaluate the TME. As the biomarker for cancers, TME can be used to predict the responses to treatment and prognoses of patients. Based on our correlation analyses, AGTRAP shows the relationship with TME that makes the precise prediction for drug responses and outcomes possible.
In order to further reveal the mechanism of AGTRAP in influencing cancer progression, we analyzed the correlation between AGTRAP and some vital genes including immune checkpoint, RNA methylation, immunoregulator, and MMR. Studies have demonstrated that MMR plays vital roles in maintaining stability and integrity of the genome in cells (Li, 2008). Abnormality of MMR can augment the mutation frequency of related genes which would promote cancer growth and affect therapy responses (Green et al., 2011). The chemokine receptor, MHC, has been widely confirmed to increase tumorous growth and metastasis (Zlotnik et al., 2011; Nagarsheth et al., 2017; Miles et al., 2021). Bruchard et al. (2022) found that CCL20 can boost ILC3-dependent anticancer immunity and increase cancer sensitivity to immunotherapy. Dai et al. (2021) demonstrated the efficacy of CXCL13+CD8+T-cell infiltration for predicting the prognosis and being an immunotherapeutic target in clear cell renal cell carcinoma. RNA methylation, a common chemical modification in eukaryotic RNAs, is related with tumorous development acting through being catalyzed by writers, being removed by erasers, and interacting with readers (Lan et al., 2019). For instance, METTL14, a writer of m(6)A methylation, demonstrated downregulation in gastric cancer tissues and poor prognosis in gastric cancer patients (Fan et al., 2022). In addition, TRMT6/TRMT61A-mediated m(1)A methylation can drive liver tumorigenesis through cholesterol metabolism (Wang et al., 2021) and NSUN2-mediated and YBX1-mediated m(5)C methylation confirmed driving bladder carcinoma tumorigenesis (Chen et al., 2019). PD-l/PD-L1 is the most common immune checkpoint and of great significance for cancer, but others should not be ignored. Analyses of cancer-related genes combining with immune checkpoints have huge potentials for researching tumorigenesis and treatments. According to our study, AGTRAP shows great relationships with immune checkpoint, RNA methylation, immunoregulator, and MMR genes in pan-cancer, which means that there possibly exist some interactions between them influencing cancer progression. However, this field is almost blank and needs further exploration.
To further verify the protein expression of AGTRAP in pan-cancer, Western blot was also conducted on breast cancer, pancreatic cancer, and gastric cancer. Our results demonstrated that AGTRAP was highly expressed in these cancer types, which is consistent with the analysis across CPTAC dataset. Subsequently, the AKT/mTOR axis was detected both in breast cancer and gastric cancer. The results revealed that p-AKT and p-mTOR protein expression was significantly highly expressed, meaning that the AKT/mTOR axis was the possible downstream of AGTRAP. Previous studies reported that AKT and mTOR proteins function by phosphorylation (Wang et al., 2022; Zhu et al., 2022), and AGTRAP can promote phosphorylation of AKT and mTOR, initiating the AKT/mTOR axis. Accumulating evidence has demonstrated the functions of p-AKT and p-mTOR in promoting the progression of breast cancer and gastric cancer (Xu et al., 2021; Khan et al., 2022). Our study reveals that a novel biomarker, AGTRAP, can regulate the AKT/mTOR axis.
Nevertheless, the study still has some limitations. First, although with a global analysis for pan-cancer, details about each cancer type are insufficient. Thus, we will conduct deep analysis for specific cancer in our future studies. In addition, as a pan-cancer bioinformatic study, although with reasonable logic and comprehensive analysis for diverse cancers, combining with the verification of AGTRAP expression in three cancer types, deeper mechanisms of AGTRAP in each cancer type still need to be explored in the future.
In conclusion, our analysis is the first to explore the expression difference of AGTRAP between pan-cancer samples and their normal tissues, analyzing the prognostic impact and potential mechanisms of AGTRAP and its related genes in different kinds of cancers and revealing the correlation between AGTRAP and TME, TMB, MSI, neoantigen, and immune-related genes, which can assist predicting immunotherapy responses. In this study, we showed that AGTRAP and its related genes show abnormal expressions in different types of cancers that might affect tumor progression through regulating immune-related and metabolism-related pathways. Moreover, the abnormal expression of AGTRAP is associated with several immune-related biomarkers and could effectively predict drug responses, meaning that AGTRAP has the potential to become a target for antitumor immunotherapy.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
KH, YZ, and YYG designed the study and drafted the manuscript. KH, YZ, LY, JZ, and YYG wrote the manuscript. JZ, XS, LS, and YG searched the publications and collected the data. KH and YYG analyzed the data. The final manuscript was reviewed by all authors and approved for publication.
Acknowledgments
The authors thank all the members for contributing to this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.913684/full#supplementary-material
References
An, Y., and Duan, H. (2022). The role of m6A RNA methylation in cancer metabolism. Mol. Cancer 21 (1), 14. doi:10.1186/s12943-022-01500-4
Annamalai, R. T., Turner, P. A., Carson, W. F., Levi, B., Kunkel, S., and Stegemann, J. P. (2018). Harnessing macrophage-mediated degradation of gelatin microspheres for spatiotemporal control of BMP2 release. Biomaterials 161, 216–227. doi:10.1016/j.biomaterials.2018.01.040
Arendse, L. B., Danser, A. H. J., Poglitsch, M., Touyz, R. M., Burnett, J. C., Llorens-Cortes, C., et al. (2019). Novel therapeutic approaches targeting the renin-angiotensin system and associated peptides in hypertension and heart failure. Pharmacol. Rev. 71 (4), 539–570. doi:10.1124/pr.118.017129
Ashkenazi, A., Fairbrother, W. J., Leverson, J. D., and Souers, A. J. (2017). From basic apoptosis discoveries to advanced selective BCL-2 family inhibitors. Nat. Rev. Drug Discov. 16 (4), 273–284. doi:10.1038/nrd.2016.253
Barneh, F., Mirzaie, M., Nickchi, P., Tan, T. Z., Thiery, J. P., Piran, M., et al. (2019). Integrated use of bioinformatic resources reveals that co-targeting of histone deacetylases, IKBK and SRC inhibits epithelial-mesenchymal transition in cancer. Brief. Bioinform 20 (2), 717–731. doi:10.1093/bib/bby030
Bennett, J. A., Fleming, G. F., Kurnit, K. C., Mills, K. A., and vanWeelden, W. J. (2022). Tumor board presentation of a woman with metastatic, hormone receptor-positive, mismatch repair-deficient endometrial cancer. CA Cancer J. Clin. 72, 102. doi:10.3322/caac.21715
Bremnes, R. M., Al-Shibli, K., Donnem, T., Sirera, R., Al-Saad, S., Andersen, S., et al. (2011). The role of tumor-infiltrating immune cells and chronic inflammation at the tumor site on cancer development, progression, and prognosis: emphasis on non-small cell lung cancer. J. Thorac. Oncol. 6 (4), 824–833. doi:10.1097/jto.0b013e3182037b76
Bruchard, M., Geindreau, M., Perrichet, A., Truntzer, C., Ballot, E., Boidot, R., et al. (2022). Recruitment and activation of type 3 innate lymphoid cells promote antitumor immune responses. Nat. Immunol. 23 (2), 262–274. doi:10.1038/s41590-021-01120-y
Chandrashekar, D. S., Bashel, B., Balasubramanya, S. A. H., Creighton, C. J., Ponce-Rodriguez, I., Chakravarthi, B. V. S. K., et al. (2017). UALCAN: A portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia 19 (8), 649–658. doi:10.1016/j.neo.2017.05.002
Chen, Y., Zhang, S., Wang, Q., and Zhang, X. (2017). Tumor-recruited M2 macrophages promote gastric and breast cancer metastasis via M2 macrophage-secreted CHI3L1 protein. J. Hematol. Oncol. 10 (1), 36. doi:10.1186/s13045-017-0408-0
Chen, X., Li, A., Sun, B. F., Yang, Y., Han, Y. N., Yuan, X., et al. (2019). 5-methylcytosine promotes pathogenesis of bladder cancer through stabilizing mRNAs. Nat. Cell Biol. 21 (8), 978–990. doi:10.1038/s41556-019-0361-y
Chowell, D., Yoo, S. K., Valero, C., Pastore, A., Krishna, C., Lee, M., et al. (2021). Improved prediction of immune checkpoint blockade efficacy across multiple cancer types. Nat. Biotechnol. 40, 499. doi:10.1038/s41587-021-01070-8
Clark, B. J. (2020). The START-domain proteins in intracellular lipid transport and beyond. Mol. Cell. Endocrinol. 504, 110704. doi:10.1016/j.mce.2020.110704
Commisso, C., Davidson, S. M., Soydaner-Azeloglu, R. G., Parker, S. J., Kamphorst, J. J., Hackett, S., et al. (2013). Macropinocytosis of protein is an amino acid supply route in Ras-transformed cells. Nature 497 (7451), 633–637. doi:10.1038/nature12138
Dai, S., Zeng, H., Liu, Z., Jin, K., Jiang, W., Wang, Z., et al. (2021). Intratumoral CXCL13+CD8+T cell infiltration determines poor clinical outcomes and immunoevasive contexture in patients with clear cell renal cell carcinoma. J. Immunother. Cancer 9 (2), e001823. doi:10.1136/jitc-2020-001823
Darvin, P., Toor, S. M., Sasidharan Nair, V., and Elkord, E. (2018). Immune checkpoint inhibitors: recent progress and potential biomarkers. Exp. Mol. Med. 50 (12), 1–11. doi:10.1038/s12276-018-0191-1
Daviet, L., Lehtonen, J. Y., Tamura, K., Griese, D. P., Horiuchi, M., and Dzau, V. J. (1999). Cloning and characterization of ATRAP, a novel protein that interacts with the angiotensin II type 1 receptor. J. Biol. Chem. 274 (24), 17058–17062. doi:10.1074/jbc.274.24.17058
Downs-Canner, S. M., Meier, J., Vincent, B. G., and Serody, J. S. (2022). B cell function in the tumor microenvironment. Annu. Rev. Immunol. 40, 169. doi:10.1146/annurev-immunol-101220-015603
Fan, H. N., Chen, Z. Y., Chen, X. Y., Chen, M., Yi, Y. C., Zhu, J. S., et al. (2022). METTL14-mediated m6A modification of circORC5 suppresses gastric cancer progression by regulating miR-30c-2-3p/AKT1S1 axis. Mol. Cancer 21 (1), 51. doi:10.1186/s12943-022-01521-z
Figueiredo, C., Camargo, M. C., Leite, M., Fuentes-Pananá, E. M., Rabkin, C. S., and Machado, J. C. (2017). Pathogenesis of gastric cancer: Genetics and molecular classification. Curr. Top. Microbiol. Immunol. 400, 277–304. doi:10.1007/978-3-319-50520-6_12
Fuji, K., Miwa, K., and Fujiwara, T. (2009). The intracellular transport of transporters: membrane trafficking of mineral transporters. Curr. Opin. Plant Biol. 12 (6), 699–704. doi:10.1016/j.pbi.2009.09.006
Fumet, J. D., Truntzer, C., Yarchoan, M., and Ghiringhelli, F. (2020). Tumour mutational burden as a biomarker for immunotherapy: Current data and emerging concepts. Eur. J. Cancer 131, 40–50. doi:10.1016/j.ejca.2020.02.038
Gan, X. X., Li, Y. Y., Li, S. J., Mo, S. S., Feng, J. H., Shen, F., et al. (2021). Significance of DMBT1 in papillary thyroid carcinoma concurrent with Hashimoto's thyroiditis. Front. Oncol. 11, 680873. doi:10.3389/fonc.2021.680873
Gonzalez, H., Hagerling, C., and Werb, Z. (2018). Roles of the immune system in cancer: from tumor initiation to metastatic progression. Genes Dev. 32 (19-20), 1267–1284. doi:10.1101/gad.314617.118
Green, B., Belcheva, A., Nepal, R. M., Boulianne, B., and Martin, A. (2011). The mismatch repair pathway functions normally at a non-AID target in germinal center B cells. Blood 118 (11), 3013–3018. doi:10.1182/blood-2011-03-345991
Hetz, C., and Saxena, S. (2017). ER stress and the unfolded protein response in neurodegeneration. Nat. Rev. Neurol. 13 (8), 477–491. doi:10.1038/nrneurol.2017.99
Hetz, C., Zhang, K., and Kaufman, R. J. (2020). Mechanisms, regulation and functions of the unfolded protein response. Nat. Rev. Mol. Cell Biol. 21 (8), 421–438. doi:10.1038/s41580-020-0250-z
Kanehisa, M., and Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28 (1), 27–30. doi:10.1093/nar/28.1.27
Kanzler, H., Barrat, F. J., Hessel, E. M., and Coffman, R. L. (2007). Therapeutic targeting of innate immunity with Toll-like receptor agonists and antagonists. Nat. Med. 13 (5), 552–559. doi:10.1038/nm1589
Kent, W. J., Sugnet, C. W., Furey, T. S., Roskin, K. M., Pringle, T. H., Zahler, A. M., et al. (2002). The human genome browser at UCSC. Genome Res. 12 (6), 996–1006. doi:10.1101/gr.229102
Khan, N. H., Wang, D., Wang, W., Shahid, M., Khattak, S., Ngowi, E. E., et al. (2022). Pharmacological inhibition of endogenous hydrogen sulfide attenuates breast cancer progression. Molecules 27 (13), 4049. doi:10.3390/molecules27134049
Kim, H., and Chung, J. H. (2019). PD-L1 testing in non-small cell lung cancer: Past, present, and Future. J. Pathol. Transl. Med. 5354 (42), 199196–199206. doi:10.4132/jptm.2019.04.24
Kim, J. Y., Banerjee, T., Vinckevicius, A., Luo, Q., Parker, J. B., Baker, M. R., et al. (2014). A role for WDR5 in integrating threonine 11 phosphorylation to lysine 4 methylation on histone H3 during androgen signaling and in prostate cancer. Mol. Cell 54 (4), 613–625. doi:10.1016/j.molcel.2014.03.043
Kim, S., Rhee, J. K., Yoo, H. J., Lee, H. J., Lee, E. J., Lee, J. W., et al. (2015). Bioinformatic and metabolomic analysis reveals miR-155 regulates thiamine level in breast cancer. Cancer Lett. 357 (2), 488–497. doi:10.1016/j.canlet.2014.11.058
Kim, T. W., Lee, S. Y., Kim, M., Cheon, C., and Ko, S. G. (2018). Kaempferol induces autophagic cell death via IRE1-JNK-CHOP pathway and inhibition of G9a in gastric cancer cells. Cell Death Dis. 9 (9), 875. doi:10.1038/s41419-018-0930-1
Kraehenbuehl, L., Weng, C. H., Eghbali, S., Wolchok, J. D., and Merghoub, T. (2022). Enhancing immunotherapy in cancer by targeting emerging immunomodulatory pathways. Nat. Rev. Clin. Oncol. 19 (1), 37–50. doi:10.1038/s41571-021-00552-7
Lam, H., McNeil, L. K., Starobinets, H., DeVault, V. L., Cohen, R. B., Twardowski, P., et al. (2021). An empirical antigen selection method identifies neoantigens that either elicit broad antitumor T-cell responses or drive tumor growth. Cancer Discov. 11 (3), 696–713. doi:10.1158/2159-8290.cd-20-0377
Lan, Q., Liu, P. Y., Haase, J., Bell, J. L., Hüttelmaier, S., and Liu, T. (2019). The critical role of RNA m6A methylation in cancer. Cancer Res. 79 (7), 1285–1292. doi:10.1158/0008-5472.can-18-2965
Latham, A., Srinivasan, P., Kemel, Y., Shia, J., Bandlamudi, C., Mandelker, D., et al. (2019). Microsatellite instability is associated with the presence of lynch syndrome pan-cancer. J. Clin. Oncol. 37 (4), 286–295. doi:10.1200/jco.18.00283
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307
Li, G. M. (2008). Mechanisms and functions of DNA mismatch repair. Cell Res. 18 (1), 85–98. doi:10.1038/cr.2007.115
Lin, Y. H., Wu, M. H., Huang, Y. H., Yeh, C. T., Cheng, M. L., Chi, H. C., et al. (2018). Taurine up-regulated gene 1 functions as a master regulator to coordinate glycolysis and metastasis in hepatocellular carcinoma. Hepatology 67 (1), 188–203. doi:10.1002/hep.29462
Lin, T. C., Chien, W. C., Hu, J. M., Tzeng, N. S., Chung, C. H., Pu, T. W., et al. (2020). Risk of colorectal cancer in patients with alcoholism: a nationwide, population-based nested case-control study. PLoS One 15 (5), e0232740. doi:10.1371/journal.pone.0232740
Liu, D., Keijzers, G., and Rasmussen, L. J. (2017). DNA mismatch repair and its many roles in eukaryotic cells. Mutat. Research/Reviews Mutat. Res. 773, 174–187. doi:10.1016/j.mrrev.2017.07.001
Liu, S., Zhao, W., Li, X., Zhang, L., Gao, Y., Peng, Q., et al. (2021). AGTRAP is a prognostic biomarker correlated with immune infiltration in hepatocellular carcinoma. Front. Oncol. 11, 713017. doi:10.3389/fonc.2021.713017
Locasale, J. W. (2013). Serine, glycine and one-carbon units: cancer metabolism in full circle. Nat. Rev. Cancer 13 (8), 572–583. doi:10.1038/nrc3557
Loftus, S. K., Baxter, L. L., Cronin, J. C., Fufa, T. D., and Pavan, W. J. (2017). Hypoxia-induced HIF1αtargets in melanocytes reveal a molecular profile associated with poor melanoma prognosis. Pigment. Cell Melanoma Res. 30 (3), 339–352. doi:10.1111/pcmr.12579
Loo, Y. M., and Gale, M. (2011). Immune signaling by RIG-I-like receptors. Immunity 34 (5), 680–692. doi:10.1016/j.immuni.2011.05.003
Lopez-Ilasaca, M., Liu, X., Tamura, K., and Dzau, V. J. (2003). The angiotensin II type I receptor-associated protein, ATRAP, is a transmembrane protein and a modulator of angiotensin II signaling. MBoC 14 (12), 5038–5050. doi:10.1091/mbc.e03-06-0383
Lunt, S. Y., and Vander Heiden, M. G. (2011). Aerobic glycolysis: meeting the metabolic requirements of cell proliferation. Annu. Rev. Cell Dev. Biol. 27, 441–464. doi:10.1146/annurev-cellbio-092910-154237
Maeda, A., Tamura, K., Wakui, H., Dejima, T., Ohsawa, M., Azushima, K., et al. (2013). Angiotensin receptor-binding protein ATRAP/Agtrap inhibits metabolic dysfunction with visceral obesity. J. Am. Heart Assoc. 2 (4), e000312. doi:10.1161/JAHA.113.000312
Malmberg, K. J., Carlsten, M., Björklund, A., Sohlberg, E., Bryceson, Y. T., and Ljunggren, H. G. (2017). Natural killer cell-mediated immunosurveillance of human cancer. Seminars Immunol. 31, 20–29. doi:10.1016/j.smim.2017.08.002
Mandal, R., Samstein, R. M., Lee, K. W., Havel, J. J., Wang, H., Krishna, C., et al. (2019). Genetic diversity of tumors with mismatch repair deficiency influences anti-PD-1 immunotherapy response. Science 364 (6439), 485–491. doi:10.1126/science.aau0447
Mao, X., Xu, J., Wang, W., Liang, C., Hua, J., Liu, J., et al. (2021). Crosstalk between cancer-associated fibroblasts and immune cells in the tumor microenvironment: new findings and future perspectives. Mol. Cancer 20 (1), 131. doi:10.1186/s12943-021-01428-1
Mas-Ponte, D., and Supek, F. (2020). DNA mismatch repair promotes APOBEC3-mediated diffuse hypermutation in human cancers. Nat. Genet. 52 (9), 958–968. doi:10.1038/s41588-020-0674-6
Mayakonda, A., Lin, D. C., Assenov, Y., Plass, C., and Koeffler, H. P. (2018). Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 28 (11), 1747–1756. doi:10.1101/gr.239244.118
Michael-Robinson, J. M., Biemer-Hüttmann, A., Purdie, D. M., Walsh, M. D., Simms, L. A., Biden, K. G., et al. (2001). Tumour infiltrating lymphocytes and apoptosis are independent features in colorectal cancer stratified according to microsatellite instability status. Gut 48 (3), 360–366. doi:10.1136/gut.48.3.360
Miles, L. A., Bowman, R. L., Merlinsky, T. R., Csete, I. S., Ooi, A. T., and Durruthy-Durruthy, R. (2021). PCSK9 causes MHC class I lysosomal degradation to promote tumor growth. Cancer Discov. 11 (1), 13. doi:10.1158/2159-8290.CD-RW2020-167
Myers, L., and Sirois, M. J. (2004). “Spearman correlation coefficients, differences between,” in Encyclopedia of statistical sciences. doi:10.1002/0471667196.ess5050
Nagarsheth, N., Wicha, M. S., and Zou, W. (2017). Chemokines in the cancer microenvironment and their relevance in cancer immunotherapy. Nat. Rev. Immunol. 17 (9), 559–572. doi:10.1038/nri.2017.49
Nishimura, T., and Stefan, C. J. (2020). Specialized ER membrane domains for lipid metabolism and transport. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1865 (1), 158492. doi:10.1016/j.bbalip.2019.07.001
Ohsawa, M., Tamura, K., Wakui, H., Maeda, A., Dejima, T., Kanaoka, T., et al. (2014). Deletion of the angiotensin II type 1 receptor-associated protein enhances renal sodium reabsorption and exacerbates angiotensin II-mediated hypertension. Kidney Int. 86 (3), 570–581. doi:10.1038/ki.2014.95
Pan, Y., Yu, Y., Wang, X., and Zhang, T. (2020). Tumor-associated macrophages in tumor immunity. Front. Immunol. 11, 583084. doi:10.3389/fimmu.2020.583084
Patel, S., Kim, J., Herrera, M., Mukherjee, A., Kabanov, A. V., and Sahay, G. (2019). Brief update on endocytosis of nanomedicines. Adv. Drug Deliv. Rev. 144, 90–111. doi:10.1016/j.addr.2019.08.004
Pestana, R. C., Beal, J. R., Parkes, A., Hamerschlak, N., and Subbiah, V. (2022). Impact of tissue-agnostic approvals for patients with sarcoma. Trends Cancer 8 (2), 135–144. doi:10.1016/j.trecan.2021.11.007
Pugh-Toole, M., Nicolela, A. P., Nersesian, S., Leung, B. M., and Boudreau, J. E. (2022). Natural killer cells: the missing link in effective treatment for high-grade serous ovarian carcinoma. Curr. Treat. Options Oncol. 23, 210. doi:10.1007/s11864-021-00929-x
Reck, M., Rodríguez-Abreu, D., Robinson, A. G., Hui, R., Csőszi, T., Fülöp, A., et al. (2019). Updated analysis of KEYNOTE-024: Pembrolizumab versus platinum-based chemotherapy for advanced non-small-cell lung cancer with PD-L1 tumor proportion score of 50% or greater. J. Clin. Oncol. 37 (7), 537–546. doi:10.1200/JCO.18.00149
Revelle, W., and Revelle, M. W. (2015). Package ‘psych’. The comprehensive R archive network, 337–338.
Sahin, I. H., Akce, M., Alese, O., Shaib, W., Lesinski, G. B., El-Rayes, B., et al. (2019). Immune checkpoint inhibitors for the treatment of MSI-H/MMR-D colorectal cancer and a perspective on resistance mechanisms. Br. J. Cancer 121 (10), 809–818. doi:10.1038/s41416-019-0599-y
Sanz-Pamplona, R., Gil-Hoyos, R., López-Doriga, A., Alonso, M. H., Aussó, S., Molleví, D. G., et al. (2016). Mutanome and expression of immune response genes in microsatellite stable colon cancer. Oncotarget 7 (14), 17711–17725. doi:10.18632/oncotarget.7293
Scharping, N. E., Rivadeneira, D. B., Menk, A. V., Vignali, P. D. A., Ford, B. R., Rittenhouse, N. L., et al. (2021). Mitochondrial stress induced by continuous stimulation under hypoxia rapidly drives T cell exhaustion. Nat. Immunol. 22 (2), 205–215. doi:10.1038/s41590-020-00834-9
Schmid, P., Adams, S., Rugo, H. S., Schneeweiss, A., Barrios, C. H., Iwata, H., et al. (2018). Atezolizumab and nab-paclitaxel in advanced triple-negative breast cancer. N. Engl. J. Med. 379 (22), 2108–2121. doi:10.1056/nejmoa1809615
Schumacher, T. N., and Schreiber, R. D. (2015). Neoantigens in cancer immunotherapy. Science 348 (6230), 69–74. doi:10.1126/science.aaa4971
Song, B. N., Kim, S. K., and Chu, I. S. (2017). Bioinformatic identification of prognostic signature defined by copy number alteration and expression of CCNE1 in non-muscle invasive bladder cancer. Exp. Mol. Med. 49 (1), e282. doi:10.1038/emm.2016.120
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Sun, T., Yang, W., Toprani, S. M., Guo, W., He, L., DeLeo, A. B., et al. (2020). Induction of immunogenic cell death in radiation-resistant breast cancer stem cells by repurposing anti-alcoholism drug disulfiram. Cell Commun. Signal 18 (1), 36. doi:10.1186/s12964-019-0507-3
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Tamura, K., Wakui, H., Maeda, A., Dejima, T., Ohsawa, M., Azushima, K., et al. (2013). The physiology and pathophysiology of a novel angiotensin receptor-binding protein ATRAP/Agtrap. Cpd 19 (17), 3043–3048. doi:10.2174/1381612811319170010
Tamura, K., Azushima, K., Kinguchi, S., Wakui, H., and Yamaji, T. (2022). ATRAP, a receptor-interacting modulator of kidney physiology, as a novel player in blood pressure and beyond. Hypertens. Res. 45 (1), 32–39. doi:10.1038/s41440-021-00776-1
Tan, Y., Dourdin, N., Wu, C., De Veyra, T., Elce, J. S., and Greer, P. A. (2006). Ubiquitous calpains promote caspase-12 and JNK activation during endoplasmic reticulum stress-induced apoptosis. J. Biol. Chem. 281 (23), 16016–16024. doi:10.1074/jbc.m601299200
Tang, Z., Kang, B., Li, C., Chen, T., and Zhang, Z. (2019). GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 47 (W1), W556–W560. doi:10.1093/nar/gkz430
The Gene Ontology Consortium (2019). The gene Ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 47 (D1), D330–D338. doi:10.1093/nar/gky1055
Thul, P. J., and Lindskog, C. (2018). The human protein atlas: a spatial map of the human proteome. Protein Sci. 27 (1), 233–244. doi:10.1002/pro.3307
Veneziani, I., Alicata, C., Pelosi, A., Landolina, N., Ricci, B., D'Oria, V., et al. (2022). Toll-like receptor 8 agonists improve NK-cell function primarily targeting CD56brightCD16− subset. J. Immunother. Cancer 10 (1), e003385. doi:10.1136/jitc-2021-003385
Wang, Q., Wang, H., Jia, Y., Pan, H., and Ding, H. (2017). Luteolin induces apoptosis by ROS/ER stress and mitochondrial dysfunction in gliomablastoma. Cancer Chemother. Pharmacol. 79 (5), 1031–1041. doi:10.1007/s00280-017-3299-4
Wang, X., Qiu, Y., Wang, M., Zhang, C., Zhang, T., Zhou, H., et al. (2020). Endocytosis and organelle targeting of nanomedicines in cancer therapy. Int. J. Nanomedicine 15, 9447–9467. doi:10.2147/ijn.s274289
Wang, Y., Wang, J., Li, X., Xiong, X., Wang, J., Zhou, Z., et al. (2021). N1-methyladenosine methylation in tRNA drives liver tumourigenesis by regulating cholesterol metabolism. Nat. Commun. 12 (1), 6314. doi:10.1038/s41467-021-26718-6
Wang, Z., Li, D., Mo, L., Liang, S., Liao, X., Guo, S., et al. (2022). Low-dose cadmium exposure promotes osteoclastogenesis by enhancing autophagy via inhibiting the mTOR/p70S6K1 signaling pathway. Toxicol. Lett. 367 (22), 900150–900183. doi:10.1016/j.toxlet.2022.07.005
Xiao, F., Li, J., Huang, K., Li, X., Xiong, Y., Wu, M., et al. (2021). Macropinocytosis: mechanism and targeted therapy in cancers. Am. J. Cancer Res. 11 (1), 14–30.
Xiong, Z., Xiong, Y., Liu, H., Li, C., and Li, X. (2020). Identification of purity and prognosis related gene signature by network analysis and survival analysis in brain lower grade glioma. J. Cell. Mol. Med. 24 (19), 11607–11612. doi:10.1111/jcmm.15805
Xu, L., Chen, X., Jiang, H., Xu, J., Wang, L., and Sun, Y. (2021). NDUFC1 is upregulated in gastric cancer and regulates cell proliferation, apoptosis, cycle and migration. Front. Oncol. 11, 709044. doi:10.3389/fonc.2021.709044
Yarchoan, M., Johnson, B. A., Lutz, E. R., Laheru, D. A., and Jaffee, E. M. (2017). Targeting neoantigens to augment antitumour immunity. Nat. Rev. Cancer 17 (4), 209–222. doi:10.1038/nrc.2016.154
Yin, M., Li, X., Tan, S., Zhou, H. J., Ji, W., Bellone, S., et al. (2016). Tumor-associated macrophages drive spheroid formation during early transcoelomic metastasis of ovarian cancer. J. Clin. Invest. 126 (11), 4157–4173. doi:10.1172/jci87252
You, K., Wang, L., Chou, C. H., Liu, K., Nakata, T., Jaiswal, A., et al. (2021). QRICH1 dictates the outcome of ER stress through transcriptional control of proteostasis. Science 371 (6524), eabb6896. doi:10.1126/science.abb6896
Yu, G., Wang, L. G., Han, Y., and He, Q. Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Zeng, R., Lyu, Y., Zhang, G., Shou, T., Wang, K., Niu, H., et al. (2018). Positive effect of RORγt on the prognosis of thyroid papillary carcinoma patients combined with Hashimoto's thyroiditis. Am. J. Transl. Res. 10 (10), 3011–3024.
Zeng, H., Li, H., Zhao, Y., Chen, L., and Ma, X. (2019). Transcripto‐based network analysis reveals a model of gene activation in tongue squamous cell carcinomas. Head Neck 41 (12), 4098–4110. doi:10.1002/hed.25952
Zeng, D., Ye, Z., Shen, R., Yu, G., Wu, J., Xiong, Y., et al. (2021). IOBR: Multi-Omics immuno-oncology biological research to decode tumor microenvironment and signatures. Front. Immunol. 12, 687975. doi:10.3389/fimmu.2021.687975
Zhang, X., Hu, P., Ding, S. Y., Sun, T., Liu, L., Han, S., et al. (2019). Induction of autophagy-dependent apoptosis in cancer cells through activation of ER stress: an uncovered anti-cancer mechanism by anti-alcoholism drug disulfiram. Am. J. Cancer Res. 9 (6), 1266–1281.
Zhu, C., Mustafa, D., Zheng, P. P., van der Weiden, M., Sacchetti, A., Brandt, M., et al. (2017). Activation of CECR1 in M2-like TAMs promotes paracrine stimulation-mediated glial tumor progression. Neuro Oncol. 19 (5), 648–659. doi:10.1093/neuonc/now251
Zhu, L., Liu, S., Liao, Y. F., Sheng, Y. M., He, J. C., Cai, Z. X., et al. (2022). Calycosin suppresses colorectal cancer progression by targeting ERβ, upregulating PTEN, and inhibiting PI3K/Akt signal pathway. Cell Biol. Int. 46 (9), 1367–1377. doi:10.1002/cbin.11840
Zlotnik, A., Burkhardt, A. M., and Homey, B. (2011). Homeostatic chemokine receptors and organ-specific metastasis. Nat. Rev. Immunol. 11 (9), 597–606. doi:10.1038/nri3049
Glossary
AGTRAP angiotensin II receptor-associated protein
ALL acute lymphoblastic leukemia
AT1R Ang II type 1 receptor
BLCA bladder urothelial carcinoma
BRCA breast invasive carcinoma
CESC cervical squamous cell carcinoma and endo-cervical adenocarcinoma
CHOL cholangiocarcinoma
CIs confidence intervals
COAD colon adenocarcinoma
CPTAC Clinical Proteomic Tumor Analysis Consortium
DFS disease-free survival
DSS disease-specific survival
EPCAM epithelial cell adhesion molecule
ER endoplasmic reticulum
ESCA esophageal carcinoma
FISH fluorescence in situ hybridization
GBM glioblastoma multiforme,
GEPIA2 Gene Expression Profiling Interactive Analysis, version 2
GO Gene Ontology
GSEA Gene set enrichment analysis
GTEx Genotype-Tissue Expression
HNSC head and neck squamous cell carcinoma
HRs hazard ratios
ICIs immune checkpoint inhibitors
KEGG Kyoto Encyclopedia of Genes and Genomes
KICH kidney chromophobe
KIPAN Pan-kidney
KIRC kidney renal clear cell carcinoma
KIRP kidney renal papillary cell carcinoma
LAML acute myeloid leukemia
LGG brain lower grade glioma
LIHC liver hepatocellular carcinoma
LUAD lung adenocarcinoma
LUSC lung squamous cell carcinoma
m(1)A N1-methyladenosine
m(5)C 5-methylcytosine
m(6)A N6-methyladenosine
MAPK mitogen-activated protein kinase
MESO mesothelioma
MHC-1 Class I major histocompatibility complex
MLH1 MutL protein homolog 1
MMR mismatch repair
MSH2 MutS protein homolog 2
MSH6 MutS homolog 6
MSI-H high microsatellite instability
NB neuroblastoma
NF-KB factor-kappa B subunit 1
NK natural killer
NMs nanomedicines
OS overall survival
OV ovarian serous cystadenocarcinoma
PAAD pancreatic adenocarcinoma
PCPG pheochromocytoma and paraganglioma
PFS progression-free survival
PMS2 PMS1 homolog 2
PPI protein–protein interaction
PRAD prostate adenocarcinoma
RAS renin–angiotensin system
READ rectum adenocarcinoma
SARC sarcoma
SKCM skin cutaneous melanoma
STAD stomach adenocarcinoma
STES stomach and esophageal carcinoma
TAMs tumor-associated macrophages
TARGET Therapeutically Applicable Research to Generate Effective Treatments
TCGA The Cancer Genome Atlas
TCGA-ACC The Cancer Genome Atlas Adrenocortical Carcinoma
TCGA-USC The Cancer Genome Atlas Uterine Serous Carcinoma
TGCT testicular germ cell tumors
THCA thyroid carcinoma
THPA The Human Protein Atlas
TILs tumor-infiltrating lymphocytes
TIMER Tumor Immune Estimation Resource
TMB tumor burden
TME tumor microenvironment
TPM transcripts per million
UCEC uterine corpus endometrial carcinoma
UCS uterine carcinosarcoma
UCSC University of California, Santa Cruz
UPR unfolded protein response
UVM uveal melanoma
WT Wilms’ tumor
Keywords: AGTRAP, pan-cancer, prognosis, tumor microenvironment, immunity
Citation: Hong K, Zhang Y, Yao L, Zhang J, Sheng X, Song L, Guo Y and Guo Y (2022) Pan-cancer analysis of the angiotensin II receptor-associated protein as a prognostic and immunological gene predicting immunotherapy responses in pan-cancer. Front. Cell Dev. Biol. 10:913684. doi: 10.3389/fcell.2022.913684
Received: 06 April 2022; Accepted: 21 July 2022;
Published: 19 August 2022.
Edited by:
Jinhui Liu, Nanjing Medical University, ChinaReviewed by:
Sihan Chen, First Affiliated Hospital of Anhui Medical University, ChinaHongjun Yu, Harbin Medical University, China
Manish Muhuri, Biogen Idec, United States
Copyright © 2022 Hong, Zhang, Yao, Zhang, Sheng, Song, Guo and Guo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yu Guo, guoyu308@163.com; Yangyang Guo, 984054863@qq.com