Neuronal vulnerability and multilineage diversity in multiple sclerosis.

Multiple sclerosis (MS) is a neuroinflammatory disease with a relapsing-remitting disease course at early stages, distinct lesion characteristics in cortical grey versus subcortical white matter and neurodegeneration at chronic stages. Here we used single-nucleus RNA sequencing to assess changes in expression in multiple cell lineages in MS lesions and validated the results using multiplex in situ hybridization. We found selective vulnerability and loss of excitatory CUX2-expressing projection neurons in upper-cortical layers underlying meningeal inflammation; such MS neuron populations exhibited upregulation of stress pathway genes and long non-coding RNAs. Signatures of stressed oligodendrocytes, reactive astrocytes and activated microglia mapped most strongly to the rim of MS plaques. Notably, single-nucleus RNA sequencing identified phagocytosing microglia and/or macrophages by their ingestion and perinuclear import of myelin transcripts, confirmed by functional mouse and human culture assays. Our findings indicate lineage- and region-specific transcriptomic changes associated with selective cortical neuron damage and glial activation contributing to progression of MS lesions.

Direct Determination of Pseudouridine in RNA by Mass Spectrometry Coupled with Stable Isotope Labeling.

Pseudouridine (Ψ) is the only "mass-silent" nucleoside produced by post-transcriptional RNA modification. We developed a mass spectrometry (MS)-based technique coupled with in vivo deuterium (D) labeling of uridines for direct determination of Ψs in cellular RNA and applied it to the comprehensive analysis of post-transcriptional modifications in human ribosomal RNAs. The method utilizes human TK6/mouse FM3A cells deficient in uridine monophosphate synthase using a CRISPR-Cas9 technique to turn off de novo uridine synthesis and fully labels uridines with D at uracil positions 5 and 6 by cultivating the cells in a medium containing uridine-5,6-D2. The pseudouridylation reaction in those cells results in the exchange of the D at the C5 of uracil with hydrogen from solvent, which produces a -1 Da mass shift, thus allowing MS-based determination of RNA Ψs. We present here the experimental details of this method and show that it allows the identification of all Ψs in human major nuclear and nucleolar RNAs, including several previously unknown Ψs. Because the method allows direct determination of Ψs at the femtomole level of RNA, it will serve as a useful tool for structure/function studies of a wide variety of noncoding RNAs.

Differential long-term stability of microRNAs and RNU6B snRNA in 12-20 year old archived formalin-fixed paraffin-embedded specimens.

BACKGROUND: The quantitative analysis of microRNA (miRNA) gene expression in archived formalin-fixed, paraffin embedded (FFPE) tissues has been instrumental to identifying their potential roles in cancer biology, diagnosis, and prognosis. However, it remains unclear whether miRNAs remain stable in FFPE tissues stored for long periods of time. METHODS: Here we report Taqman real-time RT-PCR quantification of miR-21, miR-141, miR-221, and RNU6B small nuclear RNA (snRNA) levels from 92 radical prostatectomy specimens stored for 12-20 years in FFPE blocks. The relative stability of each transcript over time was assessed using general linear models. The correlation between transcript quantities, sample age, and RNA integrity number (RIN) were determined utilizing Spearman rank correlation. RESULTS: All transcript levels linearly decreased with sample age, demonstrating a clear loss of miRNA stability and RNU6B snRNA stability over time. The most rapid rates of degradation were observed for RNU6B and miR-21, while miR-141 and miR-221 were more stable. RNA quality was not correlated with sample age or with miR-21, miR-221, or RNU6B snRNA levels. Conversely, miR-141 levels increased with RNA quality. CONCLUSIONS: MiRNA and snRNA levels gradually decreased over an eight year period in FFPE tissue blocks. Sample age was the most consistent feature associated with miRNA stability. The reference snRNA, RUN6B, was more rapidly degraded when compared to miR-141 and miR-221 miRNAs. Various miRNAs demonstrated differential rates of degradation. Quantitative miRNA studies from long-term archived FFPE tissues may therefore benefit from epidemiologic study design or statistical analysis methods that take into account differential storage-dependent transcript degradation.

[Correlation between RNA Expression Level and Early PMI in Human Brain Tissue].

OBJECTIVES: To explore the correlation between the expression levels of several RNA markers in human brain tissue and early postmortem interval （PMI）. METHODS: Twelve individuals with known PMI （range from 4.3 to 22.5 h） were selected and total RNA was extracted from brain tissue. Eight commonly used RNA markers were chosen including ß-actin, GAPDH, RPS29, 18S rRNA, 5S rRNA, U6 snRNA, miRNA-9 and miRNA-125b, and the expression levels were detected in brain tissue by real-time fluorescent quantitative PCR. The internal reference markers with stable expression in early PMI were screened using geNorm software and the relationship between its expression level and some relevant factors such as age, gender and cause of death were analyzed. RNA markers normalized by internal reference were inserted into the mathematic model established by previous research for PMI estimation using R software. Model quality was judged by the error rate calculated with estimated PMI. RESULTS: 5S rRNA, miRNA-9 and miRNA-125b showed quite stable expression and their expression levels had no relation with age, gender and cause of death. The error rate of estimated PMI using ß-actin was 24.6%, while GAPDH was 41.0%. CONCLUSIONS: 5S rRNA, miRNA-9 and miRNA-125b are suitable as internal reference markers of human brain tissue owing to their stable expression in early PMI. The expression level of ß-actin correlates well with PMI, which can be used as an additional index for early PMI estimation.

Isolation of an active step I spliceosome and composition of its RNP core.

Formation of catalytically active RNA structures within the spliceosome requires the assistance of proteins. However, little is known about the number and nature of proteins needed to establish and maintain the spliceosome's active site. Here we affinity-purified human spliceosomal C complexes and show that they catalyse exon ligation in the absence of added factors. Comparisons of the composition of the precatalytic versus the catalytic spliceosome revealed a marked exchange of proteins during the transition from the B to the C complex, with apparent stabilization of Prp19-CDC5 complex proteins and destabilization of SF3a/b proteins. Disruption of purified C complexes led to the isolation of a salt-stable ribonucleoprotein (RNP) core that contained both splicing intermediates and U2, U5 and U6 small nuclear RNA plus predominantly U5 and human Prp19-CDC5 proteins and Prp19-related factors. Our data provide insights into the spliceosome's catalytic RNP domain and indicate a central role for the aforementioned proteins in sustaining its catalytically active structure.

A homozygous mutation in RNU4ATAC as a cause of microcephalic osteodysplastic primordial dwarfism type I (MOPD I) with associated pigmentary disorder.

The designation microcephalic osteodysplastic primordial dwarfism (MOPD) refers to a group of autosomal recessive disorders, comprising microcephaly, growth retardation, and a skeletal dysplasia. The different types of MOPD have been delineated on the basis of clinical, radiological, and genetic criteria. We describe two brothers, born to healthy, consanguineous parents, with intrauterine and postnatal growth retardation, microcephaly with abnormal gyral pattern and partial agenesis of corpus callosum, and skeletal anomalies reminiscent of those described in MOPD type I. This was confirmed by the identification of the homozygous g.55G > A mutation of RNU4ATAC encoding U4atac snRNA. The sibs had yellowish-gray hair, fair skin, and deficient retinal pigmentation. Skin biopsy showed abnormal melanin function but OCA genes were normal. The older sib had an intracranial hemorrhage at 1 week after birth, the younger developed chilblains-like lesions at the age 2½ years old but analysis of the SAMHD1 and TREX1 genes did not show any mutations. To the best of our knowledge, vasculopathy and pigmentary disorders have not been reported in MOPD I.

snRNA-specific role of SMN in trypanosome snRNP biogenesis in vivo.

Pre-mRNA splicing in trypanosomes requires the SMN-mediated assembly of small nuclear ribonucleoproteins (snRNPs). In contrast to higher eukaryotes, the cellular localization of snRNP biogenesis and the involvement of nuclear-cytoplasmic trafficking in trypanosomes are controversial. By using RNAi knockdown of SMN in T. brucei to investigate its functional role in snRNP assembly, we found dramatic changes in the steady-state levels of snRNAs and snRNPs: The SL RNA accumulates, whereas U1, U4, and U5 snRNA levels decrease, and Sm core assembly in particular of the SL RNA is strongly reduced. In addition, SMN depletion blocks U4/U6 di-snRNP formation; the variant Sm core of the U2 snRNP, however, still forms efficiently after SMN knockdown. Concerning the longstanding question, whether nuclear-cytoplasmic trafficking is involved in trypanosomal snRNP biogenesis, fluorescence in situ hybridization (FISH) and immunofluorescence assays revealed that the SL RNA genes and transcripts colocalize with SMN. Remarkably, SMN silencing leads to a nucleoplasmic accumulation of both SL RNA and the Sm proteins. In sum, our data demonstrate an essential and snRNA-selective role of SMN in snRNP biogenesis in vivo and strongly argue for a nucleoplasmic Sm core assembly of the SL RNP.

Minor spliceosome components are predominantly localized in the nucleus.

Recently, it has been reported that there is a differential subcellular distribution of components of the minor U12-dependent and major U2-dependent spliceosome, and further that the minor spliceosome functions in the cytoplasm. To study the subcellular localization of the snRNA components of both the major and minor spliceosomes, we performed in situ hybridizations with mouse tissues and human cells. In both cases, all spliceosomal snRNAs were nearly exclusively detected in the nucleus, and the minor U11 and U12 snRNAs were further shown to colocalize with U4 and U2, respectively, in human cells. Additionally, we examined the distribution of several spliceosomal snRNAs and proteins in nuclear and cytoplasmic fractions isolated from human cells. These studies revealed an identical subcellular distribution of components of both the U12- and U2-dependent spliceosomes. Thus, our data, combined with several earlier publications, establish that, like the major spliceosome, components of the U12-dependent spliceosome are localized predominantly in the nucleus.

Large-scale integration of single-cell transcriptomic data captures transitional progenitor states in mouse skeletal muscle regeneration.

Skeletal muscle repair is driven by the coordinated self-renewal and fusion of myogenic stem and progenitor cells. Single-cell gene expression analyses of myogenesis have been hampered by the poor sampling of rare and transient cell states that are critical for muscle repair, and do not inform the spatial context that is important for myogenic differentiation. Here, we demonstrate how large-scale integration of single-cell and spatial transcriptomic data can overcome these limitations. We created a single-cell transcriptomic dataset of mouse skeletal muscle by integration, consensus annotation, and analysis of 23 newly collected scRNAseq datasets and 88 publicly available single-cell (scRNAseq) and single-nucleus (snRNAseq) RNA-sequencing datasets. The resulting dataset includes more than 365,000 cells and spans a wide range of ages, injury, and repair conditions. Together, these data enabled identification of the predominant cell types in skeletal muscle, and resolved cell subtypes, including endothelial subtypes distinguished by vessel-type of origin, fibro-adipogenic progenitors defined by functional roles, and many distinct immune populations. The representation of different experimental conditions and the depth of transcriptome coverage enabled robust profiling of sparsely expressed genes. We built a densely sampled transcriptomic model of myogenesis, from stem cell quiescence to myofiber maturation, and identified rare, transitional states of progenitor commitment and fusion that are poorly represented in individual datasets. We performed spatial RNA sequencing of mouse muscle at three time points after injury and used the integrated dataset as a reference to achieve a high-resolution, local deconvolution of cell subtypes. We also used the integrated dataset to explore ligand-receptor co-expression patterns and identify dynamic cell-cell interactions in muscle injury response. We provide a public web tool to enable interactive exploration and visualization of the data. Our work supports the utility of large-scale integration of single-cell transcriptomic data as a tool for biological discovery.

In-gel digestion for mass spectrometric characterization of RNA from fluorescently stained polyacrylamide gels.

Although current mass spectrometry-based proteomics technology allows for high-throughput analysis of protein components in functional ribonucleoprotein complexes, this technology has had limited application to studies of RNA itself. Here we present a protocol for RNA analysis using polyacrylamide gel electrophoresis coupled with liquid chromatography-tandem mass spectrometry. Specifically, RNAs of interest are subjected to polyacrylamide gel electrophoresis and stained with a fluorescent dye, and RNAs in gel bands are digested with nuclease and then analyzed directly liquid chromatography-mass spectrometry, resulting in highly accurate mass values and reliable information on post-transcriptional modifications. We demonstrate that the method can be applied to the identification and chemical analysis of small RNAs in mouse embryonic stem cell extracts and of small RNAs in the spliceosomal ribonucleoprotein complex pulled down from yeast cells using a tagged protein cofactor as bait. The protocol is relatively simple and allowed us to identify not only three novel methylated nucleotide residues of RNase P RNA, U6 snRNA, and 7SL RNA prepared from mouse ES cells but also various 3'-end forms of U4, U5S, and U6 snRNAs isolated from the yeast spliceosome at the femtomole level. The method is thus a convenient tool for direct analysis of RNAs in various cellular ribonucleoprotein complexes, particularly for the analysis of post-transcriptional modifications and metabolic processing of RNA.

Identification of a new microRNA expression profile as a potential cancer screening tool.

PURPOSE: Small non-coding microRNAs (miRNAs) are key components of cancer development and are considered as potential biomarkers for cancer diagnosis and treatment monitoring. This study investigated miRNA expression profiles of human cancer cells in order to develop a screening method for lung cancer. METHODS: A series of lung cancer related miRNAs (miR-21, miR-145, miR-155, miR-205, miR-210, miR-92, miR-17-5p, miR-143, miR-182, miR-372, let-7a) were selected as candidates for miRNA expression profiles of human lung cancer cell lines (A549, SK-mes-1). MicroRNA u6 was the endogenous control. Cancer cell lines for positive controls; breast MCF-7, prostate Du-145, and glioblastoma U118. The negative control was normal lung fibroblast cell line MRC-5. RT-PCR was performed on StepOnePlus (Applied Biosystem, USA). MiRNA expressions of malignant cells were compared with normal fibroblast cells as well as endogenous control (u6) using the thermal cycle at threshold. Assessment of miRNA expression profiles were then performed using agglomerative hierarchical cluster analysis software (SPSS13, USA). RESULTS: We demonstrated that miR-21, miR-182 and let7-5a were over-expressed, and miR-145 and miR-155 were under-expressed in all cancer cell lines. Combined with the cluster analysis we were able to clearly distinguish cell lines for normal fibroblasts, breast cancer, prostate cancer, glioblastoma, and lung cancer. CONCLUSION: There is potential utility of screening for lung cancer with miRNA expression profiles. Future work will focus on the sensitivity of such miRNA expression profiles in screening sputum for lung cancer, which can be performed in real time.

Evaluating the effect of body fluid mixture on the relative expression ratio of blood-specific RNA markers.

The estimation of the time elapsed since a biological stain was deposited at a crime scene can provide crucial information to a forensic investigation, indicating either when a crime was committed, or whether the biological evidence was deposited at the time of a known crime event. This would enable the investigators to limit the number of suspects and to assess alibis. The relative expression ratios (RERs) of body fluid-specific RNA markers are promising molecular tools for indicating the age of biological stains. However, the nature of some forensic samples found at crime scenes could be challenging, as they frequently occur in a mixture of different body fluid types. The research presented here has utilised reverse transcription quantitative PCR (RT-qPCR) to explore the impact of bloodstains being present in mixtures with other body fluids (saliva or semen) on the resulting RERs of blood-specific markers. The expression level of three blood-specific markers (HBA, HBB and miR16) along with two reference genes (18S and U6) were analysed across multiple ageing time points in pure and mixed bloodstains. For some markers, no significant differences were found when comparing RERs in pure and mixed bloodstains, however some RERs were altered in mixed stains. This indicates that the presence of body fluid mixtures may have a significant effect on the RERs of some blood-specific markers. This should therefore be considered when selecting markers for estimating the age of stains, particularly when multiple body fluids are thought to be present.

Nucleoplasmic organization of small nuclear ribonucleoproteins in cultured human cells.

The organization of eight small nuclear ribonucleoproteins (the U1, U2, U4, U5, and U6 RNAs previously studied by others and three additional snRNAs, U11, U12, and 7SK) has been investigated in cultured human cells by fluorescence in situ hybridization with antisense DNA and 2'-O-Me RNA oligonucleotides. Using highly sensitive digital imaging microscopy we demonstrate that all of these snRNAs are widespread throughout the nucleoplasm, but they are excluded from the nucleoli. In addition, the U2, U4, U5, U6, and U12 snRNAs are concentrated in discrete nuclear foci, known as coiled bodies, but U1 and 7SK are not. In addition to coiled bodies, a classic speckled pattern was observed in the nucleoplasm of monolayer-grown HeLa cells, whereas suspension-grown HeLa cells revealed a more diffuse nucleoplasmic labeling. Immunofluorescence staining using various snRNP-specific antisera shows complete agreement with that of their antisense snRNA oligonucleotide counterparts. Although U2 RNA is concentrated in coiled bodies, quantitation of the fluorescence signals from the U2 antisense probe reveals that the bulk of the U2 snRNP is located in the nucleoplasm. Furthermore, simultaneous visualization of the U2 snRNAs and the tandemly repeated U2 genes demonstrates that coiled bodies are not the sites of U2 transcription.

Sea urchin maternal and embryonic U1 RNAs are spatially segregated in early embryos.

We have used in situ hybridization and cell fractionation methods to follow the distribution of U1 RNA and immunofluorescence microscopy to follow the distribution of snRNP proteins in oocytes, eggs, and embryos of several sea urchin species. U1 RNA and U1-specific snRNP antigens are concentrated in germinal vesicles of oocytes. Both appear to relocate after oocyte maturation because they are found primarily, if not exclusively, in the cytoplasm of mature unfertilized eggs. This cytoplasmic residence is maintained during early cleavage and U1 RNA is first detectable in nuclei of micromeres at the 16-cell stage. Between morula and gastrula stages the steady-state concentrations of both RNA and antigens gradually increase in nuclei and decrease in cytoplasm. Surprisingly, analysis of the distribution of newly synthesized U1 RNA shows that it does not equilibrate with the maternal pool. Instead new transcripts are confined to nuclei, while cytoplasmic U1 RNAs are of maternal origin. This lack of equilibration and the conversion of maternal U1 RNAs from nuclear species in oocytes to cytoplasmic in embryos suggests that these RNPs (or RNAs) are structurally altered when released to the cytoplasm at oocyte maturation.

Small nuclear ribonucleoproteins and heterogeneous nuclear ribonucleoproteins in the amphibian germinal vesicle: loops, spheres, and snurposomes.

We have examined the distribution of snRNPs in the germinal vesicle (GV) of frogs and salamanders by immunofluorescent staining and in situ nucleic acid hybridization. The major snRNAs involved in pre-mRNA splicing (U1, U2, U4, U5, and U6) occur together in nearly all loops of the lampbrush chromosomes, and in hundreds to thousands of small granules (1-4 microns diameter) suspended in the nucleoplasm. The loops and granules also contain several antigens that are regularly associated with snRNAs or spliceosomes (the Sm antigen, U1- and U2-specific antigens, and the splicing factor SC35). A second type of granule, often distinguishable by morphology, contains only U1 snRNA and associated antigens. We propose the term "snurposome" to describe the granules that contain snRNPs ("snurps"). Those that contain only U1 snRNA are A snurposomes, whereas those that contain all the splicing snRNAs are B snurposomes. GVs contain a third type of snRNP granule, which we call the C snurposome. C snurposomes range in size from less than 1 micron to giant structures greater than 20 microns in diameter. Usually, although not invariably, they have B snurposomes on their surface. They may also contain from one to hundreds of inclusions. Because of their remarkably spherical shape, C snurposomes with their associated B snurposomes have long been referred to as spheres or sphere organelles. Most spheres are free in the nucleoplasm, but a few are attached to chromosomes at specific chromosome loci, the sphere organizers (SOs). The relationship of sphere organelles to other snRNP-containing structures in the GV is obscure. We show by immunofluorescent staining that the lampbrush loops and B snurposomes also react with antibodies against heterogeneous nuclear ribonucleoproteins (hnRNPs). Transcription units on the loops are uniformly stained by anti-hnRNP and anti-snRNP antibodies, suggesting that nascent transcripts are associated with hnRNPs and snRNPs along their entire length, perhaps in the form of a unitary hnRNP/snRNP particle. That B snurposomes contain so many components involved in pre-mRNA packaging and processing suggests that they may serve as sites for assembly and storage of hnRNP/snRNP complexes destined for transport to the nascent transcripts on the lampbrush chromosome loops.

Saccharomyces cerevisiae has a U1-like small nuclear RNA with unexpected properties.

Previous experiments indicated that only a small subset of the approximately equal to 24 small nuclear RNAs (snRNAs) in Saccharomyces cerevisiae have binding sites for the Sm antigen, a hallmark of metazoan small nuclear ribonucleoproteins (snRNPs) involved in pre-messenger RNA splicing. Antibodies from human serum to Sm proteins were used to show that four snRNAs (snR7, snR14, snR19, and snR20) can be immunoprecipitated from yeast extracts. Three of these four, snR7, snR14, and snR20, have been shown to be analogs of mammalian U5, U4, and U2, respectively. Several regions of significant homology to U1 (164 nucleotides) have now been found in cloned and sequenced snR19 (568 nucleotides). These include ten out of ten matches to the 5' end of U1, the site known to interact with the 5' splice site of mammalian introns. Surprisingly, the precise conservation of this sequence precludes perfect complementarity between snR19 and the invariant yeast 5' junction (GTATGT), which differs from the mammalian consensus at the fourth position (GTPuAGT).

A computational screen for methylation guide snoRNAs in yeast.

Small nucleolar RNAs (snoRNAs) are required for ribose 2'-O-methylation of eukaryotic ribosomal RNA. Many of the genes for this snoRNA family have remained unidentified in Saccharomyces cerevisiae, despite the availability of a complete genome sequence. Probabilistic modeling methods akin to those used in speech recognition and computational linguistics were used to computationally screen the yeast genome and identify 22 methylation guide snoRNAs, snR50 to snR71. Gene disruptions and other experimental characterization confirmed their methylation guide function. In total, 51 of the 55 ribose methylated sites in yeast ribosomal RNA were assigned to 41 different guide snoRNAs.

The Lsm2-8 complex determines nuclear localization of the spliceosomal U6 snRNA.

Lsm proteins are ubiquitous, multifunctional proteins that are involved in the processing and/or turnover of many, if not all, RNAs in eukaryotes. They generally interact only transiently with their substrate RNAs, in keeping with their likely roles as RNA chaperones. The spliceosomal U6 snRNA is an exception, being stably associated with the Lsm2-8 complex. The U6 snRNA is generally considered to be intrinsically nuclear but the mechanism of its nuclear retention has not been demonstrated, although La protein has been implicated. We show here that the complete Lsm2-8 complex is required for nuclear accumulation of U6 snRNA in yeast. Therefore, just as Sm proteins effect nuclear localization of the other spliceosomal snRNPs, the Lsm proteins mediate U6 snRNP localization except that nuclear retention is the likely mechanism for the U6 snRNP. La protein, which binds only transiently to the nascent U6 transcript, has a smaller, apparently indirect, effect on U6 localization that is compatible with its proposed role as a chaperone in facilitating U6 snRNP assembly.

Proteomic analysis of in vivo-assembled pre-mRNA splicing complexes expands the catalog of participating factors.

Previous compositional studies of pre-mRNA processing complexes have been performed in vitro on synthetic pre-mRNAs containing a single intron. To provide a more comprehensive list of polypeptides associated with the pre-mRNA splicing apparatus, we have determined the composition of the bulk pre-mRNA processing machinery in living cells. We purified endogenous nuclear pre-mRNA processing complexes from human and chicken cells comprising the massive (>200S) supraspliceosomes (a.k.a. polyspliceosomes). As expected, RNA components include a heterogeneous mixture of pre-mRNAs and the five spliceosomal snRNAs. In addition to known pre-mRNA splicing factors, 5' end binding factors, 3' end processing factors, mRNA export factors, hnRNPs and other RNA binding proteins, the protein components identified by mass spectrometry include RNA adenosine deaminases and several novel factors. Intriguingly, our purified supraspliceosomes also contain a number of structural proteins, nucleoporins, chromatin remodeling factors and several novel proteins that were absent from splicing complexes assembled in vitro. These in vivo analyses bring the total number of factors associated with pre-mRNA to well over 300, and represent the most comprehensive analysis of the pre-mRNA processing machinery to date.

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.