Gene expression profiling reveals U1 snRNA regulates cancer gene expression

Small non-coding RNAs and pancreatic ductal adenocarcinoma: Linking diagnosis, pathogenesis, drug resistance, and therapeutic potential

This review comprehensively investigates the intricate interplay between small non-coding RNAs (sncRNAs) and pancreatic ductal adenocarcinoma (PDAC), a devastating malignancy with limited therapeutic options. Our analysis reveals the pivotal roles of sncRNAs in various facets of PDAC biology, spanning diagnosis, pathogenesis, drug resistance, and therapeutic strategies. sncRNAs have emerged as promising biomarkers for PDAC, demonstrating distinct expression profiles in diseased tissues. sncRNA differential expression patterns, often detectable in bodily fluids, hold potential for early and minimally invasive diagnostic approaches. Furthermore, sncRNAs exhibit intricate involvement in PDAC pathogenesis, regulating critical cellular processes such as proliferation, apoptosis, and metastasis. Additionally, mechanistic insights into sncRNA-mediated pathogenic pathways illuminate novel therapeutic targets and interventions. A significant focus of this review is dedicated to unraveling sncRNA mechanisms underlying drug resistance in PDAC. Understanding these mechanisms at the molecular level is imperative for devising strategies to overcome drug resistance. Exploring the therapeutic landscape, we discuss the potential of sncRNAs as therapeutic agents themselves as their ability to modulate gene expression with high specificity renders them attractive candidates for targeted therapy. In summary, this review integrates current knowledge on sncRNAs in PDAC, offering a holistic perspective on their diagnostic, pathogenic, and therapeutic relevance. By elucidating the roles of sncRNAs in PDAC biology, this review provides valuable insights for the development of novel diagnostic tools and targeted therapeutic approaches, crucial for improving the prognosis of PDAC patients.

Expression signature of ten small nuclear RNAs serves as novel biomarker for prognosis prediction of acute myeloid leukemia

This study aimed to screen for small nuclear RNAs (snRNAs) associated with the prognosis of acute myeloid leukemia (AML) by using The Cancer Genome Atlas (TCGA) whole-transcriptome sequencing dataset. A total of 130 AML patients from TCGA cohort with complete prognostic information and transcriptome data were enrolled in the current study. Comprehensive survival and functional enrichment analyses were performed to explore the prognostic value and potential biological functions of prognostic snRNAs in AML patients. In the current study, we screened 72 snRNAs that were notably associated with the clinical outcome of AML and developed an expression signature consist of ten snRNAs, that can be accurately applied to assess the overall survival of AML patients. Functional mechanism analysis revealed that this expression signature may be strongly linked to some classical tumor-associated pathways, such as Notch and Wnt pathways, as well as being closely related to B and T cell receptor pathways. Furthermore, we screened six compounds (chicago sky blue 6 B, 5230742, clorsulon, nefopam, nicardipine, and streptomycin) that may serve as targeted therapeutic drugs for AML using connectivity maps. Tumor immunoassays indicated significant differences in the immune microenvironment of the bone marrow tissue between high-risk and low-risk AML patients. Immune infiltration analysis also revealed significant differences in the abundance of multiple immune cells in the bone marrow of the two groups of AML patients groups. In conclusion, our results revealed a novel prognostic expression signature of AML consisting of ten snRNAs, and we conducted a preliminary exploration of its potential biological functions and tumor immunity.

RNA-Based Liquid Biopsy in Head and Neck Cancer

Head and neck cancer (HNC) is a prevalent and diverse group of malignancies with substantial morbidity and mortality rates. Early detection and monitoring of HNC are crucial for improving patient outcomes. Liquid biopsy, a non-invasive diagnostic approach, has emerged as a promising tool for cancer detection and monitoring. In this article, we review the application of RNA-based liquid biopsy in HNC. Various types of RNA, including messenger RNA (mRNA), microRNA (miRNA), long non-coding RNA (lncRNA), small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), circular RNA (circRNA) and PIWI-interacting RNA (piRNA), are explored as potential biomarkers in HNC liquid-based diagnostics. The roles of RNAs in HNC diagnosis, metastasis, tumor resistance to radio and chemotherapy, and overall prognosis are discussed. RNA-based liquid biopsy holds great promise for the early detection, prognosis, and personalized treatment of HNC. Further research and validation are necessary to translate these findings into clinical practice and improve patient outcomes.

RNA Sequencing of Intraoperative Peritumoral Tissues Reveals Potential Pathways Involved in Glioma-Related Seizures

Tumor-induced changes in the peritumoral neocortex play a crucial role in generation of seizures. This study aimed to investigate the molecular mechanisms potentially involved in peritumoral epilepsy in low-grade gliomas (LGGs). Intraoperative peritumoral brain tissues resected from LGG patients with seizures (pGRS) or without seizures (pGNS) were used for RNA sequencing (RNA-seq). Comparative transcriptomics was performed to identify differentially expressed genes (DEGs) in pGRS compared to pGNS using deseq2 and edgeR packages (R). Gene set enrichment analysis (GSEA) using Gene Ontology terms and Kyoto Encyclopedia of Genes & Genomes (KEGG) pathways was performed using the clusterProfiler package (R). The expression of key genes was validated at the transcript and protein levels in the peritumoral region using real-time PCR and immunohistochemistry, respectively. A total of 1073 DEGs were identified in pGRS compared to pGNS, of which 559 genes were upregulated and 514 genes were downregulated (log2 fold-change ≥ 2, padj < 0.001). The DEGs in pGRS were highly enriched in the "Glutamatergic Synapse" and "Spliceosome" pathways, with increased expression of GRIN2A (NR2A), GRIN2B (NR2B), GRIA1 (GLUR1), GRIA3 (GLUR3), GRM5, CACNA1C, CACNA1A, and ITPR2. Moreover, increased immunoreactivity was observed for NR2A, NR2B, and GLUR1 proteins in the peritumoral tissues of GRS. These findings suggest that altered glutamatergic signaling and perturbed Ca2+ homeostasis may be potential causes of peritumoral epilepsy in gliomas. This explorative study identifies important genes/pathways that merit further characterization for their potential involvement in glioma-related seizures.

Tumor educated platelet: the novel BioSource for cancer detection

Platelets, involved in the whole process of tumorigenesis and development, constantly absorb and enrich tumor-specific substances in the circulation during their life span, thus called "Tumor Educated Platelets" (TEPs). The alterations of platelet mRNA profiles have been identified as tumor markers due to the regulatory mechanism of post-transcriptional splicing. Small nuclear RNAs (SnRNAs), the important spliceosome components in platelets, dominate platelet RNA splicing and regulate the splicing intensity of pre-mRNA. Endogenous variation at the snRNA levels leads to widespread differences in alternative splicing, thereby driving the development and progression of neoplastic diseases. This review systematically expounds the bidirectional tumor-platelets interactions, especially the tumor induced alternative splicing in TEP, and further explores whether molecules related to alternative splicing such as snRNAs can serve as novel biomarkers for cancer diagnostics.

Benchmark study for evaluating the quality of reference genomes and gene annotations in 114 species

Introduction

For reference genomes and gene annotations are key materials that can determine the limits of the molecular biology research of a species; however, systematic research on their quality assessment remains insufficient.

Methods

We collected reference assemblies, gene annotations, and 3,420 RNA-sequencing (RNA-seq) data from 114 species and selected effective indicators to simultaneously evaluate the reference genome quality of various species, including statistics that can be obtained empirically during the mapping process of short reads. Furthermore, we newly presented and applied transcript diversity and quantification success rates that can relatively evaluate the quality of gene annotations of various species. Finally, we proposed a next-generation sequencing (NGS) applicability index by integrating a total of 10 effective indicators that can evaluate the genome and gene annotation of a specific species.

Results and discussion

Based on these effective evaluation indicators, we successfully evaluated and demonstrated the relative accessibility of NGS applications in all species, which will directly contribute to determining the technological boundaries in each species. Simultaneously, we expect that it will be a key indicator to examine the direction of future development through relative quality evaluation of genomes and gene annotations in each species, including countless organisms whose genomes and gene annotations will be constructed in the future.

Comprehensive analysis of molecular mechanism and a novel prognostic signature based on small nuclear RNA biomarkers in gastric cancer patients

Small nuclear RNAs (snRNAs) are rarely reported in cancer. This study is based on The Cancer Genome Atlas genome-wide data set to explore the prognostic value and molecular mechanism of snRNAs in gastric cancer (GC). Gene ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analysis were used to explore the molecular mechanism of snRNAs. A total of 351 patients were included in the survival analysis, and 14 prognostic snRNAs were identified using multivariate survival analysis. We constructed a prognostic signature containing nine snRNAs, which can signally classify patients into high- and low-risk phenotypes (adjusted P < 0.0001, hazard ratio = 2.671, 95% confidence interval = 1.850–3.858). Combining the molecular mechanisms obtained by the three functional enrichment approaches, we concluded that this prognostic signature snRNAs participated in classical tumor-related signaling pathways, including Notch, PI3K, toll-like receptor, etc.; cell adhesion; cell cycle; cell proliferation; and other biological processes that affect the biological phenotype of cancer cells. We also found significant downregulation of the abundance of immune cell infiltrates and immune microenvironment scores for high-risk phenotypes of GC patients. In conclusion, this study has identified 14 prognostic snRNAs signally associated with GC overall survival and also constructed a novel prognostic signature containing nine prognostic snRNAs.

Identification of heterogeneous nuclear ribonucleoprotein as a candidate biomarker for diagnosis and prognosis of hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common type of liver cancer with a high mortality rate. However, spliceosomal genes are still lacking in the diagnosis and prognosis of HCC.

Methods

Identification of differentially expressed genes (DEGs) was performed using the limma package in R software. Modules highly related to HCC were obtained by weighted gene co-expression network analysis (WGCNA), and the module genes were analyzed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The biomarker for diagnosing HCC was determined by receiver operating characteristic (ROC) curve analysis, and the effect of the biomarker in the diagnosis of HCC was evaluated by performing five-fold cross-validation with logistic regression. HCC specimens from preoperatively treated patients were tested for biomarker by real-time quantitative polymerase chain reaction (RT-qPCR). Kaplan-Meier analysis was used to assess the relationship between biomarker and patient survival. The role of biomarker was evaluated using ESTIMATE analysis in the tumor microenvironment.

Results

In this study, 389 DEGs were screened out from three Gene Expression Omnibus (GEO) datasets. We also found that the turquoise module of 123 genes from The Cancer Genome Atlas (TCGA) data was the key module with the highest correlation with HCC traits. Then, 123 genes were analyzed using the KEGG enrichment pathway, and eight genes were found to be most significantly related to the spliceosome pathway. We selected 8 genes and 389 DEGs shared genes, and finally got the only gene, heterogeneous nuclear ribonucleoprotein (hnRNPU). The high expression of hnRNPU was associated with poor prognosis of HCC, and hnRNPU was a biomarker for diagnosing HCC. In the tissues of patients with excellent HCC treatment hnRNPU messenger RNA (mRNA) was lower than in the tissues of patients with poor HCC treatment. High expression of hnRNPU was significantly increased in HCC patients with low stromal (P<0.05), low immune (P<0.05), and low estimation scores (P<0.05), and with high tumor purity (P<0.05) and high malignant progression (P<0.05) of the HCC.

Conclusions

The hnRNPU gene identified in this study may become a new biomarker for the diagnosis and prognosis of HCC.

Breast Cancer and the Other Non-Coding RNAs

Breast cancer is very heterogenous and the most common gynaecological cancer, with various factors affecting its development. While its impact on human lives and national health budgets is still rising in almost all global areas, many molecular mechanisms affecting its onset and development remain unclear. Conventional treatments still prove inadequate in some aspects, and appropriate molecular therapeutic targets are required for improved outcomes. Recent scientific interest has therefore focused on the non-coding RNAs roles in tumour development and their potential as therapeutic targets. These RNAs comprise the majority of the human transcript and their broad action mechanisms range from gene silencing to chromatin remodelling. Many non-coding RNAs also have altered expression in breast cancer cell lines and tissues, and this is often connected with increased proliferation, a degraded extracellular environment, and higher endothelial to mesenchymal transition. Herein, we summarise the known abnormalities in the function and expression of long non-coding RNAs, Piwi interacting RNAs, small nucleolar RNAs and small nuclear RNAs in breast cancer, and how these abnormalities affect the development of this deadly disease. Finally, the use of RNA interference to suppress breast cancer growth is summarised.

Small non-coding RNAs in human cancer: function, clinical utility, and characterization

Small non-coding RNAs (sncRNAs) play critical roles in multiple regulatory processes, including transcription, post-transcription, and translation. Emerging evidence reveals the critical roles of sncRNAs in cancer development and their potential role as biomarkers and/or therapeutic targets. In this paper, we review recent research on four sncRNA species with functional significance in cancer: small nucleolar RNAs, transfer RNA, small nuclear RNAs, and piwi-interacting RNAs. We introduce their functional roles in tumorigenesis and discuss the potential utility of sncRNAs as prognostic and diagnostic biomarkers and therapeutic targets. We further summarize approaches to characterize sncRNAs in a high-throughput manner, including the specific library construction and computational framework. Our review provides a perspective of the functions, clinical utility, and characterization of sncRNAs in cancer.

The Role of RNA Splicing Factors in Cancer: Regulation of Viral and Human Gene Expression in Human Papillomavirus-Related Cervical Cancer

The spliceosomal complex components, together with the heterogeneous nuclear ribonucleoproteins (hnRNPs) and serine/arginine-rich (SR) proteins, regulate the process of constitutive and alternative splicing, the latter leading to the production of mRNA isoforms coding multiple proteins from a single pre-mRNA molecule. The expression of splicing factors is frequently deregulated in different cancer types causing the generation of oncogenic proteins involved in cancer hallmarks. Cervical cancer is caused by persistent infection with oncogenic human papillomaviruses (HPVs) and constitutive expression of viral oncogenes. The aberrant activity of hnRNPs and SR proteins in cervical neoplasia has been shown to trigger the production of oncoproteins through the processing of pre-mRNA transcripts either derived from human genes or HPV genomes. Indeed, hnRNP and SR splicing factors have been shown to regulate the production of viral oncoprotein isoforms necessary for the completion of viral life cycle and for cell transformation. Target-therapy strategies against hnRNPs and SR proteins, causing simultaneous reduction of oncogenic factors and inhibition of HPV replication, are under development. In this review, we describe the current knowledge of the functional link between RNA splicing factors and deregulated cellular as well as viral RNA maturation in cervical cancer and the opportunity of new therapeutic strategies.

U1 snRNP regulates cancer cell migration and invasion in vitro

Stimulated cells and cancer cells have widespread shortening of mRNA 3'-untranslated regions (3'UTRs) and switches to shorter mRNA isoforms due to usage of more proximal polyadenylation signals (PASs) in introns and last exons. U1 snRNP (U1), vertebrates' most abundant non-coding (spliceosomal) small nuclear RNA, silences proximal PASs and its inhibition with antisense morpholino oligonucleotides (U1 AMO) triggers widespread premature transcription termination and mRNA shortening. Here we show that low U1 AMO doses increase cancer cells' migration and invasion in vitro by up to 500%, whereas U1 over-expression has the opposite effect. In addition to 3'UTR length, numerous transcriptome changes that could contribute to this phenotype are observed, including alternative splicing, and mRNA expression levels of proto-oncogenes and tumor suppressors. These findings reveal an unexpected role for U1 homeostasis (available U1 relative to transcription) in oncogenic and activated cell states, and suggest U1 as a potential target for their modulation.

Aberrantly expressed snoRNA, snRNA, piRNA and tRFs in canine melanoma

Among small non-coding RNAs (sncRNAs/sRNAs), the functional regulation of microRNAs (miRNAs) has been studied in canine oral melanoma (COM). However, the expression level of other sncRNAs, like small nucleolar RNAs (snoRNAs), small nuclear RNAs (snRNAs), transfer RNA-derived fragments (tRFs) and PIWI-interacting RNAs (piRNAs), in COM is unknown. The aim of this study was to investigate sncRNAs other than miRNAs in COM from our small RNA sequencing project (PRJNA516252). We found that several snRNAs and piRNAs were upregulated, whereas tRFs and snoRNAs were downregulated in COM. Upregulation of U1 snRNA and piR-972, and downregulation of tRNA-ser (1) and snoRA24 was confirmed in dog melanoma tissue and cell lines by quantitative reverse transcription PCR. Consistently, the expression of tRNA-ser (1) and snoRA24 in plasma of COM cases was also decreased. Finally, we found a similar expression trend of U1 and snoRA24 in the human cutaneous melanoma cell line, MEWO, compared with human epidermal melanocyte cells (HEMa-Lp). In our study, snRNA, snoRNA, tRFs and piRNA were dysregulated during melanoma progression. Moreover, the melanoma-associated expression of U1 and snoRA24 was similar in human and dog melanoma.

RNA components of the spliceosome regulate tissue- and cancer-specific alternative splicing

Alternative splicing of pre-mRNAs plays a pivotal role during the establishment and maintenance of human cell types. Characterizing the trans-acting regulatory proteins that control alternative splicing has therefore been the focus of much research. Recent work has established that even core protein components of the spliceosome, which are required for splicing to proceed, can nonetheless contribute to splicing regulation by modulating splice site choice. We here show that the RNA components of the spliceosome likewise influence alternative splicing decisions. Although these small nuclear RNAs (snRNAs), termed U1, U2, U4, U5, and U6 snRNA, are present in equal stoichiometry within the spliceosome, we found that their relative levels vary by an order of magnitude during development, across tissues, and across cancer samples. Physiologically relevant perturbation of individual snRNAs drove widespread gene-specific differences in alternative splicing but not transcriptome-wide splicing failure. Genes that were particularly sensitive to variations in snRNA abundance in a breast cancer cell line model were likewise preferentially misspliced within a clinically diverse cohort of invasive breast ductal carcinomas. As aberrant mRNA splicing is prevalent in many cancers, we propose that a full understanding of such dysregulated pre-mRNA processing requires study of snRNAs, as well as protein splicing factors. Together, our data show that the RNA components of the spliceosome are not merely basal factors, as has long been assumed. Instead, these noncoding RNAs constitute a previously uncharacterized layer of regulation of alternative splicing, and contribute to the establishment of global splicing programs in both healthy and malignant cells.

Metformin inhibits Aβ25-35 -induced apoptotic cell death in SH-SY5Y cells

Metformin, a first-line drug for type-2 diabetes, plays a potentially protective role in preventing Alzheimer's disease (AD), but its underlying mechanism is unclear. In this study, Aβ_25-35 -treated SH-SY5Y cells were used as a cell model of AD to investigate the neuroprotective effect of metformin, as well as its underlying mechanisms. We found that metformin decreased the cell apoptosis rate and death, ratio of Bcl-2/Bax, and expression of NR2A and NR2B, and increased the expression of LC3 in Aβ_25-35 -treated SH-SY5Y cells. Metformin also reduced intracellular and extracellular Glu concentrations, as well as the intracellular concentration of Ca²⁺ and ROS in Aβ_25-35 -treated SH-SY5Y cells. These findings suggest that metformin inhibits Aβ_25-35 -treated SH-SY5Y cell death by inhibiting apoptosis, decreasing intracellular Ca²⁺ and ROS by reducing neurotoxicity of excitatory amino acids, and by possibly reversing autophagy disorder via regulating autophagy process.

Transcriptomic Sequencing of Airway Epithelial Cell NCI-H292 Induced by Synthetic Cationic Polypeptides

Eosinophil asthma is characterized by the infiltration of eosinophils to the bronchial epithelium. The toxic cationic protein released by eosinophils, mainly major basic protein (MBP), is one of the most important causative factors of epithelium damage. Poly-L-Arginine (PLA) is a kind of synthetic cationic polypeptides, which is widely used to mimic the effects of MBP on epithelial cells in vitro. However, little is known about the changes of differentially expressed genes (DEGs) and transcriptome profiles in cationic protein stimulated epithelial cells. In this study, we compared the expression of DEGs and transcriptome profiles between PLA-treated airway epithelial cells NCI-H292 and control. The results showed that there were a total of 230 DEGs, of which 86 were upregulated and 144 were downregulated. These DEGs were further analyzed using gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The results showed that the upregulated DEGs were involved in cholesterol synthesis, protein binding, and composition of cellular membranes, mainly enriched in metabolic and biosynthesis pathways. While downregulated DEGs were implicated in cell adhesion, extracellular matrix (ECM) composition and cytoskeleton and were enriched in ECM pathway. In conclusion, our research provided the mechanism of the cationic polypeptides acting on the airway epithelial cells on the basis of transcriptomic profile, and this could be regarded as important indications in unveiling the pathologic role of natural cationic proteins in the damage to epithelial cells of asthmatics.

U1 snRNA over-expression affects neural oscillations and short-term memory deficits in mice

Small nuclear RNAs (snRNAs) and other RNA spliceosomal components are involved in neurological and psychiatric disorders. U1 snRNA has recently been demonstrated to be altered in pathology in some neurodegenerative diseases, but whether it has a causative role is not clear. Here we have studied this by overexpressing U1 snRNA in mice and measured their hippocampal oscillatory patterns and brain functions. Novel object recognition test showed that the recognition index was significantly decreased in the U1 snRNA over-expression mice compared to that in the C57BL mice. U1 snRNA over-expression regulated not only the pattern of neural oscillations but also the expression of neuron excitatory and inhibitory proteins. Here we show that U1 snRNA over-expression contains the shrinkage distribution of theta-power, theta-phase lock synchronization, and theta and low-gamma cross-frequency coupling in the hippocampus. The alternations of neuron receptors by the U1 snRNA overexpression also modulated the decreasing of recognition index, the energy distribution of theta power spectrum with the reductions of theta phase synchronization and phase-amplitude coupling between theta and low-gamma. Linking these all together, our results suggest that U1 snRNA overexpression particularly causes a deficit in short-term memory. These findings make a bedrock of our research that U1 snRNA bridges the gap about the mechanism behind short-term memory based on the molecular and mesoscopic level.

Using Pan RNA-Seq Analysis to Reveal the Ubiquitous Existence of 5' and 3' End Small RNAs

In this study, we used pan RNA-seq analysis to reveal the ubiquitous existence of both 5′ and 3′ end small RNAs (5′ and 3′ sRNAs). 5′ and 3′ sRNAs alone can be used to annotate nuclear non-coding and mitochondrial genes at 1-bp resolution and identify new steady RNAs, which are usually transcribed from functional genes. Then, we provided a simple and cost effective way for the annotation of nuclear non-coding and mitochondrial genes and the identification of new steady RNAs, particularly long non-coding RNAs (lncRNAs). Using 5′ and 3′ sRNAs, the annotation of human mitochondrial was corrected and a novel ncRNA named non-coding mitochondrial RNA 1 (ncMT1) was reported for the first time in this study. We also found that most of human tRNA genes have downstream lncRNA genes as lncTRS-TGA1-1 and corrected the misunderstanding of them in previous studies. Using 5′, 3′, and intronic sRNAs, we reported for the first time that enzymatic double-stranded RNA (dsRNA) cleavage and RNA interference (RNAi) might be involved in the RNA degradation and gene expression regulation of U1 snRNA in human. We provided a different perspective on the regulation of gene expression in U1 snRNA. We also provided a novel view on cancer and virus-induced diseases, leading to find diagnostics or therapy targets from the ribonuclease III (RNase III) family and its related pathways. Our findings pave the way toward a rediscovery of dsRNA cleavage and RNAi, challenging classical theories.

Prognostic value of small nuclear RNAs (snRNAs) for digestive tract pan- adenocarcinomas identified by RNA sequencing data

Malignant tumors of the digestive tract include esophageal, gastric, and colorectal carcinomas, which all have high global mortality rates. A clinical role for small nuclear RNA (snRNA), a type of small non-coding RNA, has not yet been documented for digestive tract pan-adenocarcinomas. Therefore, the aim of the study was to identify differentially expressed snRNAs and to explore their prognostic implications in pan-adenocarcinomas from the esophagus, stomach, colon, and rectum. The pan-carcinoma RNA-sequencing data of four types of digestive tract cancers with 1, 102 cases obtained from The Cancer Genome Atlas (TCGA) project were analyzed and the differentially expressed snRNAs were evaluated using the edgeR package. The prognostic value of each of the selected snRNAs was determined by univariate and multivariate Cox regression analyses. All the digestive tract pan-adenocarcinomas showed differential expression of three snRNAs: the up-regulated RNU1-106 P and RNU6-850 P and the down-regulated RNU6-529 P. Interestingly, RNU6-101 P appeared to be a risk factor for esophageal adenocarcinoma (ESAD) and RNVU1-4 was potentially a protective factor for stomach adenocarcinoma (STAD) survival. This consistent finding of differential expression of all three snRNAs in all four types of digestive system cancers suggests potential roles for these snRNAs in the tumorigenesis of digestive system cancers. RNU6-101 P could play a pivotal role in the progression of ESAD and RNVU1-4 could perform a protective role in STAD. However, since the current findings were based on RNA-sequencing data mining, more studies are needed for verification.