- Open Access
N6-methyladenosine (m6A) RNA methylation regulator SNRPC is a prognostic biomarker and is correlated with immunotherapy in hepatocellular carcinoma
World Journal of Surgical Oncology volume 19, Article number: 241 (2021)
Hepatocellular carcinoma (HCC) is one of the most common malignancies in the world, and due to its complex pathogenic factors, its prognosis is poor. N6-methyladenosine (m6A) RNA methylation plays an important role in the tumorigenesis, progression, and prognosis of many tumors. The m6A RNA methylation regulator small nuclear ribonucleoprotein polypeptide C (SNRPC), which encodes one of the specific protein components of the U1 small nuclear ribonucleoprotein (snRNP) particle, has been proven to be related to the prognosis of patients with HCC. However, the effect of SNRPC on the tumor microenvironment and immunotherapy in HCC remains unclear.
The HCC RNA-seq profiles in The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases, including 421 LIHC and 440 LIRI-JP samples, respectively, were used in this study. Both the expression of SNRPC in HCC was upregulated in the TCGA and ICGC databases compared to normal tissues. Next, the expression of SNRPC was validated as a risk factor for prognosis by Kaplan-Meier analysis and employed to establish a nomogram with T pathologic stage. By gene set variation (GSVA) analysis and gene set enrichment (GSEA) analysis, we found that SNRPC was mainly related to protein metabolism and the immune process. Furthermore, the estimation of stromal and immune cells in malignant tumor tissues using expression (ESTIMATE), microenvironment cell population counter (MCP-counter), and single sample GSEA (ssGSEA) algorithms revealed that the high-SNRPC group had a lower stromal score, lower abundance of endothelial cells and fibroblasts, and lower immune infiltration. Ultimately, a tumor immune dysfunction and exclusion (TIDE) analysis revealed that patients in the low-SNRPC group may be more sensitive to immune checkpoint inhibitor therapy.
SNRPC could serve as a promising prognostic and immunotherapeutic marker in HCC and might contribute to new directions and strategies for HCC treatment.
Liver cancer was the third leading cause of cancer death worldwide in 2020 . Hepatocellular carcinoma (HCC) is the predominant form of primary liver cancer and accounts for 75–85% of cases. Chronic infection with hepatitis B virus (HBV) or hepatitis C virus (HCV), heavy alcohol intake, type 2 diabetes, and smoking have been related to the occurrence of HCC . More than 1 million liver cancer deaths are estimated by 2030 . Although surgical resection of early HCC is still the preferred treatment, the recurrence rate at 5 years can be as high as 70% . In addition, even if liver transplantation can reduce the recurrence rate, the number of available donor organs is always insufficient relative to the demand . Moreover, the majority (> 70%) of patients with advanced disease do not benefit from surgery . For these patients, immunotherapy is a potential option. In recent years, emerging immunotherapies, such as immune checkpoint blockade and chimeric antigen receptor t cell (CAR-T) immunotherapy, have shown better therapeutic outcomes for some tumors [6, 7]. Nevertheless, responses to immunotherapy typically occur in a minority of patients, ranging from 20 to 50% depending on the tumor type . Thus, it is especially important to identify candidate biomarkers to target patients who have the greatest likelihood of benefiting from immunotherapy . In 2017, the Food and Drug Administration (FDA) approved the PD-1 inhibitor nivolumab (OPDIVO®), which marked the beginning of the immunotherapy era for HCC . Compared with certain tumors, immunotherapy for liver cancer is still in its infancy. Additionally, knowledge of the biomarkers for predicting the responses or primary resistance to immunotherapeutic agents is crucial for establishing immunotherapy criteria for HCC patients in the future.
Genome-wide analysis of mRNA expression profiles has been used to screen prognostic and treatment-related markers in many cancers in recent years [10,11,12,13]. The N6-methyladenosine (m6A) methylation modification occurs at the 6th N atom of adenine, and it is the most common posttranscriptional modification on RNA and mediates > 60% of RNA methylation [14, 15]. Studies have shown that m6A RNA methylation plays an important role in tumor occurrence, development, and prognosis . As an important m6A RNA methylation regulator, small nuclear ribonucleoprotein polypeptide C (SNRPC) encodes one of the specific protein components of the U1 small nuclear ribonucleoprotein (snRNP) particle, which is required for the formation of the spliceosome . A recent study showed that SNRPC has the potential to promote HCC cell motility by inducing epithelial-mesenchymal transition . In addition, SNRPC also contributes to sex bias in systemic autoimmune diseases and plays a role in spinal muscular atrophy pathogenesis [19, 20]. However, to our knowledge, the interaction between SNRPC and immunotherapy in HCC has not been reported. Therefore, in this study, we aimed to analyze the potential of SNRPC as a prognostic biomarker and predictor of immunotherapy responses in HCC to provide reference data for the clinical treatment of HCC.
In this study, we found that SNRPC expression was upregulated in HCC samples compared to normal samples. Next, the expression of SNRPC was validated as a risk factor for prognosis by a Kaplan-Meier analysis and employed to establish a nomogram with T pathologic stage. Furthermore, we investigated the correlation between SNRPC and immunotherapeutic reactions through a tumor immune dysfunction and exclusion (TIDE) analysis. These findings may contribute to new directions and strategies for HCC treatment.
Materials and methods
Public transcriptome profiling
In this study, we used 861 HCC cases from two public databases. After excluding samples without clinical annotation and survival information, 421 LIHC samples (50 normal samples labeled with 11A and 364 HCC samples labeled with 01A) were obtained from The Cancer Genome Atlas (TCGA) database and used as a training set. Additionally, 440 LIRI-JP cases (197 normal and 243 HCC) were obtained from the International Cancer Genome Consortium (ICGC) database and used as a validation set. The detailed clinical data of these samples are listed in Table 1. Moreover, both the expression data from the TCGA and ICGC databases were standardized by the log2 function in R.
Validation of the expression of SNRPC in HCC based on the TCGA and ICGC databases
To further investigate the expression of SNRPC in HCC, the Wilcoxon test was used to verify the significant difference in mRNA expression levels of SNRPC between the tumor and normal tissue samples based on the TCGA and ICGC databases. Moreover, the expression of SNRPC was further examined between 50 HCC tissue samples and 50 paired paracancerous tissue samples from the TCGA database (the ICGC database does not include paired cancer tissue and paracancerous tissue). Finally, the protein expression of SNRPC between HCC tissues and normal tissues was verified in the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/).
Survival analysis of SNRPC
To evaluate whether SNRPC could impact the survival probability of HCC, all of the HCC samples in the TCGA and ICGC databases were divided into a high-SNRPC group and a low-SNRPC group based on the median expression value of SNRPC. Subsequently, survival probabilities were analyzed by a Kaplan-Meier analysis and log-rank test , with a p value < 0.05 considered statistically significant. In addition, time-dependent receiver operating characteristic (ROC) curves were plotted to assess the prognostic prediction performance of SNRPC, and the areas under the curve (AUCs) for overall survival (OS) were calculated using the survival ROC R package .
Construction of a nomogram for prognostic prediction
For better clinical application of SNRPC, we established a nomogram by univariate and multivariate Cox regression analysis to evaluate the probability of 1-, 3-, and 5-year OS for HCC patients via the rms R package  by integrating the clinical factors and SNRPC expression.
To investigate the potential biological function of SNRPC in HCC, all of the differentially expressed genes between the high-SNRPC group and the low-SNRPC group in the TCGA database were ranked based on log2FC values. Then, a gene set enrichment analysis (GSEA) was implemented using the clusterProfiler package  to estimate the pathway activity changes in the above two groups. Similarly, the gene set variation (GSVA) package  was utilized to find the predominant pathways associated with SNRPC. Terms with values of FDR < 0.25 and p < 0.05 were retained.
The Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression (ESTIMATE) algorithm can calculate immune scores and stromal scores via the estimate package in R , allowing us to quantify the immune and stromal components in a tumor. Therefore, the differences in immune scores, stromal scores, and tumor purity between the high-SNRPC group and the low-SNRPC group in the TCGA database were assessed by the Wilcoxon test.
Analysis of the tumor microenvironment
To gain further insights into the composition of the tumor microenvironment between the high-SNRPC group and low-SNRPC group in the TCGA database, a microenvironment cell population counter (MCP-counter) was used to quantify the numbers of immune cells, fibroblasts, and epithelial cells from each HCC sample based on marker genes . Then, the single sample GSEA (ssGSEA) algorithm was used to calculate the immune infiltration of each sample, and significant differences in immune cell numbers were identified by analysis of variance (ANOVA).
Tumor immune dysfunction and exclusion analysis
Immune checkpoint pathways are associated with tumor immune escape. Therefore, immune checkpoint inhibitors can enhance anticancer immunity. Thus, we also employed the tumor immune dysfunction and exclusion (TIDE) algorithm  and subclass mapping  to predict the clinical response to immune checkpoint inhibitors between the high-SNRPC group and the low-SNRPC group. Moreover, we also calculated the correlation between SNRPC expression and TIDE scores.
All statistical tests were performed using R (version 3.5.2) software. Continuous data were analyzed using the Wilcoxon test or ANOVA. A p value of less than 0.05 was considered statistically significant in all analyses.
SNRPC is upregulated in HCC
To further investigate the expression of SNRPC in HCC, SNRPC expression between the HCC and normal samples was analyzed using data from the TCGA and ICGC databases. Interestingly, the results of both the TCGA and ICGC databases suggested that SNRPC expression was upregulated in patients with HCC compared to normal samples (Fig. 1A, B). Moreover, the analysis of paired HCC and paracancerous tissues also revealed that SNRPC expression was upregulated in HCC tissues (Fig. 1C). Finally, a deeper staining level of SNRPC in HCC tissues compared to normal tissues revealed a higher protein expression of SNRPC in HCC (Fig. 1D).
SNRPC is related to the survival of HCC
A Kaplan-Meier survival curve was drawn, and the results showed that in both the TCGA and ICGC databases, the low-SNRPC groups had a high survival probability relative to the high-SNRPC groups (p < 0.05) (Fig. 2A, B). In addition, as shown in the ROC curves, the AUCs of SNRPC were 0.614 and 0.729 in the TCGA and ICGC databases, respectively (Fig. 2C, D).
Construction of a nomogram for predicting HCC prognosis
The univariate and multivariate Cox regression analysis results suggested that the pathological T stage and SNRPC were significantly correlated with the prognosis of HCC (Fig. 2E, F). Therefore, a nomogram including T clinical factors and SNRPC expression values was constructed for predicting the 1-, 3-, and 5-year survival probabilities of HCC using samples from the TCGA database (Fig. 2G). A calibration curve demonstrated that our nomogram could accurately estimate the mortality of HCC (Fig. 2H).
Identification of SNRPC-related functional annotation
To gain insights into the potential functions related to SNRPC, we first identified 64 differentially expressed genes (DEGs) between the HCC samples of the high-SNRPC group and the low-SNRPC group (Fig. 3A). We found that catalytic activity-related biological processes, including catalytic_activity_acting_on_DNA, catalytic_step_2_spliceosome, and catalytic_activity_acting_on_RNA, were mainly activated in the high-SNRPC group (Table 2). Conversely, metabolism-related biological processes, including alcohol metabolism, amino acid metabolism, lipid metabolism, etc., were mainly activated in the low-SNRPC group (Table 2). Among the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, the cell cycle-, spliceosome-, ribosome-, and complement-related pathways were mainly activated in the high-SNRPC group while the metabolism-related pathways were mainly activated in the high-SNRPC group (Table 3). In addition, the GSVA results suggested that a large number of immune-related processes were disordered. For example, antigen processing-related biological processes were mainly upregulated in the high-SNRPC group while immune response-related biological processes were mainly upregulated in the low-SNRPC group (Fig. 3B, Table 4).
In our study, we found that except for the immune scores, the stromal scores, ESTIMATE scores, and tumor purity were significantly different between the high- and low-SNRPC groups (Fig. 4). We also found that a lower stromal score and higher tumor purity were associated with poor survival in HCC patients.
Correlation between SNRPC expression and tumor microenvironment
MCP-counter and ssGSEA analyses were conducted to clarify the abundance of various cells in the tumor microenvironment (TME). Immune cells, such as T cells, CD8+T cells, neutrophils, and stromal cells, such as endothelial cells and fibroblasts, were significantly different between the high- and low-SNRPC groups. Among these cell types, neutrophils, endothelial cells, and fibroblasts were increased in the low-SNRPC group while other types were decreased (Fig. 5A, B). Through the ssGSEA, additional subtypes of immune cells were counted. In the low-SNRPC group, only activated CD4+T cells were significantly decreased while CD56dim natural killer cells, central memory CD8+T cells, effector memory CD8+T cells, eosinophils, immature dendritic cells, macrophages, memory B cells, monocytes, natural killer T cells, natural killer T cells, plasmacytoid dendritic cells, regulatory T cells, type 1 T helper cells, and type 17 T helper cells were all increased (Fig. 5C, D).
SNRPC expression could predict the immunotherapeutic response of HCC
The TIDE score is used to predict the efficacy of anti-PD1 and anti-CTLA4 treatments . As shown in Fig. 6A, B, both the high-SNRPC groups showed higher TIDE scores in the TCGA and ICGC databases. Consistently, both the high-SNRPC groups showed no response to anti-PD1 therapy in the TCGA and ICGC databases by subclass mapping analysis (Fig. 6C, D). Furthermore, both the low-SNRPC groups presented more patients with the response to immune checkpoint inhibitor therapy in the TCGA and ICGC databases (Fig. 6E, F). More importantly, the correlation analysis suggested that the expression of SNRPC was positively correlated with the TIDE scores (Fig. 6G). Thus, these findings revealed that patients in the low-SNRPC group may be more sensitive to immune checkpoint inhibitor therapy.
HCC accounts for 75–85% of liver cancers and is also the third leading cause of cancer-related death globally, which has a poor prognosis . Studies have shown that multiple genes play an important role in HCC progression [30,31,32,33,34]. Currently, immunotherapy has already become a viable option for advanced patients. However, not all patients show complete responses to such immunotherapies . As a result, a reliable biomarker for predicting HCC prognosis and immunotherapy response is urgently needed.
In the present study, we further analyzed the role of SNRPC in HCC. SNRPC encodes one of the specific protein components of U1snRNP, which plays an essential role in the formation of the spliceosome. U1snRNP is reported to suppress the migration and invasion of multiple types of tumor cells . Additionally, genes of the snRNP polypeptide family also exert an important role in the tumorigenesis and prognosis of cancers [37, 38]. For example, the high expression of snRNP polypeptide G, which is another indispensable component in the formation of snRNPs, is positively associated with the occurrence, progression and severity of breast, lung, and colon cancers [39,40,41]. Moreover, the snRNP polypeptide N is highly expressed in colorectal cancer tissues and involved in the progression of colorectal cancer . A recent study demonstrated that SNRPC has the potential to enhance HCC cell motility by promoting epithelial-mesenchymal transition, which might contribute to a poor prognosis . However, few studies have focused on the role of SNRPC in the immunotherapy of tumors. Therefore, this study aimed to reveal the effect of SNRPC on the prognosis and immunotherapy of HCC. In our study, we first found that SNRPC was overexpressed in HCC and that higher expression of SNRPC was associated with a poorer prognosis. Then, we constructed a nomogram containing T clinical factors and SNRPC expression values, which can precisely predict 1-, 3-, and 5-year survival probabilities. To gain detailed insights into the mechanisms of SNRPC, we performed functional analyses and found that SNRPC was mainly related to protein metabolism and the immune process. Interestingly, a TME analysis indicated that natural killer cells, CD8+T cells, macrophages, neutrophils, endothelial cells, and fibroblasts were less infiltrated in tissues with high SNRPC expression. These results implied that both stromal cells and immune cells were abundant, which was beneficial for HCC prognosis. Ultimately, we employed the TIDE analysis method to predict the efficacy of anti-PD1 and anti-CTLA4 treatments. The results showed that the high-SNRPC group failed to respond to immune checkpoint inhibitor (ICI) therapy, especially to a PD-1 inhibitor.
At present, many clinical prognostic biomarkers are available for HCC, such as alpha-fetoprotein (AFP), alpha-fetoprotein lens culinaris agglutin-3 (AFP-L3), and des-γ-carboxy prothrombin (DCP). The relationship between AFP and prognosis has been extensively studied . Increasing AFP values are associated with lower survival and a higher tumor recurrence rate in patients at very early or early stages as well as poor prognosis in patients undergoing advanced HCC . By comparison, overexpressed SNRPC is correlated with low infiltration of immune cells and stromal cells, poor prognosis, and poor immunotherapy response. However, SNRPC expression remains relatively complex to detect compared with AFP, which can be detected in serum. We hope that in the future, SNRPC could be used in combination with other biomarkers for the clinical staging system and prognosis prediction of HCC.
At present, numerous biomarkers have been employed to select candidates for immunotherapy. PD-L1 expression is the most widely adopted predictor of the immunotherapeutic response rate. High PD-L1 expression is associated with an increased response rate and clinical benefit in anti-PD-1/anti-PD-L1 therapy . However, PD-L1 expression is not associated with treatment response in some patients. Additionally, the TIDE score can also directly guide the use of ICIs. Patients with higher TIDE scores have higher antitumor immune escape opportunities, resulting in lower ICI treatment response rates . Moreover, the TIDE score has been shown to have higher accuracy than the PD-L1 expression level in predicting the survival outcome of cancer patients treated with ICI agents [45,46,47]. In the current study, the TIDE analysis revealed that patients in the low-SNRPC group may be more sensitive to immune checkpoint inhibitor therapy. Therefore, we speculated that SNRPC might be associated with the immunotherapeutic response of HCC.
Recently, novel strategies have been identified to enhance the effect of anti-PD-1/PD-L1 by promoting T cell infiltration, such as an anti-TGF-β/PD-L1 bispecific antibody and anti-PD-1 combining anti-CTLA-4 [48, 49]. Combining immunotherapy with anti-VEGF agents also displayed a good therapeutic effect in some HCC cases. Consequently, SNRPC might contribute to improving the effectiveness of anti-PD-1/PD-L1 by regulating the expression of SNRPC.
In this study, we showed that the HCC biomarker SNRPC was overexpressed in tumor tissues and linked to a worse prognosis in patients with HCC. Functional analysis showed that SNRPC was mainly related to protein metabolism and the immune process. In addition, we found that upregulated SNRPC may be associated with a less sensitive response to immune checkpoint inhibitor therapy. Therefore, SNRPC may be a promising prognostic and immunotherapeutic marker for HCC. Nevertheless, we only included anti-PD1 and anti-CTLA4 treatments in our TIDE analysis rather than considering more promising immunotherapy methods such as combination therapies. In addition, we did not choose a proposed cutoff for predicting poor prognosis and immunotherapeutic response. Hence, follow-up experiments and clinical validation are extremely necessary.
Hepatitis B virus
Hepatitis C virus
Chimeric antigen receptor T cell
Food and Drug Administration
- SNRPC :
Small nuclear ribonucleoprotein polypeptide C
Small nuclear ribonucleoprotein
Tumor immune dysfunction and exclusion
The Cancer Genome Atlas
International Cancer Genome Consortium
Human Protein Atlas
Receiver operating characteristic
Areas under the curve
Gene set enrichment analysis
Gene set variation
Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression
Microenvironment Cell Population counter
Single sample GSEA
Analysis of variance
Differential expression genes
Kyoto Encyclopedia of Genes and Genomes
Alpha-fetoprotein lens culinaris agglutin-3
Immune checkpoint inhibitors
Sung H, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660
Karb, D.B. and S.N. Sclair, Hepatocellular carcinoma, in liver disease. 2019. 141-154. Springer International Publishing, Cham.
Llovet JM, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Vivarelli M, et al. Liver transplantation for hepatocellular carcinoma under calcineurin inhibitors: reassessment of risk factors for tumor recurrence. Ann Surg. 2008;248(5):857–62.
Thomas MB, et al. Hepatocellular carcinoma: consensus recommendations of the National Cancer Institute Clinical Trials Planning Meeting. J Clin Oncol. 2010;28(25):3994–4005.
Fournier C, et al. Trial Watch: adoptively transferred cells for anticancer immunotherapy. Oncoimmunology. 2017;6(11):e1363139.
Smyth MJ, et al. Combination cancer immunotherapies tailored to the tumour microenvironment. Nat Rev Clin Oncol. 2016;13(3):143–58.
Llovet JM, et al. Molecular therapies and precision medicine for hepatocellular carcinoma. Nat Rev Clin Oncol. 2018;15(10):599–616.
Johnston MP, Khakoo SI. Immunotherapy for hepatocellular carcinoma: current and future. World J Gastroenterol. 2019;25(24):2977–89.
Boyault S, et al. Transcriptome classification of HCC is related to gene alterations and to new therapeutic targets. Hepatology. 2007;45(1):42–52.
Chiang DY, et al. Focal gains of VEGFA and molecular classification of hepatocellular carcinoma. Cancer Res. 2008;68(16):6779–88.
Hoshida Y, et al. Integrative transcriptome analysis reveals common molecular subclasses of human hepatocellular carcinoma. Cancer Res. 2009;69(18):7385–92.
Lee JS, et al. Classification and prediction of survival in hepatocellular carcinoma by gene expression profiling. Hepatology. 2004;40(3):667–76.
Du K, et al. m(6)A RNA methylation controls neural development and is involved in human diseases. Mol Neurobiol. 2019;56(3):1596–606.
Wang S, et al. Roles of RNA methylation by means of N(6)-methyladenosine (m(6)A) in human cancers. Cancer Lett. 2017;408:112–20.
Chen J, et al. Identification of a m(6)A RNA methylation regulators-based signature for predicting the prognosis of clear cell renal carcinoma. Cancer Cell Int. 2020;20:157.
Libri D, et al. A role for the Psi-U mismatch in the recognition of the 5' splice site of yeast introns by the U1 small nuclear ribonucleoprotein particle. J Biol Chem. 2002;277(20):18173–81.
Zhang Y, et al. SNRPC promotes hepatocellular carcinoma cell motility by inducing epithelial-mesenchymal transition. FEBS Open Bio. 2021;11(6):1757–70. https://doi.org/10.1002/2211-5463.13175
Linden M, et al. Sex influences eQTL effects of SLE and Sjogren's syndrome-associated genetic polymorphisms. Biol Sex Differ. 2017;8(1):34.
Cheng D, et al. The arginine methyltransferase CARM1 regulates the coupling of transcription and mRNA processing. Mol Cell. 2007;25(1):71–83.
Heagerty PJ, Saha-Chaudhuri P, Saha-Chaudhuri MP. Package ‘survivalROC’ [J]; 2013.
Heagerty PJ, Zheng Y. Survival model predictive accuracy and ROC curves. Biometrics. 2005;61(1):92–105.
Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. 2008;26(8):1364–70. https://doi.org/10.1200/JCO.2007.12.9791
Yu G, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Hänzelmann SH, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-Seq data. BMC Bioinformatics. 2013;14:7. Published 2013 Jan 16. https://doi.org/10.1186/1471-2105-14-7
Yoshihara K, et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612. https://doi.org/10.1038/ncomms3612
Becht E, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016;17(1):218.
Jiang P, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. 2018;24(10):1550–8.
Hoshida Y, et al. Subclass mapping: identifying common subtypes in independent disease data sets. PLoS One. 2007;2(11):e1195.
Gao S, et al. Reduced expression of DACT2 promotes hepatocellular carcinoma progression: involvement of methylation-mediated gene silencing. World J Surg Oncol. 2013;11:57.
Wu Y, et al. Clinical significance of mitofusin-2 and its signaling pathways in hepatocellular carcinoma. World J Surg Oncol. 2016;14(1):179.
Yu M, et al. Identification and validation of three core genes in p53 signaling pathway in hepatitis B virus-related hepatocellular carcinoma. World J Surg Oncol. 2021;19(1):66.
Luo W, et al. The role of kinesin KIF18A in the invasion and metastasis of hepatocellular carcinoma. World J Surg Oncol. 2018;16(1):36.
Chen S, et al. Low expression of organic anion-transporting polypeptide 1B3 predicts a poor prognosis in hepatocellular carcinoma. World J Surg Oncol. 2020;18(1):127.
Weber JS, et al. Toxicities of Immunotherapy for the Practitioner. J Clin Oncol. 2015;33(18):2092–9.
Oh JM, et al. U1 snRNP regulates cancer cell migration and invasion in vitro. Nat Commun. 2020;11(1):1.
Mabonga L, Kappo AP. The oncogenic potential of small nuclear ribonucleoprotein polypeptide G: a comprehensive and perspective view. Am J Transl Res. 2019;11(11):6702–16. Published 2019 Nov 15.
Tu J, et al. Bioinformatics analysis of molecular genetic targets and key pathways for hepatocellular carcinoma. Onco Targets Ther. 2019;12:5153–62.
Li L, et al. PACT is a negative regulator of p53 and essential for cell growth and embryonic development. Proc Natl Acad Sci U S A. 2007;104(19):7951–6.
Hull R, et al. The Drosophila retinoblastoma binding protein 6 family member has two isoforms and is potentially involved in embryonic patterning. Int J Mol Sci. 2015;16(5):10242–66.
Ezkurdia I, et al. Multiple evidence strands suggest that there may be as few as 19,000 human protein-coding genes. Hum Mol Genet. 2014;23(22):5866–78.
Ji M, et al. Small nuclear ribonucleoprotein polypeptide N accelerates malignant progression and poor prognosis in colorectal cancer transcriptionally regulated by E2F8. Front Oncol. 2020;10:561287. Published 2020 Nov 2. https://doi.org/10.3389/fonc.2020.561287
Pinero F, Dirchwolf M, Pessoa MG. Biomarkers in hepatocellular carcinoma: diagnosis, prognosis and treatment response assessment. Cells. 2020;9(6):1370. Published 2020 Jun 1. https://doi.org/10.3390/cells9061370
Yi M, et al. Biomarkers for predicting efficacy of PD-1/PD-L1 inhibitors. Mol Cancer. 2018;17(1):129.
Kaderbhai C, Tharin Z, Ghiringhelli F. The role of molecular profiling to predict the response to immune checkpoint inhibitors in lung cancer. Cancers (Basel). 2019;11(2):201. Published 2019 Feb 10. https://doi.org/10.3390/cancers11020201
Keenan TE, Burke KP, Van Allen EM. Genomic correlates of response to immune checkpoint blockade. Nat Med. 2019;25(3):389–402.
Wang S, et al. Antigen presentation and tumor immunogenicity in cancer immunotherapy response prediction. Elife. 2019;8:e49020. Published 2019 Nov 26. https://doi.org/10.7554/eLife.49020
Yi M, et al. The construction, expression, and enhanced anti-tumor activity of YM101: a bispecific antibody simultaneously targeting TGF-beta and PD-L1. J Hematol Oncol. 2021;14(1):27.
Wu K, et al. The efficacy and safety of combination of PD-1 and CTLA-4 inhibitors: a meta-analysis. Exp Hematol Oncol. 2019;8:26.
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Cai, J., Zhou, M. & Xu, J. N6-methyladenosine (m6A) RNA methylation regulator SNRPC is a prognostic biomarker and is correlated with immunotherapy in hepatocellular carcinoma. World J Surg Onc 19, 241 (2021). https://doi.org/10.1186/s12957-021-02354-8
- Hepatocellular carcinoma
- Immune checkpoint