Abstract
Background:
As one of the most common cancer, colorectal cancer (CRC) is with high morbidity and mortality. Peritoneal metastasis (PM) is a fatal state of CRC, and few patients may benefit from traditional therapies. There is a complex interaction between PM and immune cell infiltration. Therefore, we aimed to determine biomarkers associated with colorectal cancer peritoneal metastasis (CRCPM) and their relationship with immune cell infiltration.
Methods:
By informatic analysis, differently expressed genes (DEGs) were selected and hub genes were screened out. RAB13, one of the hub genes, was identificated from public databases and validated in CRC tissues. The ESTIMATE, CEBERSORT and TIMER algorithms were applied to analyze the correlation between RAB13 and immune infiltration in CRC. RAB13’s expression in different cells were analyzed at the single-cell level in scRNA-Seq. The Gene Set Enrichment Analysis (GSEA) was performed for RAB13 enrichment and further confirmed. Using oncoPredict algorithm, RAB13’s impact on drug sensitivity was evaluated.
Results:
High RAB13 expression was identified in public databases and led to a poor prognosis. RAB13 was found to be positively correlated with the macrophages and other immune cells infiltration and from scRNA-Seq, RAB13 was found to be located in CRC cells and macrophages. GSEA revealed that high RAB13 expression enriched in a various of biological signaling, and oncoPredict algorithm showed that RAB13 expression was correlated with paclitaxel sensitivity.
Conclusion:
Our study indicated clinical role of RAB13 in CRC-PM, suggesting its potential as a therapeutic target in the future.
1 Introduction
Colorectal cancer (CRC) is with about 2 million new cases and 1 million deaths from the disease worldwide every year (), making it one of the most common malignancies. Most of CRCs are poorly differentiated and with distant metastasis, and for unresectable patients, more than 1/3 of the patients died of disease progression within 5 years (). Surgery, chemoradiotherapy, targeted therapy and immuno-therapy are the conventional treatment for CRC, but the effect is limited and the side effects are huge.
Patients with advanced CRC may experience metastasis to different organs, including liver, peritoneum and lung (). Patients with CRC-PM often mean that they cannot achieve no evidence of disease (NED) through local treatments such as surgery, which is different from liver or lung metastasis. PM is discovered in approximately 13% of CRC cases, most of which are peritoneal and combined with other organs (). The average survival of CRC-PM patients was shorter than patients with liver-isolated metastasis. Therefore, it is important to explore the pathogenesis of CRC-PM, and to identify new diagnostic markers.
The Ras, Rho, Ran, Arf, and Rab families constitute the small GTPases of the Ras superfamily, which are core regulators of intracellular trafficking signals as scaffold proteins and regulate molecules and downstream effectors (, ). One such Rab GTPases, namely RAB13, controls cellular functions that are often altered in cancer, and found to be up-regulated in different kinds of cancers. For example, in ovarian cancer, RAB13 was found as a marker of epithelial-mesenchymal transition (EMT) and could promote metastasis by regulating cytoskeleton and tight junctions (). However, the role of RAB13 in CRC is still not clear.
As a branch of artificial intelligence, machine learning (ML) has been applied in different fields to analyze massive amounts of data. Several studies described the applicability of ML models to cancer studies (). ML combines large amount of variables by a computer algorithm, while traditional statistical methods only use selected variables for further calculations, so ML improves the forecasting accuracy (). Here, we combined bulk RNA-seq with scRNA-seq to investigate the role of RAB13 in CRC-PM and immune cell infiltration.
In this study, we detected abnormal RAB13 expression in CRC-PM. We also demonstrated RAB13 promoted CRC-PM by promoting macrophage infiltration. Therefore, RAB13 is a promising target for CRC-PM immunotherapy.
2 Methods and materials
2.1 Public data collection
From the National Center of Biotechnology Information Gene Expression Omnibus (NCBI GEO) database, we searched raw gene expression data of the CRC-PM: colorectal AND (cancer OR tumor OR carcinoma) AND peritoneal metastasis. Altogether, 94 datasets were brought into initial search. Mainly inclusion criteria included: (1) The dataset included gene expression profiling data. (2) colorectal cancer samples with or without peritoneal metastasis were involved and at least 3 samples were included each group. By this criteria, two GEO datasets, GSE21510 and GSE87211 were obtained for investigating colorectal cancer peritoneal metastasis hub genes (CRC-PM-Hub Gene).
GEO2R () was used to screen DEGs (differentially expressed genes) between samples with or without peritoneal metastasis in the GEO datasets (the cut-off criteria: P-value < 0.05 and absolute log fold-change (FC) > 2) and visualized by volcano plot. Besides, RNA-seq data from TCGA-COAD were also downloaded and analyzed. The RobustRankAggreg (RRA) method was used to integrate and analyze the three datasets (GSE21510 and GSE87211 and TCGA-COAD) to obtain the common DEGs (). R package ‘VennDiagram’ was used to identify the overlapping genes between the common DEGs and peritoneal metastasis genes obtained from GeneCards database ().
Gene Ontology (GO) is a standardized language and classification system used to describe gene function and genomic data (), while Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource containing functions and information of biological systems (). GO and KEGG analysis of the DEGs were performed using WebGestalt ().
2.2 The identification of CRC-PM-hub gene
Using the glmnet package of R (), Least Absolute Shrinkage and Selection Operator (LASSO) analysis was enforced based on the DEGs. LASSO algorithm are instrumental in mitigating the risk of overfitting and enhancing the accuracy of predictions (). The model or gene signature was defined as: colorectal cancer peritoneal metastasis risk score (CRC-PM RS) = ∑Exp(i)*coef(i) (i: genes; exp: expression levels; coef: LASSO coefficients). The reliability and stability of the model were proven in the GSE17538 and GSE37892 datasets.
We further validated CRC-PM genes from the LASSO coefficients by two different algorithms in GSE17538 dataset. Random forest is based on multiple decision trees (ensemble), which employed bootstrap aggregation and predictor randomization to attain notable predictive accuracy (). Support Vector Machine-Recursive Feature Elimination (SVM-RFE) trains a subset of features to reduce the feature set and find the most predictive features (). In SVM-RFE, the feature genes were the genes corresponding to the minimum point of cross-validation error.
2.3 Expression, mutation, methylation and phosphorylation of RAB13 in TCGA-COAD
(1) Expression: TCGA data was collected to analyze the RAB13 mRNA expression level and the survival package was used to estimate the survival rate of TCGA_COAD patients with different RAB13 expression.
(2) Mutation: R package ‘maftools’ was used for the mutation analysis to identify variations of mutated genes in RAB13 high and low expression groups (, ).
(3) Methylation: We assessed the survival time of patients with CpG methylation in RAB13 based on the MethSurv database, a web tool to investigate correlations between survival time and CpG methylation patterns ().
(4) Phosphorylation: The phosphorylation levels of RAB13 were analyzed with the CPTAC database, a central repository for the public dissemination of proteomic sequence dataset ().
2.4 Patients tissue samples and immunohistochemistry
The research was agreed by Ethics Committee of Liaoning Provincial People’s Hospital. Informed consent were signed and complianced in sample collection in all involved patients. A total of 92 CRC tissues were obtained from surgical resection.
Hematoxylin and eosin (H&E) staining and immunohistochemical staining were done to assess RAB13 expression. Immunoreactivity was semi-quantitatively as we described before ().
2.5 Evaluation of immune cell type components
The Estimation of Stromal and Immune Cells (ESTIMATE) algorithm was performed to calculate the ImmuneScore (proportion of immune ingredient) by R package ‘estimate’ (). By using CIBERSORT (), we investigated the correlation of RAB13 expression with diverse immune infiltrating cell types. Single-Sample Gene Set Enrichment Analysis (ssGSEA) is to investigate the correlation of tumor-infiltrating immune cells with their associated functional pathways (). We performed ssGSEA analyze immune signatures using the GSVA package. Immune infiltration from TCGA-COAD was analyzed in Tumor Immune Estimation Resource database (TIMER 2.0).
2.6 Analysis of RAB13 in CRC using scRNA-seq data
Single-cell RNA sequencing (scRNA-seq) datasets (GSE136394, GSE161277 and GSE188741) were downloaded from GEO and we analyzed the scRNA-seq data using Seurat R package (). Cell type annotation was performed by CellMarker database and RAB13’s expression in different cells was visualized by violinplot. Lastly, the R package ‘CellChat’ () was employed to analyze cell-to-cell communication networks and interactions.
2.7 Gene set enrichment analysis
GSEA is conducted by using GSEA v2.2.2 to identify RAB13 associated gene sets in GSE17538 and GSE37892 ().
2.8 Statistical methods
We analyzed the experimental data by SPSS software (Chicago, USA), who was presented as the mean ± SD. Differences were considered significant when p-values < 0.05.
3 Results
3.1 Identification of CRC-PM-DEGs in public database
Following careful screening of the content, the two datasets with the largest sample size (GSE21510 and GSE87211) were obtained. The GSE21510 dataset consisted of 123 CRC and 23 normal samples () and the GSE87211 dataset contained 230 CRC and 133 normal samples (). Besides, we also downloaded gene expression records of 619 CRC patients from TCGA-COAD. We identified 157, 78 and 113 upregulated genes in the GSE21510, GSE87211 and TCGA-COAD datasets, respectively, and visualized them in a volcano plot (Figure 1A). The overlapping DEGs were screened by the RRA method and 161 integrated DEGs, comprising 74 upregulated genes and 87 downregulated genes, were identified, as shown in Supplementary Table 1 and Figure 1B. Next, we intersected the 74 upregulated genes with PM genes obtained from the GeneCards database and obtained 30 genes (Figure 1C). The GO and KEGG analysis results showed that DEGs may have multiple biological functions: the collagen fibril organization was the most significant biological process (BP) term, extracellular space was the cellular component (CC) term with highest significance, and the most significant (MF) term was metalloendopeptidase activity (Figures 2A–C). The KEGG pathway showed that the DEGs were involved in TNF signaling pathways (Figure 2D).
Figure 1
Figure 2
3.2 Establishment and validation of a signature based on PM-DEGs
To further validate the prognostic value of CRC-PM-DEGs, LASSO regression was conducted to construct a signature that predict the prognosis of patients in the training cohort, the GSE17538 datasets. Six genes (THBS2, ANXA2, RAB13, SYK, USP11, PITX2) were finally identified to establish the signature (Figure 3A). The risk coefficients were calculated from the correlation coefficients of the six hub CRC-PM-DEGs with the following formula (Figure 3B): CRC-PM-RS = THBS2*0.27+ANXA2*0.22+RAB13*0.18+SYK*0.06-USP11**0.5-PITX2*0.02.
Figure 3
Based on the CRC-PM-RS, the patients in the GSE17538 dataset were divided into the high- and low-risk subgroups (Figure 3C). Survival status scatterplot and Kaplain-Meier survival analyses suggested that patients with high-Risk score had poorer survival than those with low-Risk score (Figure 3D). To validate the model, we used the GSE37892 dataset as an external validation cohort. The analysis showed a statistically significant difference in overall survival between the low- and high-risk groups (Figures 3E–F).
We next integrated another two different algorithms in GSE17538 dataset to increase diagnostic effectiveness while reducing noise information. By random forest algorithm, the most valuable top 2 genes in CRC-PM was CLDN1 and RAB13 (Figure 3G). In addition, the SVM-RFE results show that the best prediction performance can be obtained by selecting 7 feature variables (RAB13, CLDN1, ANXA2, SYK, TPM3, IL1A, USP11, Figure 3H). In the CRC-PM-RS formula, RAB13 has a relative high coefficient, which also rank at the front in random forest and SVM-RFE algorithms, so we defined RAB13 as the CRC-PM hub gene for further analysis.
3.3 Expression, prognosis, mutation, methylation and phosphorylation of RAB13 in TCGA-COAD
To further understanding the biological function of RAB13, we employed TCGA-COAD to evaluate RAB13 expression levels. It revealed that RAB13 mRNA level was significantly elevated in CRC tissues (Figure 4A). Kaplan-Meier analysis showed that high expression of RAB13 has a tendency to prolong survival, but there is no statistical significance (p = 0.09, Figure 4B). Expression of RAB13 was also correlated with T stage in TCGA-COAD (Figure 4C). CRC is a kind of severe cancer with high heterogeneity and genetic factors (), so we explored the association between RAB13 expression and somatic mutations in TCGA-COAD by maftools in R. The results showed that somatic mutations differed between RAB13 high and low groups. For example, as shown in Figure 4D, 52% of TP53 mutation, which was with the highest mutation frequency, was observed in RAB13 low group, while 46% in RAB13 high group. Pathogenic mutations, such as CTNNB1 () are more common in the RAB13 high expression group (Figure 4E). Treatment approaches discovered driver events in colorectal cancer, such as mutations in EGFR, led to profound improvements in clinical outcomes (). The difference of the driver events between RAB13 high and low groups was also explored. As shown in Figure 4F, EGFR mutation was higher in RAB13 low group. That means targeted therapy may be more effective in patients with low RAB13 expression.
Figure 4
DNA methylation is an important epigenetic characteristic, and DNA methylation levels of RAB13 with the prognostic value of each single CpG were investigated using the MethSurv tool. One CpG site’s methylation level, namely cg26336059, was strongly associated with prognosis (Figure 4G). In the CPTAC dataset, we also found that the phosphorylation levels of S111 and S178 of RAB13 CRC were increased (Figures 4H, I).
To further understanding the expression of the CRC-PM hub gene, RAB13, the expression of RAB13 was detected in 92 CRC tissues by immunohistochemistry. The expression of RAB13 was upregulated in CRC tissues for 64.1% of the patients (59/92), and RAB13 was mostly expressed in the cytoplasm (Supplementary Figures 1A, B).
3.4 RAB13 correlated with tumor immunity in CRC
We next evaluated the association between RAB13 expression and tumor immunity using ESTIMATE algorithm. As shown in Figures 5A, B, the RAB13 expression was significantly positively correlated with Immune Score (GSE21510: Pearson’s r = 0.15, P = 0.03, GSE87211: Pearson’s r = 0.26, P = 5.29e-3) and ESTIMATE Score (GSE21510: Pearson’s r = 0.16, P = 0.02, GSE87211: Pearson’s r = 0.24, P = 7.87e-3). The R-value of correlation analysis was among -0.3~0.3, which was marginally statistical significant, so we employed different algorithms to confirm the conclusion. First, we calculated tumor-infiltrating immune cell proportions in the two GEO datasets employing CIBERSORT algorithm. In GSE21510 dataset, M2 macrophages and resting NK cells were decreased in RAB13 high group (Figure 5C), while in GSE87211 datasets, memory B cells and mast cells were changed in RAB13 high group (Figure 5D). Next, the tumor-infiltrating immune cells in CRC was further analyzed by ssGSEA in the two GEO datasets, which showed significantly positive correlations between RAB13 expression and Monocyto-Macropahges or Mast cell (Figures 5E, F). At last, the TIMER database also showed that RAB13 expression was correlated with different immune cells infiltration in TCGA-COAD datasets (Figure 5G).
Figure 5
Next, we explored RAB13 expression at single level by scRNA-seq data in GSE136394 (comprising 6 types of cells and 13 clusters), GSE161277 (comprising 6 types of cells and 11 clusters) and GSE188741 (comprising 7 types of cells and 11 clusters) datasets. Single-cell subgroup analysis isolated different cell populations and RAB13 high-expression cell populations were identified. As shown in Figures 6A–C and Supplementary Figure 2, RAB13 expression was exclusively prominent in malignant cells and macrophages, indicating that RAB13 may play a crucial role in macrophage infiltration, which warrants further investigation. Finally, we assessed interactions between different cell subpopulations in CRC patients by CellChat and close intercellular communication was observed between macrophage and CRC cells (Figures 6D, E).
Figure 6
3.5 GSEA and prediction of potentially sensitive drugs
In order to obtain understanding into the biological pathways and processes affected by RAB13, GSEA was conducted in GSE21510 and GSE87211 datasets (Figures 7A, B). The results showed that COMPLRMENT AND COAGULATION DESCADES was enriched in RAB13 high expression group in the two datasets collectively (GSE21510: NES = -0.60, adjust P = 1.38E-06, Figure 7C left; GSE87211: NES = -0.52, adjust P = 1.08E-03, Figure 7C right). Coagulant substances such as von Willebrand factor (vWF) are stored in Weibel Palade bodies (WPB), a specialized secretory granules, in a ready-to-be-used form and the biogenesis of WPB are regulated by RAB13 (). Besides, in GSE21510 datasets, CELL CYCLE was enriched in RAB13 high expression group (Figure 7D).
Figure 7
Next, we employed the R package ‘oncoPredict’ to explore the common chemotherapy drugs sensitivity with different RAB13 expression. The results showed that the IC50 values of paclitaxel was positively correlated with RAB13 expression (GSE21510: Pearson’s r = 0.14, P = 0.04, Figure 7E left; GSE87211: Pearson’s r = 0.20, P = 0.03, Figure 7E right). These results suggested that RAB13 impacted paclitaxel sensitivity in CRC.
4 Discussion
Most of CRC patients is died of metastasis, which is developed in nearly half of CRC patients when diagnosed and leading to a poor prognosis (). Patients have a heterogeneous prognosis, as there are great differences in metastasis patterns and tumor characteristics among CRC (). Therefore, it is necessary to identify CRC patients with high metastatic risk for treatment decision-making, long-term survival assessment and follow-up frequency. In this study, by employing public datasets, we filtered and identified 74 CRC-PM DEGs, and GO and KEGG analyses were performed on these 74 DEGs.
It is indispensable to find new diagnostic biomarkers and therapeutic targets to improve the prognosis and quality of life of CRC patients. Comprehensive studies have suggested that abnormal gene expression and gene variation may be involved in the colorectal carciongenesis, however, the molecular pathogenesis of CRC is not fully understood. Rab proteins, members of the RAS GTPase superfamily, had crucial roles in the membrane trafficking and fusion events (). RAB13’s expression has been reported to be aberrant in many cancers, including hepatocellular carcinoma () and gastric cancer (), but there is rare study on the RAB13 expression and clinical significance in CRC patients. Our study found that RAB13 was up-regulated in CRC and high expression of RAB13 indicated shorter progress free survival.
At the genetic level, CRC is a clonal process where somatic gene mutations accumulate to result in abnormal growth of normal intestinal epithelial cells. Tens of drugs targeting receptor tyrosine kinase (RTK) in CRC patients have been uncovered. In this study, we found that druggable mutations, such as EGFR, had higher mutation frequencies in RAB13-low expression group. That indicated CRC patients with low RAB13 expression may get benefit from targeted therapy.
Tumor-infiltrating immune cells were closely related to the tumor biological behavior and prognosis (). However, the role of tumor-infiltrating immune cells in CRC-PM is still not clear. In our study, by different algorithms, we found RAB13 expression was positively correlated with macrophages, and other immune cells. Zhu et al’s study in breast cancer also showed that RAB13 controlled the C-X-C chemokine receptor type 1/2 (CXCR1/2) membrane translocation to make tumor cells to interact with tumor-associated macrophages (), indicating RAB13 mediated interaction of tumor cells and macrophages in different cancers. The scRNA-seq technologies dissected the heterogeneity of tumor-infiltrating immune cells and their interactions with the diverse cell types in tumors (). By scRNA analysis in GEO datasets, we found that RAB13 was observed in tumor cells and macrophages, and cell chat analysis showed that osteoclasts interacted with tumor cells in CRC.
In this study, we showed RAB13’s role and its potential value in CRC-PM, nevertheless, some limitations existed in our research. Firstly, the mechanisms of RAB13 in CRC-PM are necessary to be validated. Secondly, we only used several GEO cohorts for verification, which may have some systematic biases. At last, CRC-PM is a complex process but we only considered the association between RAB13 and CRC-PM, and not consider the relationship between genes and genes or environment (). We will implement these limitations in future studies.
5 Conclusion
In conclusion, this study provided a better understanding of RAB13 in CRC-PM and for further investigation of CRC-PM diagnosis and treatment.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: NCBI GEO database (https://www.ncbi.nlm.nih.gov/gds/); GeneCards database (https://www.genecards.org/); TCGA (https://www.cancer.gov/ccg/research/genome-sequencing/tcga); MethSurv database (https://biit.cs.ut.ee/methsurv); CPTAC database (https://proteomics.cancer.gov/programs/cptac).
Ethics statement
The studies involving humans were approved by Ethics Committee of Liaoning Provincial People’s Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
FL: Writing – original draft, Writing – review & editing, Investigation, Methodology, Data curation, Validation. XL: Data curation, Writing – original draft, Writing – review & editing, Conceptualization, Formal analysis, Investigation, Visualization. ZW: Validation, Writing – review & editing, Data curation. XW: Investigation, Validation, Writing – review & editing, Project administration, Data curation. JL: Conceptualization, Data curation, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
The authors acknowledge Zhenguang Du and Mingming Xiao, both from People’s Hospital of Liaoning Province, for their help in tissue sample collection and immunohistochemistry.
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/fimmu.2024.1403008/full#supplementary-material
References
1
SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660
2
YeoHBetelDAbelsonJSZhengXEYantissRShahMA. Early-onset colorectal cancer is distinct from traditional colorectal cancer. Clin Colorectal Cancer. (2017) 16:293–9. doi: 10.1016/j.clcc.2017.06.002
3
LuoDLiuQYuWMaYZhuJLianPet al. Prognostic value of distant metastasis sites and surgery in stage IV colorectal cancer: a population-based study. Int J Colorectal Dis. (2018) 33:1241–9. doi: 10.1007/s00384-018-3091-x
4
FrankoJShiQMeyersJPMaughanTSAdamsRASeymourMTet al. Prognosis of patients with peritoneal metastatic colorectal cancer given systemic therapy: an analysis of individual patient data from prospective randomised trials from the Analysis and Research in Cancers of the Digestive System (ARCAD) database. Lancet Oncol. (2016) 17:1709–19. doi: 10.1016/S1470-2045(16)30500-9
5
CherfilsJZeghoufM. Regulation of small GTPases by GEFs, GAPs, and GDIs. Physiol Rev. (2013) 93:269–309. doi: 10.1152/physrev.00003.2012
6
SchwartzSLCaoCPylypenkoORakAWandinger-NessA. Rab GTPases at a glance. J Cell Sci. (2007) 120:3905–10. doi: 10.1242/jcs.015909
7
GuoJHanXLiJLiZYiJGaoYet al. Single-cell transcriptomics in ovarian cancer identify a metastasis-associated cell cluster overexpressed RAB13. J Transl Med. (2023) 21:254. doi: 10.1186/s12967-023-04094-7
8
HamamotoRKomatsuMTakasawaKAsadaKKanekoS. Epigenetics analysis and integrated analysis of multiomics data, including epigenetic data, using artificial intelligence in the era of precision medicine. Biomolecules. (2019) 10:62. doi: 10.3390/biom10010062
9
LiYMLiZLChenFLiuQPengYChenM. A LASSO-derived risk model for long-term mortality in Chinese patients with acute coronary syndrome. J Transl Med. (2020) 18:157. doi: 10.1186/s12967-020-02319-7
10
DibounIWernischLOrengoCAKoltzenburgM. Microarray analysis after RNA amplification can detect pronounced differences in gene expression using limma. BMC Genomics. (2006) 7:252. doi: 10.1186/1471-2164-7-252
11
KoldeRLaurSAdlerPViloJ. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. (2012) 28:573–80. doi: 10.1093/bioinformatics/btr709
12
LinYHuZ. Bioinformatics analysis of candidate genes involved in ethanol-induced microtia pathogenesis based on a human genome database: GeneCards. Int J Pediatr Otorhinolaryngol. (2021) 142:110595. doi: 10.1016/j.ijporl.2020.110595
13
AshburnerMBallCABlakeJABotsteinDButlerHCherryJMet al. Gene ontology: tool for the unification of biology. Gene Ontol Consortium Nat Genet. (2000) 25:25–9. doi: 10.1038/75556
14
KanehisaMGotoS. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
15
ZhangBKirovSSnoddyJ. WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. (2005) 33:W741–8. doi: 10.1093/nar/gki475
16
EngebretsenSBohlinJ. Statistical predictions with glmnet. Clin Epigenet. (2019) 11:123. doi: 10.1186/s13148-019-0730-1
17
RafiqueRIslamSKaziJU. Machine learning in the prediction of cancer therapy. Comput Struct Biotechnol J. (2021) 19:4003–17. doi: 10.1016/j.csbj.2021.07.003
18
IshwaranHKogalurUB. Consistency of random survival forests. Stat Probab Lett. (2010) 80:1056–64. doi: 10.1016/j.spl.2010.02.020
19
HuangMLHungYHLeeWMLiRKJiangBR. SVM-RFE based feature selection and Taguchi parameters optimization for multiclass SVM classifier. Sci World J. (2014) 2014:795624. doi: 10.1155/2014/795624
20
NgPKLiJJeongKJShaoSChenHTsangYHet al. Systematic functional annotation of somatic mutations in cancer. Cancer Cell. (2018) 33:450–62. doi: 10.1016/j.ccell.2018.01.021
21
MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. (2018) 28:1747–56. doi: 10.1101/gr.239244.118
22
ModhukurVIljasenkoTMetsaluTLokkKLaisk-PodarTViloJ. MethSurv: a web tool to perform multivariable survival analysis using DNA methylation data. Epigenomics-uk. (2018) 10:277–88. doi: 10.2217/epi-2017-0118
23
WhiteakerJRHalusaGNHoofnagleANSharmaVMacLeanBYanPet al. CPTAC Assay Portal: a repository of targeted proteomic assays. Nat Methods. (2014) 11:703–4. doi: 10.1038/nmeth.3002
24
LvFLiXWangYHaoL. MAGP1 maintains tumorigenicity and angiogenesis of laryngeal cancer by activating Wnt/beta-catenin/MMP7 pathway. Carcinogenesis. (2024) 45:220–34. doi: 10.1093/carcin/bgad003
25
YoshiharaKShahmoradgoliMMartinezEVegesnaRKimHTorres-GarciaWet al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. (2013) 4:2612. doi: 10.1038/ncomms3612
26
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
27
SubramanianATamayoPMoothaVKMukherjeeSEbertBLGilletteMAet al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
28
GribovASillMLuckSRuckerFDohnerKBullingerLet al. SEURAT: visual analytics for the integrated analysis of microarray data. BMC Med Genomics. (2010) 3:21. doi: 10.1186/1755-8794-3-21
29
JinSGuerrero-JuarezCFZhangLChangIRamosRKuanCHet al. Inference and analysis of cell-cell communication using CellChat. Nat Commun. (2021) 12:1088. doi: 10.1038/s41467-021-21246-9
30
TsukamotoSIshikawaTIidaSIshiguroMMogushiKMizushimaHet al. Clinical significance of osteoprotegerin expression in human colorectal cancer. Clin Cancer Res. (2011) 17:2444–50. doi: 10.1158/1078-0432.CCR-10-2884
31
HuYGaedckeJEmonsGBeissbarthTGradeMJoPet al. Colorectal cancer susceptibility loci as predictive markers of rectal cancer prognosis after surgery. Genes Chromosomes Cancer. (2018) 57:140–9. doi: 10.1002/gcc.22512
32
BaiZZhouYYeZXiongJLanHWangF. Tumor-infiltrating lymphocytes in colorectal cancer: the fundamental indication and application on immunotherapy. Front Immunol. (2021) 12:808964. doi: 10.3389/fimmu.2021.808964
33
KikuchiAKishidaSYamamotoH. Regulation of Wnt signaling by protein-protein interaction and post-translational modifications. Exp Mol Med. (2006) 38:1–10. doi: 10.1038/emm.2006.1
34
YuanMHuangLLChenJHWuJXuQ. The emerging treatment landscape of targeted therapy in non-small-cell lung cancer. Signal Transduct Target Ther. (2019) 4:61. doi: 10.1038/s41392-019-0099-9
35
NassJTerglaneJGerkeV. Weibel palade bodies: unique secretory organelles of endothelial cells that control blood vessel homeostasis. Front Cell Dev Biol. (2021) 9:813995. doi: 10.3389/fcell.2021.813995
36
AbedizadehRMajidiFKhorasaniHRAbediHSabourD. Colorectal cancer: a comprehensive review of carcinogenesis, diagnosis, and novel strategies for classified treatments. Cancer Metastasis Rev. (2023) 43:729. doi: 10.1007/s10555-023-10158-3
37
KnudsenMDWangKWangLPolychronidisGBerstadPWuKet al. Development and validation of a risk prediction model for post-polypectomy colorectal cancer in the USA: a prospective cohort study. EClinicalMedicine. (2023) 62:102139. doi: 10.1016/j.eclinm.2023.102139
38
StenmarkH. Rab GTPases as coordinators of vesicle traffic. Nat Rev Mol Cell Biol. (2009) 10:513–25. doi: 10.1038/nrm2728
39
JiangCLiuZYuanJWuZKongLYangJet al. Construction of two independent RAB family-based scoring systems based on machine learning algorithms and definition of RAB13 as a novel therapeutic target for hepatocellular carcinoma. Int J Mol Sci. (2023) 24:4335. doi: 10.3390/ijms24054335
40
ChenPChenGWangCMaoC. RAB13 as a novel prognosis marker promotes proliferation and chemotherapeutic resistance in gastric cancer. Biochem Biophys Res Commun. (2019) 519:113–20. doi: 10.1016/j.bbrc.2019.08.141
41
FridmanWHPagesFSautes-FridmanCGalonJ. The immune contexture in human tumours: impact on clinical outcome. Nat Rev Cancer. (2012) 12:298–306. doi: 10.1038/nrc3245
42
WangHXuHChenWChengMZouLYangQet al. Rab13 sustains breast cancer stem cells by supporting tumor-stroma cross-talk. Cancer Res. (2022) 82:2124–40. doi: 10.1158/0008-5472.CAN-21-4097
43
RenXZhangLZhangYLiZSiemersNZhangZ. Insights gained from single-cell analysis of immune cells in the tumor microenvironment. Annu Rev Immunol. (2021) 39:583–609. doi: 10.1146/annurev-immunol-110519-071134
44
RischNMerikangasK. The future of genetic studies of complex human diseases. Science. (1996) 273:1516–7. doi: 10.1126/science.273.5281.1516
Summary
Keywords
colorectal cancer, peritoneal metastasis, Rab13, macrophage, scRNA-seq
Citation
Lv F, Li X, Wang Z, Wang X and Liu J (2024) Identification and validation of Rab GTPases RAB13 as biomarkers for peritoneal metastasis and immune cell infiltration in colorectal cancer patients. Front. Immunol. 15:1403008. doi: 10.3389/fimmu.2024.1403008
Received
18 March 2024
Accepted
29 July 2024
Published
13 August 2024
Volume
15 - 2024
Edited by
Jieqiong Liu, Sun Yat-sen University, China
Reviewed by
Chen Xiaofeng, Nanjing Medical University, China
Qin Liu, Nanjing University, China
Updates
Copyright
© 2024 Lv, Li, Wang, Wang and Liu.
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: Jing Liu, liujing23@vip.163.com; Xiaobo Wang, 121737072@qq.com
†These authors have contributed equally to this work and share first authorship
Disclaimer
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.