Next Article in Journal
Method for B Cell Receptor Enrichment in Malignant B Cells
Previous Article in Journal
Identification and Application of Emerging Biomarkers in Treatment of Non-Small-Cell Lung Cancer: Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unprocessed snRNAs Are a Prognostic Biomarker and Correlate with a Poorer Prognosis in Colorectal Cancer

1
Oncology Area, Center for Biomedical Research of La Rioja—CIBIR, 26006 Logrono, Spain
2
Scientific Computing Group (GRUCACI), University of La Rioja, 26006 Logroño, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2024, 16(13), 2340; https://doi.org/10.3390/cancers16132340
Submission received: 2 May 2024 / Revised: 19 June 2024 / Accepted: 20 June 2024 / Published: 26 June 2024
(This article belongs to the Section Cancer Biomarkers)

Abstract

:

Simple Summary

Integrator-complex deregulation is involved in 8.3% of colorectal cancer cases. Lack of Integrator-complex function, measured by an increased level of unprocessed snRNA, is a prognostic biomarker and correlates with a poorer prognosis in colorectal cancer. Our results show that lack of Integrator-complex function (measured as level of snRNA processing) and not necessarily its expression level correlates with a poorer prognosis in colorectal-cancer patients. Population screening, combined with early typing of tumors, appears to be the most efficient way to increase patient survival.

Abstract

The human Integrator complex is a set of 15 subunits that mediates processing of small nuclear RNAs (snRNAs), and which later participates in splicing messenger RNAs (mRNAs). In addition, it controls the pause and release of RNA polymerase II (RNA pol II) at specific gene promoters in response to growth factors. Mutations in Integrator-complex subunit 6 (INTS6) are associated with different types of tumors. However, the INTS6 gene product does not have a significant prognostic value as a biomarker for tumor progression. Here we show that Integrator-complex deregulation is involved in 8.3% of the colorectal cancer cases diagnosed from the population screen carried out in La Rioja (Spain) from the years 2017 to 2019. Lack of Integrator-complex function, measured by an increased level of unprocessed snRNA, is a prognostic biomarker and correlates with a poorer prognosis in colorectal-cancer patients. The transcriptomic profile of all analyzed colorectal tumors shows a strong alteration of the metabolic state of tumor cells, which compromises standard energy production through mitochondrial respiration, known as the Warburg effect. Furthermore, there is a significant upregulation of genes involved in extracellular matrix organization and collagen rearrangement. This is consistent with tissue reorganization in a growing tumor forming a polyp. Crossing the molecular data generated in this study with the follow-up of patients from population screening indicates that population screening combined with early typing of tumors appears to be the most efficient way to increase patient survival.

Graphical Abstract

1. Introduction

Colon cancer is the second leading cause of cancer-related deaths worldwide. In 2020, more than 1.9 million new cases of colorectal cancer and more than 930,000 deaths due to colorectal cancer were estimated to have occurred worldwide [1].
Colorectal cancer is a type of cancer that affects the colon (large intestine) or rectum. The risk of colorectal cancer increases with age. Most cases affect people over 50 years old. Common symptoms include diarrhea, constipation, blood in the stool, abdominal pain, unexplained weight loss, fatigue, and low iron levels. Many people will not have symptoms in the early stages of the disease. The risk of colorectal cancer can be reduced by eating a healthy diet, staying physically active, not smoking tobacco and limiting alcohol consumption. Regular screenings are crucial for early detection [2]. In fact, incidence rates of colorectal cancer have been decreasing in high-income countries, largely as a result of effective screening programs [3]. The most common screening method is the fecal occult blood test (FOBT). The FOBT detects hidden blood in the stool, which can be an indicator of colorectal cancer or polyps. If blood or abnormal findings are detected in the stool, further diagnostic procedures, such as colonoscopy, are usually recommended to confirm the presence of colorectal cancer or polyps [4].
Diagnostic methods for colorectal cancer include physical examination, imaging (such as abdominal ultrasound, computed tomography scans, and magnetic resonance imaging), examination of the inside of the colon using colonoscopy or sigmoidoscopy, taking a sample of tissue (biopsy) for histopathological examination, and molecular testing to identify specific gene mutations or biomarkers to guide the best treatment option [2].
Treatments for colorectal cancer are based on type and progression of the tumor and the person’s medical history. The primary treatment for early-stage disease (tumor limited to the bowel or local lymph nodes, with no metastatic dissemination to distant organs) is surgical removal of the tumor and nearby lymph nodes. In advanced disease, the primary treatment approach involves chemotherapy and immunotherapy. In some cases, surgery may be recommended for metastatic colorectal cancer to remove tumors that are causing symptoms or blocking the intestine. After treatment, regular follow-up visits and monitoring are essential to detect for any signs of recurrence or new cancer [5,6].
The prognosis for colorectal cancer varies depending on the stage at diagnosis. Early-stage cancers have higher overall survival rates (defined as the percentage of patients who survive the cancer during a specific period of time) than advanced-stage cancers. Timely diagnosis, appropriate treatment, and regular follow-up care are important for improving survival rates and quality of life [5,6].
Over the last two decades, Integrator complex and especially Integrator subunit 6 (INTS6/OMIM 604331) dysregulation have been established to be involved in several types of tumors, such as prostate cancer [7,8,9], esophageal squamous cell carcinoma [10], nasopharyngeal carcinoma [11], hepatocellular carcinoma [12,13] or lung carcinoma [14] among others. This suggests that all these disparately located cancers share a common mechanism of tumor transformation as a consequence of Integrator-complex dysfunction.
The Integrator complex is evolutionary well-conserved from nematodes to humans [15,16], indicating that it plays a central role in the cell biology. It has been described as a set of 15 proteins in humans [15,17,18,19] that binds to RNA pol II and is responsible for snRNA processing during their transcription. Mature snRNAs then form part of the spliceosome for mRNA processing [20]. In addition to this function, Integrator complex also plays a role in pause and release of the RNA pol II at the promoter of specific genes [21,22] and DNA repair [23,24], and in some organisms, INTS6 controls mitochondrial structure and function [25,26]. Given the multiple functions in which the Integrator complex participates, tumors that present mutations in Integrator subunits show global effects on the cell’s biology, altering key functions such as genomic stability and DNA repair, mitochondrial function, general gene expression and snRNA processing. A major consequence of Integrator mutation is the readthrough of the RNA pol II downstream of the snRNA loci [15,18] Consequently, long RNAs derived from unprocessed snRNAs are generated in the cells. In order to delve deeper into the functional characterization of tumors caused by defective Integrator, it is first necessary to establish a rapid and sensitive system for identifying their defects.
In this work, we analyze the impact that Integrator dysfunction has on colorectal tumors. To this end, we measure the level of Integrator subunit 6 mRNA and unprocessed snRNAs in a cohort of colorectal tumor samples and asses its value as a diagnostic and prognostic marker for those tumors. We conclude that snRNA processing may be a good biomarker for patient survival.

2. Materials and Methods

2.1. Patient Selection and Tissue Samples

Biopsy samples were obtained from a cohort of 90 patients showing different degrees of illness. Ten healthy biopsies were selected as control. All samples came from the San Pedro Hospital’s Digestive System Department (La Rioja, Spain). Pathological samples were obtained from lesions with suspected tumor origin observed during diagnostic colonoscopies performed between 2017 and 2019. In all cases, these samples were collected and processed with the informed consent of the patients. The samples from these biopsies were divided for routine histological analysis in the Anatomic Pathology Service and RNA extraction. Control samples came from biopsies of healthy mucosal tissue taken from a region close to the lesion in patients diagnosed with colorectal cancer of this study. This prevents biases related to the patient’s age and sex. Control samples were collected during the same exploration. The indication for having a colonoscopy is the same in all cases (preventive program for colorectal cancer, symptoms associated with the colon or radiological findings suggestive of colorectal pathology). The Mann–Whitney U test yielded a z-score of 2.49326 and a corresponding p value of 0.013. At a significance level of 0.01, the obtained p indicates that the difference in ages between the control group and the patients is not statistically significant. In other words, at a significance level of 0.01, the age distribution in control and patient samples are similar.

2.2. RNA Extraction

RNA extraction was performed with QIAGEN RNeasy® Kit (QIAGEN GmbH, Hilden, Germany) following the instructions for purification of total RNA from animal cells from the manufacturer. Briefly, samples were collected in 350 µL RLT buffer containing β-mercaptoethanol and quickly frozen. Homogenization of the tissue was performed using a conventional Polytron® Rotor-Stator homogenizer (Kinematica AG, Malters, Switzerland) pre-chilled with liquid nitrogen.

2.3. cDNA Synthesis

RNA was treated with DNase to eliminate any DNA contamination. In each sample, a total reaction of 10 µL contained 400 µg RNA, 1 µL RQ1 RNase-Free DNase (Promega Corp., Madison, WI, USA), 1× RQ1 DNase 10X Reaction Buffer and DEPC water. The reaction was incubated for 30 min at 37 °C, stopped by adding 1 µL Stop solution (Promega Corp., Madison, WI, USA) and incubated for 10 min at 65 °C. cDNA synthesis was performed using SuperScript™ III First-Strand Synthesis System for RT-PCR (Invitrogen-Life Technologies Corp., Paisley, UK) with random hexamer primers in a GeneAmp PCR System 9700, Applied Biosystems machine (Thermo Fisher Scientific, Waltham, MA, USA).

2.4. Semiquantitative PCR

INTS6 (Gene ID: 26512. Chromosome 13, NC_000013.11) RNA expression was analyzed using the primers INTS6Fw: 5′-ACAAGGAGACTGGAATAAAAGTCC-3′ and INTS6Rv 5′-TAAACTGCCTTGCACATGCT-3′. A total of 1 µL of cDNA was used for PCR amplification (GeneAmp PCR System 9700, Applied, Foster City, CA, USA) and reactions were optimized to 94 °C for 2 min, 38 amplification cycles at 94 °C for 30 s, 59 °C annealing temperature for 1 min, 72 °C for 1 min, and a final extension of 10 min at 72 °C. U2 processing was analyzed using the following primers in the U2 snRNA pseudogene RNU2-69P. (Gene ID: 106481656. Chromosome 18, NC_000018.10): U269Fw: 5′-CAAGGAGCTGAAAGGCACTGA-3′ and U269Rv: 5′-GACCTGTGCTTTCTGGGGTAG-3′. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH. Gene ID: 2597. Chromosome 12, NC_000012.12) was used as a housekeeping reference gene. GAPDH expression was analyzed using the primers GAPDHFw 5′-AAATCCCATCACCATCTTCC-3′ and GAPDHRv: 5′-GACTCCACGACGTACTCAGC-3′. PCR amplification for U2 and GAPDH were optimized to 94 °C for 2 min, 38 amplification cycles for U2 and 30 for GAPDH at 94 °C for 30 s, 59 °C annealing temperature for 1 min, 72 °C for 15 s, and a final extension of 10 min at 72 °C. A total of 7 µL of the amplified products was resolved on 1.5% agarose gels, visualized by SYBR™-safe staining (Invitrogen-Life Technologies, Carlsbad, CA, USA) and documented using the Bio-Rad Universal Hood II-GelDoc System, 1D Analysis Software: Quantity One 4.6.3, as TIFF digital images.
PCR products were quantified by measuring gel band intensity using ImageJ software (v1.52a) [27]. GelDoc images were stored as TIFF files. On the gel image, a region of a fixed area covering the size of the band on the gel was drawn using the “rectangular selection tool”. This region was moved over the different bands to measure the optical density of each. The optical density of the gel background was measured and subtracted from the value obtained for each band. In addition, we applied a correction factor based on the GAPDH loading control corresponding to each of the samples. The GAPDH band whose intensity was the highest was taken as a reference. Its value was divided by the intensity of each of the other samples to determine the correction factor that had to be applied to each one. Mean and standard deviations of the optical density values of the control samples were obtained to define whether the U2 snRNA and INTS6 mRNA levels in the tumor samples were high (greater than the mean + SD), medium (between the mean ± SD) or low (less than mean − SD).

2.5. Transcriptomic Analysis

Global gene expression profiles from healthy individuals and patients were analyzed at the Genomics and Bioinformatics Platform of the CIBIR, as previously described [16]. Briefly, sequencing libraries were prepared by following the Illumina Stranded Total RNA Prep with Ribo-Zero Plus (Illumina Inc., San Diego, CA, USA) protocol from 100 ng of the total extracted RNA. All libraries were run in a HiSeq1500 PE100 lane. After read quality analysis, trimming, and read mapping against the reference genome, the expression analysis was performed using DESeq2 (v1.34.0) [28] and EdgeR (v3.36.0) [29] algorithms through the SARTools (v1.7.4) [30] pipeline. Genes were scored as differentially expressed genes (DEGs) if their expression, compared to the control, met the following parameters: Padj < 0.05, Up: log2FC > 0, Down: log2FC < 0. The reference genome GRCh38 (patch 105) was obtained from the ENSEMBL database (https://www.ensembl.org/info/data/ftp/index.html, accessed on 3 April 2024). The up- and downregulated genes were analyzed in order to find enriched GO-terms, KEGG and Reactome pathways by using the clusterProfiler R-package (v4.10.1) [31,32]. GO-term enrichment was ordered by their adjusted p value, listing first those which were statistically more significant.

2.6. Patient Monitoring

Patients were monitored with different follow-up guidelines depending on the disease stage. For patients in stage I, surveillance consisted of an annual colonoscopy. If advanced adenoma was detected, the colonoscopy was repeated after one year and if not, it was repeated after three years and subsequently after five years. For patients in stages II and III, surveillance consisted of anamnesis, physical examination and CEA (carcinoembryonic antigen) testing every 3 months for 2 years, and then, every 6 months for 5 years. In addition, these patients underwent annual thoracic–abdominal–pelvic CT (Computed Axial Tomography) for 5 years and an annual colonoscopy. For patients in stage IV, surveillance consisted of anamnesis, physical examination and CEA testing every 3 months for 2 years, and then every 6 months for 5 years. These patients also underwent thoracic–abdominal–pelvic CT every 3–6 months for 5 years and an annual colonoscopy.

2.7. Statistical Analysis

Statistical computations were performed using the R programming language (version 4.0.5) within the RStudio integrated development environment (version 1.4.1106) [33], and specific libraries such as plotly (version 4.9.3) [34], heatmaply (version 1.2.1) [35], ggcorrplot (version 0.1.3) [36] were employed. Correlation of the variables was determined using the Spearman correlation analysis. The statistical significance of the correlation was determined using Student’s t distribution. This analysis allowed us to identify variables that had a strong correlation and quantify their influence on the target variable. All statistical tests were two-tailed and p < 0.05 was considered statistically significant. A Kaplan–Meier risk and survival analysis was then performed to investigate the occurrence of events of interest (in this case, survival) over time. Survival (version 3.2.11) [37], survminer (version 0.4.9) [38], ggplot2 (version 3.3.3) [39] and ggpubr (version 0.4.0) [40] were used as libraries for computing survival analyses and for summarizing and visualizing the results. Finally, univariate and multivariate nonparametric Cox analyses were performed to generate a predictive model for time-to-event wait data. While univariate Cox regression assessed the relationship between any clinicopathological predictor and survival outcome, multivariate Cox regression allowed us to extend this analysis by considering multiple predictors simultaneously, thereby capturing the combined effects of various factors on survival.

3. Results

3.1. Selection of a Set of Samples from Patients with Different Types of Colorectal Tumors

To study the effect that deregulation of the Integrator complex and snRNA processing have on colorectal cancer, we analyzed a cohort of 90 patients showing different degrees of illness, and 10 healthy individuals who served as controls. Patients were selected from the San Pedro Hospital’s Digestive System Department (La Rioja, Spain) based on the appearance of a lesion suspected to be of tumor origin during a diagnostic colonoscopy. The study was run from 2017 to 2019.
To evaluate the degree to which different types of tumors were represented in the selected cohort, we studied several clinicopathological characteristics of the analyzed tumors (Table 1). The group of patients included 54 males (60%) and 36 females (40%); 43 patients (47.8%) were less than 75 years old and 47 (52.2%) were more than 75 years old. This age of 75 years is important, since it is the age from which individuals are no longer included in population screening analyses in Spain. Regarding the reason why a colonoscopy was requested for the patient, 26 patients (28.9%) came from population screening, 46 (51.1%) presented a clinical picture suggestive of neoplasm in the consultation and 18 (20%) presented a previous radiological image in computed axial tomography (CT), nuclear magnetic resonance (MRI) or positron emission tomography (PET/CT) compatible with colon thickening/neoplasm. Further analysis by anatomical location showed that 29 tumors (32.2%) were located in the proximal colon, 35 (38.9%) in the distal colon and 26 (28.9%) in the rectum. Regarding the tumor stage, 17 patients (18.9%) had a grade 1 tumor, 19 patients (21.1%) had a grade 2 tumor, 34 patients (37.8%) had a grade 3 tumor, and 20 patients (22.2%) had a grade 4 tumor. Finally, most of the diagnosed colorectal cancers presented a moderate degree of differentiation at diagnosis (81 of 90, 90%). Therefore, we concluded that the multiple tumors represented in the selected cohort allowed us to study the effect of snRNA processing and INTS6 levels in different types of colorectal tumors.

3.2. Classification of Colorectal Tumors Based on INTS6 Levels and snRNA Processing

To ascertain the impact that Integrator-complex deregulation and snRNA processing have on colorectal cancer, we analyzed their levels in biopsies of colon tumor lesions suggestive of malignancy within the selected cohort of 90 colorectal-tumor patients and compared these to healthy controls.
INTS6 mRNA expression was measured by semiquantitative RT-PCR using internal primers from the gene sequence. snRNA processing was defined by measuring transcription in the region directly downstream of the U2-69 snRNA also using semi-quantitative PCR with a pair of primers in that region (as described in Material and Methods). GAPDH expression was used as an endogenous control for normalizing data (Figure 1). Mean ± standard deviation of expression in healthy control samples was used to define three categories of INTS6 expression and snRNA processing: low, normal, and high expression.
For INTS6 expression, 16 samples (17.8%) from patients showed low levels of expression (lower than the mean minus the standard deviation in healthy individuals). A total of 30 samples (30.3%) showed normal levels (those included between the values of the mean ± the standard deviation in healthy individuals) and 44 samples (48.9%) showed high levels of INTS6 expression (higher than the mean plus the standard deviation in healthy individuals). Similarly, 73 (81.1%) patient tumor samples showed normal snRNA processing (those included between the mean ± the standard deviation values of healthy individuals). A total of 17 patient tumors (18.9%) showed high levels of unprocessed snRNAs (higher than the mean plus the standard deviation of expression in healthy individuals). This group represents those samples which were defective on processing the snRNAs. Interestingly, we did not find any tumors with low levels of unprocessed snRNAs (lower than the mean plus the standard deviation of expression in healthy individuals), which would indicate more efficient snRNA processing than healthy cells. The absence of any tumors in this category indicates that snRNA processing is highly optimized in normal cells and neither a mutation nor deregulation of the Integrator complex can make it more efficient (Table 2).

3.3. Tumor Cluster in Different Transcriptional Categories Based on snRNA Processing and INTS6 mRNA Levels

To assess whether snRNA processing and INTS6 mRNA levels affect the global gene expression profile of tumors, we analyzed the transcriptome of three biopsy samples selected from each of the six categories derived by combining snRNA processing (proper or defective) and INTS6 expression (low, normal or high) data. In addition, we also included in the analysis samples from three biopsies from healthy individuals as controls. Raw sequence data generated in this study are available in the Gene Expression Omnibus (GEO) data repository (Accession number GSE261888).
Gene expression data were normalized by a negative binomial distribution model using DESeq2 and EdgeR implementations. The global-expression profile similarity of all the samples was analyzed using a principal component analysis (PCA) plot (Figure 2A) generated with DESeq2 (v1.34.0) [28]. Very close clustering of the three controls reflects a highly uniform transcriptional signature on biopsy samples from healthy individuals and indicates that colon mucosa biopsy cells express a specific and well-defined set of genes with low variability among individuals. Thus, a colon mucosa biopsy appears to be a good system in which to study gene expression alterations. Differences in gene expression patterns in samples will likely be caused by tumor transformation rather than by individual variability. Strikingly, samples from categories that have high levels of unprocessed U2 plot closer to healthy controls than samples with normal U2 snRNA processing along the principal component 1 (PC1) axis in the PCA. In addition, samples containing high levels of INTS6 mRNA separate from those with low or normal levels of INTS6 along the PC2 axis. These data illustrate a trend in which tumors with defective U2 snRNA processing (and which, therefore, accumulate unprocessed U2) show a more similar gene expression pattern to those of healthy controls than those tumors with normal processing of U2. Among all analyzed samples, those belonging to the category including normal processing of U2 and high levels of INTS6 show the most dissimilar transcriptional profile from those of healthy controls (Figure 2A).
To gain insight into the transcriptional features of each category, we performed a gene ontology analysis (GO)-term enrichment of up- and downregulated gene biological processes (BPs) in each category. The three biopsy samples of each category were considered biological replicas for the statistical analysis of gene expression in that category. Quantitative analysis of differential expression was performed as described in materials and methods in [28,29,30]. The five most significant GO terms (adjusted p < 0.002) in each group are shown in Figure 2B (generated with the function “enrichGO” of the clusterProfiler R-package v4.10.1) [31,32].
No downregulated genes in any of the analyzed groups fall within multiple categories of GO biological processes or cellular compartments. Instead, in all cases, there is a remarkable downregulation of genes involved in processes such as lipid oxidation, aerobic respiration, or metabolism. These genes code for proteins whose location is enriched in mitochondria and cellular apical membrane compartments (Tables S1 and S2). These transcriptomic changes either cause or reflect a dramatic alteration in the metabolic state of tumor cells, compromising standard energy production by mitochondrial respiration, and suggest an alternative production of energy through cytosolic fermentation. This mechanism is known as the Warburg effect [41,42,43]. Interestingly, in addition to the mentioned metabolic changes, samples of tumors with normal U2 snRNA processing and normal INTS6 mRNA levels also show a remarkable downregulation of genes involved in muscle contraction (Figure 2B).
Parallel to this gene downregulation response, gene upregulation also shows a common set of biological processes enriched in all the analyzed groups. This upregulation mainly affects genes involved in extracellular matrix organization and collagen rearrangement. This is consistent with tissue reorganization in a growing tumor forming a polyp. In addition to this general feature, categories showing a more distant transcriptional profile from the control (such as those with normal processing of U2 snRNAs and low or normal levels of INTS6 mRNA) show upregulation of ribosome biogenesis. Finally, tumors with normal processing of U2 snRNAs and a high level of INTS6 mRNA (transcriptionally, the most different from those of the healthy control samples) show a high expression of genes involved in DNA replication and cell division. This could be indicative of a high proliferation ratio in those tumors (Figure 2B, Tables S1 and S2).

3.4. Unprocessed U2 snRNAs Level Correlates with Other Tumor Features

To evaluate the value of snRNA processing and INTS6 mRNA levels as possible biomarkers for the diagnosis and prognosis of colorectal cancer, we performed a Spearman correlation analysis (Figure 3) with the rest of the clinicopathological variables of the tumors. Statistical computations were performed using the R programming language (version 4.0.5) within the RStudio integrated development environment (version 1.4.1106). The Spearman correlation coefficient is a statistical measure used to assess the relationship between two categorical or ordinal variables. To conduct the correlation analysis, the categorical variables were transformed into ordinal scales (Table S3). Specifically, each category within the variables was assigned to a numerical value representing its position in a meaningful order or ranking. This transformation ensured that the data maintained their inherent structure while enabling the use of Spearman’s rank correlation coefficient to effectively explore associations between variables (as described in [44,45]). The Spearman coefficient allows us to capture the direction and strength of the association. Figure 3 displays the Spearman correlation matrix, highlighting significant rho coefficients. The statistical significance of the correlation was determined using the Student’s t distribution, with degrees of freedom equal to the sample length minus two. Coefficients with a p value < 0.05 are shown sharply highlighted to indicate significance, while those with a p value > 0.05 are blurred to indicate insignificance. The complete tables containing all the coefficients and their corresponding p values are provided as Supplementary Material (Table S4). In addition, Table 1 encompasses all analyzed characteristics transformed into ordinal scales, and their association with U2 expression levels stratified into medium and high categories. Each relationship includes the associated number of samples, the total sample size, the p value extracted from the correlation matrix, and the correlation coefficient.
Importantly, we observed a negative correlation between unprocessed U2 snRNA and survival. Lack of snRNA processing and subsequent accumulation of unprocessed snRNAs in tumors significantly correlates with lower survival expectancy in patients (p = 0.029, rho = −0.229). In addition, there was a positive correlation between unprocessed U2 snRNA and INTS6 levels (p = 1 × 10−9, rho = 0.436). High levels of unprocessed U2 snRNA are statistically more abundant in tumors with high levels of INTS6 mRNA (12 of 44, 27.3%) than among those with low or medium levels of INTS6 mRNA (5 of 46, 10.9%). This result suggests the existence of a feedback mechanism, in which INTS6 expression could be induced upon lack of snRNA processing to revert the accumulation of unprocessed snRNAs. Moreover, we evaluated the relationship between tumor grade and levels of unprocessed U2 snRNA. We observed that in the initial stages (I and II) there are fewer patients with high levels of unprocessed U2 snRNA (5 out 36) than in more advanced stages (III and IV) (12 out 54), but this increase is not statistically significant (p = 0.594, rho = 0.050). Finally, we found no correlation between U2 snRNA processing and other variables analyzed such as the reason for requesting the colonoscopy, evaluation of revision colonoscopies, endoscopic pattern presented by the neoplasm, presence of carcinoembryonic antigen (CEA) or age. With respect to the anatomical location of the tumors, high levels of unprocessed U2 snRNA were uniformly distributed in all groups and appeared in five samples from tumors located in the proximal colon, seven in the distal colon, and five from rectal tumors.
Unlike U2 snRNA processing, INTS6 mRNA levels did not show significant correlation with any of the studied clinicopathological features of the tumors. This indicates that although the Int6 protein level has been described as a marker for certain tumors such as hepatocellular carcinomas or breast cancer, [13,46] INTS6 mRNA level cannot be considered a biomarker. This suggests a mechanism in which INTS6 might be post-transcriptionally deregulated in tumor cells, either translationally or at the level of protein stability. Alternatively, INTS6 might be involved in tumors affecting other tissues, but not in colorectal cancer.
Regarding the other tumor features, we found that the presence of tumor markers (p = 0.005, rho = −0.302), high tumor stage (p = 0.0064 × 10−4, rho = −0.530), or advanced progression (p = 0.0087 × 10−9, rho = −0.641) have a very strong negative correlation with survival time, meaning that the more developed the tumor is, the lower the life expectancy of the patient. Along those lines, we found a strong positive correlation between tumor progression and tumor stage (p = 0.0031 × 10−10, rho = 0.7), and presence of tumor markers and tumor stage (p = 0.0068 × 10−5, rho = 0.54). This indicates that high-grade tumors have progressed more and express more tumor markers than tumors that were detected early. Other significant correlations are more related to medical practice than to tumor characteristics. Thus, age shows a negative correlation with the prescription of immunohistochemistry analysis or surgical procedures. The clinical goal is to prevent elderly patients from undergoing aggressive treatments wherein recovery is more harmful than the tumor itself. In these circumstances, palliative care is considered. Finally, performing revision colonoscopies shows a positive correlation with survival since they are only performed in surviving patients. These expected results serve as an internal control for the quality of the statistical analysis.
In conclusion, levels of unprocessed U2snRNAs (but not INTS6 mRNA levels) may be a good biomarker for diagnosis/prognosis of colorectal tumors.

3.5. Unprocessed U2 snRNA Level Impacts Survival

To assess the impact that high levels of unprocessed snRNAs in tumors have on patient survival, we performed a Kaplan–Meier analysis. We divided the patient cohort into two groups (those with normal processing of U2 snRNAs and those with high levels of unprocessed U2 snRNAs) and followed their survival throughout the duration of the study (24 months). We observed that six months after diagnosis, survival in the group that had high levels of unprocessed U2 snRNA was consistently around 20% lower compared to the group with normal U2 snRNA processing. A total of 47% of patients with tumors showing high levels of unprocessed U2 snRNA died during the time the study lasted. In contrast, 28.7% of patients with tumors showing normal levels of unprocessed U2 snRNA died during the same period (Figure 4).
To deconstruct the impact that unprocessed U2 snRNA levels have on survival, not only at the global level, but under different patient clinicopathological conditions, we performed a multivariate Cox analysis. The univariate Cox analysis regarding survival did not indicate statistically significant association with U2 snRNA processing (p = 0.16) (Table S5A). However, conducting a multivariate Cox model to simultaneously observe the influence of unprocessed U2 levels, combined with multiple other clinicopathological conditions of patients revealed a significant impact of U2 processing in combination with age, tumor progression and INTS6 levels on patient survival (p of 0.047, log-rank test p = 7 × 10−5) (Table S5B). By considering the joint contribution of multiple factors, rather than assessing each variable in isolation, multivariate analysis produces more robust and clinically relevant prognostic models. We observed, for instance, that in all groups, high levels of unprocessed U2 snRNA correlated with a worse prognosis (Figure S1). This indicates that snRNA processing per se constitutes a good biomarker for evaluating the prognosis of colorectal tumors, regardless of tumor categorization based on other variables.

3.6. Population Screening Detects Colorectal Tumors in Early Stages

To contrast the social benefit that early detection provides vs. the search for clinicopathological markers that permit doctors to personalize specific tumor treatments, we crossed our results of the impact of snRNA processing in tumors with the population screening carried out in La Rioja (Spain) during the same time period (2017–2019) (Figure 5).
La Rioja is a region in northern Spain, with a population of 316,719 individuals during the years when the study was performed. The colorectal cancer screening program has been deployed in this region since 2010, using fecal occult blood tests. The population screening for early detection of colorectal cancer addressed people aged 50 to 69. (83,111 people) (https://www.ine.es/jaxiT3/Tabla.htm?t=2879&L=0) (accessed on 3 April 2024). Within this target population, an average of 36,526 analyses to detect blood in fecal samples were performed (that is, 44% of the target population, as 56% refused the test for different reasons). A total of 1586 fecal occult blood tests (FOBTs) were positive (4.34% detection rate). As a result of these analyses, 1405 patients with positive results underwent colonoscopies. As expected, a high percentage of patients (88.58%) agreed to undergo a colonoscopy exam once they had the positive blood-in-feces result. Of all the colonoscopies performed, suspicious lesions thought to have a tumor origin were detected in 512 patients. Of these, 416 were advanced adenomas (1.139 % of the people analyzed for blood in their fecal samples) and 96 were diagnosed with colorectal cancer (a rate of 2.63 0/00 of the people analyzed for blood in their fecal samples). Of the patients with colorectal cancer, 66 were grade I or II and 30 were already grade III or IV. These data indicate that colorectal cancer screening is essential to making an early diagnosis. A total of 96 tumors were diagnosed in two years and 68.7% (66/96) of them were in the initial stages of the disease, therefore allowing the cancer to be approached with a curative objective. Finally, 24 of those 96 patients with colorectal cancer came within our study. Of them, 22/24 patients (91.7%) had normal levels of unprocessed U2 snRNA and 2/24 patients (8.3%) had high levels of unprocessed U2 snRNA, which are associated with a worse prognosis (Figure 5). Data were obtained from the San Pedro Hospital’s administration (where all the cancer services in the region are centralized) based on a review of the patients’ medical history. These results suggest that the lack of snRNA processing and the consequent accumulation of non-processed snRNAs could be implicated in the quicker progression of the tumors toward higher and more advanced grades.

4. Discussion

The Integrator complex has been classically described as being primarily involved in the processing of snRNAs and the pause and release of the RNA pol II at the promoter of specific genes [22]. Mutations in Integrator-complex subunits are known to be involved in several tumors, from prostate cancer [7,8,9], esophageal squamous cell carcinoma [10], nasopharyngeal carcinoma [11], and hepatocellular carcinoma [12,13], to lung carcinoma [14]. While INTS6 is highly downregulated in most of the tumors and therefore is suggested to play a role as a tumor suppressor, in other tumors, such as colorectal tumor, INTS6 is markedly upregulated, showing, in this case, an oncogenic activity [47,48].
Lack of snRNA processing and transcriptional profile alteration caused by Integrator-complex deregulation affects the intrinsic molecular mechanism of the tumor transformation of the mutated cells and, therefore, it is independent of its tissular localization. Although some studies indicate that expression of different Integrator-complex subunits, such as INTS6, is altered in specific tumors [13], extensive comparative studies show that INTS6 expression is moderately present in all tissues, being more abundant in the testis, pancreas and liver. In addition, Int6 protein has low cancer specificity. This probably reflects its basic role in cell function, ensuring snRNA maturation and a proper gene expression pattern. Consequently, it is still unclear whether the Int6 protein is a useful tumor biomarker [49]. This apparent paradox (being involved in tumor transformation but having no consistent prognostic value in different tumors) may reflect a tissue-dependent function of Int6 [47,48]. But, additionally, it could be explained at a functional or mechanistic level. In healthy prostate cells, INTS6 induces cell cycle arrest. Mutations leading to a noticeable reduction in the amount of Int6 protein result in tumor progression. Interestingly, re-expression of INTS6 cDNA in those tumor cell lines suppresses their high mitosis rate by regulating cell cycle progression [46]. This means that loss of function of INTS6 is responsible for tumor transformation. In contrast, in colorectal cancer INTS6 expression is significantly increased when compared with normal tissues [48].
Our results show a positive correlation between levels of unprocessed snRNAs and INTS6 mRNA in colorectal tumors. This suggests the existence of a feedback mechanism in which INTS6 expression could be induced upon lack of snRNA processing to revert the accumulation of unprocessed snRNAs. Alternatively, processing of snRNAs could be a finely regulated process and both up- and downregulation of INTS6 might lead to defects in snRNA processing and readthrough of RNA pol II downstream of the snRNA loci.
In this work, we classified a cohort of colorectal tumors by measuring the degree of snRNA processing and found that the lack of snRNA processing has a better prognostic value than the amount of INTS6 mRNA. The higher the level of unprocessed snRNA, the worse the survival expectancy for the patient. Our results strengthen the hypothesis that low Integrator-complex activity, but not necessarily low INTS6 expression, is associated with a poor prognosis in colorectal-cancer patients.
In our cohort we found only tumors that contain normal activity of the Integrator complex (normal processing of snRNAs) or loss of function of the Integrator complex (lack of processing of snRNAs). The fact that we did not find any case in which more efficient processing of snRNAs occurs than in healthy cells indicates that, under normal conditions, the integrator complex is highly efficient in the processing of snRNAs and any mutation or dysregulation unequivocally leads to a loss of function, but never to a gain in function. As expected, other parameters such as high tumor progression or stage show a strong statistical correlation with a poor prognosis.
The transcriptomic changes shown in any of the analyzed tumor groups, either with high or normal processing of snRNAs and low, normal or high levels of INTS6 mRNA, reveal a notable downregulation of genes involved in standard energy production by mitochondrial respiration, known as the Warburg effect [41,42,43]. This effect was detected in all the groups, indicating that it is a major consequence of tumor transformation. However, it was more notable in the group with high levels of unprocessed snRNAs and low levels of INTS6 mRNA, suggesting that the Integrator complex might play a direct role in mitochondrial function. It has been shown that, in organisms such as C. elegans, INTS6 controls mitochondrial structure and function [25,26].
The search for biomarkers and the characterization of the molecular mechanism that underlies the tumors is essential for developing efficacious and directed treatments once the tumor is diagnosed. In our work, the accumulation of unprocessed snRNAs as a malignant tumor marker occurs in 19% of all colorectal cancers. However, specific treatments for tumors that present this specific alteration have not yet been established. As a complementary approach, the population screening program for colorectal tumors allows for 94% detection of the tumor transformations in the adenoma phase or early stages (I and II), which can be safely removed by surgical means. This approach, based on early detection, does not discern the molecular mechanism that underlies the tumor transformation. Thus, it is not efficient for treating tumors that are already in the more advanced stages. Both strategies, searching for biomarkers and understanding the molecular mechanisms that underlie the tumors, in combination with prevention and early diagnosis, are complementary for reducing deaths caused by cancer. However, in the short term, cancer-fighting practices based on prevention and early diagnosis have been overwhelmingly more efficient in terms of the broader population.

5. Conclusions

In conclusion, in this work we report the clinical significance of the Integrator-complex function in tumor severity and the use of unprocessed snRNA as a biomarker for colorectal cancer prognosis. Larger studies are required to validate the usefulness of unprocessed snRNAs as biomarkers in other types of cancer. Population screening combined with early typing of tumors oriented to identify the molecular mechanisms of their tumor transformation, and development of specific therapies, appear to be the most efficient ways to increase patient survival.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers16132340/s1, Figure S1: Graphical representation of the multivariate analysis of the age, tumor progression and INTS6 mRNA value variables with respect to the unprocessed U2 snRNA values; Table S1: Differentially Expressed Genes (DEGs) of each tumor group based on INTS6 levels and snRNA processing compared to healthy individual biopsies; Table S2: Gene Ontology (GO)-term enrichment within the Differentially Expressed Genes (DEGs) of each tumor group based on INTS6 levels and snRNA processing compared to healthy individual biopsies; Table S3: transformation of categorical variables into ordinal scales; Table S4: Values of the rho coefficients in the correlation matrix and p values in paired sample correlation (Student’s t-test distribution); Table S5: Univariate and multivariate Cox analyses.

Author Contributions

Conceptualization, J.C.; methodology, V.E. and E.G.-O.; validation, V.E., E.G.-O. and J.C.; formal analysis, V.E., E.G.-O., R.L., M.d.T. and J.C.; investigation, V.E., C.R.-A., E.G.-O., R.L., M.d.T., Á.M., B.E. and J.C.; data curation, V.E., E.G.-O., R.L. and M.d.T.; writing—original draft preparation, J.C. and V.E.; writing—review and editing, J.C.; supervision, J.C.; funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministerio de Ciencia e Innovación (grant PID2021-127388NB-I00). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

The study was approved by the Institutional Review Board; Medical Research Ethics Committee of La Rioja (CEICLAR, protocol number PI-229, 24 October 2016).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study, or their family.

Data Availability Statement

Raw sequence data generated in this study are available in the Gene Expression Omnibus (GEO) data repository (https://www.ncbi.nlm.nih.gov/geo/) (Accession number GSE261888).

Acknowledgments

We gratefully acknowledge José Luis Vizmanos, José G. Pichel and Christy Esmahan for their valuable comments and discussions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Xi, Y.; Xu, P. Global colorectal cancer burden in 2020 and projections to 2040. Transl. Oncol. 2021, 14, 101174. [Google Scholar] [CrossRef]
  2. Colorectal Cancer. Available online: https://www.who.int/news-room/fact-sheets/detail/colorectal-cancer (accessed on 26 April 2024).
  3. Jodal, H.C.; Helsingen, L.M.; Anderson, J.C.; Lytvyn, L.; Vandvik, P.O.; Emilsson, L. Colorectal Cancer Screening with Faecal Testing, Sigmoidoscopy or Colonoscopy: A Systematic Review and Network Meta-Analysis. BMJ Open 2019, 9, e032773. [Google Scholar] [CrossRef]
  4. Young, G.P.; Symonds, E.L.; Allison, J.E.; Cole, S.R.; Fraser, C.G.; Halloran, S.P.; Kuipers, E.J.; Seaman, H.E. Advances in Fecal Occult Blood Tests: The FIT Revolution. Dig. Dis. Sci. 2015, 60, 609–622. [Google Scholar] [CrossRef]
  5. Argilés, G.; Tabernero, J.; Labianca, R.; Hochhauser, D.; Salazar, R.; Iveson, T.; Laurent-Puig, P.; Quirke, P.; Yoshino, T.; Taieb, J.; et al. Localised Colon Cancer: ESMO Clinical Practice Guidelines for Diagnosis, Treatment and Follow-Up. Ann. Oncol. Off. J. Eur. Soc. Med. Oncol. 2020, 31, 1291–1305. [Google Scholar] [CrossRef]
  6. Cervantes, A.; Adam, R.; Roselló, S.; Arnold, D.; Normanno, N.; Taïeb, J.; Seligmann, J.; de Baere, T.; Osterlund, P.; Yoshino, T.; et al. Metastatic Colorectal Cancer: ESMO Clinical Practice Guideline for Diagnosis, Treatment and Follow-up ☆. Ann. Oncol. 2023, 34, 10–32. [Google Scholar] [CrossRef]
  7. Filleur, S.; Hirsch, J.; Wille, A.; Schön, M.; Sell, C.; Shearer, M.H.; Nelius, T.; Wieland, I. INTS6/DICE1 Inhibits Growth of Human Androgen-Independent Prostate Cancer Cells by Altering the Cell Cycle Profile and Wnt Signaling. Cancer Cell Int. 2009, 9, 28. [Google Scholar] [CrossRef]
  8. Röpke, A.; Buhtz, P.; Böhm, M.; Seger, J.; Wieland, I.; Allhoff, E.P.; Wieacker, P.F. Promoter CpG Hypermethylation and Downregulation of DICE1 Expression in Prostate Cancer. Oncogene 2005, 24, 6667–6675. [Google Scholar] [CrossRef]
  9. Bohm, M.; Maier, C.; Kufer, R.; Ropke, A.; Vogel, W.; Wieland, I. Genetic Variants of DICE1/INTS6 in German Prostate Cancer Families with Linkage to 13q14. Urol. Int. 2015, 95, 386–389. [Google Scholar] [CrossRef]
  10. Li, W.J.; Hu, N.; Su, H.; Wang, C.; Goldstein, A.M.; Wang, Y.; Emmert-Buck, M.R.; Roth, M.J.; Guo, W.J.; Taylor, P.R. Allelic Loss on Chromosome 13q14 and Mutation in Deleted in Cancer 1 Gene in Esophageal Squamous Cell Carcinoma. Oncogene 2003, 22, 314–318. [Google Scholar] [CrossRef]
  11. Lei, T.; Yuen, K.S.; Xu, R.; Tsao, S.W.; Chen, H.; Li, M.; Kok, K.H.; Jin, D.Y. Targeting of DICE1 Tumor Suppressor by Epstein-Barr Virus-Encoded MiR-BART3* MicroRNA in Nasopharyngeal Carcinoma. Int. J. Cancer 2013, 133, 79–87. [Google Scholar] [CrossRef]
  12. Peng, H.; Ishida, M.; Li, L.; Saito, A.; Kamiya, A.; Hamilton, J.P.; Fu, R.; Olaru, A.V.; An, F.; Popescu, I.; et al. Pseudogene INTS6P1 Regulates Its Cognate Gene INTS6 through Competitive Binding of MiR-17-5p in Hepatocellular Carcinoma. Oncotarget 2015, 6, 5666–5677. [Google Scholar] [CrossRef] [PubMed]
  13. Lui, K.Y.; Zhao, H.; Qiu, C.; Li, C.; Zhang, Z.; Peng, H.; Fu, R.; Chen, H.; Lu, M. Integrator Complex Subunit 6 (INTS6) Inhibits Hepatocellular Carcinoma Growth by Wnt Pathway and Serve as a Prognostic Marker. BMC Cancer 2017, 17, 644. [Google Scholar] [CrossRef] [PubMed]
  14. Wieland, I.; Arden, K.C.; Michels, D.; Klein-Hitpass, L.; Böhm, M.; Viars, C.S.; Weidle, U.H. Isolation of DICE1: A Gene Frequently Affected by LOH and Downregulated in Lung Carcinomas. Oncogene 1999, 18, 4530–4537. [Google Scholar] [CrossRef] [PubMed]
  15. Baillat, D.; Hakimi, M.A.; Näär, A.M.; Shilatifard, A.; Cooch, N.; Shiekhattar, R. Integrator, a Multiprotein Mediator of Small Nuclear RNA Processing, Associates with the C-Terminal Repeat of RNA Polymerase II. Cell 2005, 123, 265–276. [Google Scholar] [CrossRef] [PubMed]
  16. Gómez-Orte, E.; Sáenz-Narciso, B.; Zheleva, A.; Ezcurra, B.; de Toro, M.; López, R.; Gastaca, I.; Nilsen, H.; Sacristán, M.P.; Schnabel, R.; et al. Disruption of the Caenorhabditis Elegans Integrator Complex Triggers a Non-Conventional Transcriptional Mechanism beyond SnRNA Genes. PLoS Genet. 2019, 15, e1007981. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, J.; Ezzeddine, N.; Waltenspiel, B.; Albrecht, T.R.; Warren, W.D.; Marzluff, W.F.; Wagner, E.J. An RNAi Screen Identifies Additional Members of the Drosophila Integrator Complex and a Requirement for Cyclin C/Cdk8 in SnRNA 3’-End Formation. RNA 2012, 18, 2148–2156. [Google Scholar] [CrossRef]
  18. Welsh, S.A.; Gardini, A. Genomic Regulation of Transcription and RNA Processing by the Multitasking Integrator Complex. Nat. Rev. Mol. Cell Biol. 2023, 24, 204. [Google Scholar] [CrossRef] [PubMed]
  19. Offley, S.R.; Pfleiderer, M.M.; Zucco, A.; Fraudeau, A.; Welsh, S.A.; Razew, M.; Galej, W.P.; Gardini, A. A Combinatorial Approach to Uncover an Additional Integrator Subunit. Cell Rep. 2023, 42, 112244. [Google Scholar] [CrossRef] [PubMed]
  20. Matera, A.G.; Wang, Z. A Day in the Life of the Spliceosome. Nat. Rev. Mol. Cell Biol. 2014, 15, 108–121. [Google Scholar] [CrossRef]
  21. Gardini, A.; Baillat, D.; Cesaroni, M.; Hu, D.; Marinis, J.M.; Wagner, E.J.; Lazar, M.A.; Shilatifard, A.; Shiekhattar, R. Integrator Regulates Transcriptional Initiation and Pause Release Following Activation. Mol. Cell 2014, 56, 128–139. [Google Scholar] [CrossRef]
  22. Wagner, E.J.; Tong, L.; Adelman, K. Integrator Is a Global Promoter-Proximal Termination Complex. Mol. Cell 2023, 83, 416–427. [Google Scholar] [CrossRef]
  23. Skaar, J.R.; Richard, D.J.; Saraf, A.; Toschi, A.; Bolderson, E.; Florens, L.; Washburn, M.P.; Khanna, K.K.; Pagano, M. INTS3 Controls the HSSB1-Mediated DNA Damage Response. J. Cell Biol. 2009, 187, 25–32. [Google Scholar] [CrossRef]
  24. Zhang, F.; Ma, T.; Yu, X. A Core HSSB1-INTS Complex Participates in the DNA Damage Response. J. Cell Sci. 2013, 126, 4850–4855. [Google Scholar] [CrossRef]
  25. Lee, T.H.; Mun, J.Y.; Han, S.M.; Yoon, G.; Han, S.S.; Koo, H.S. DIC-1 over-Expression Enhances Respiratory Activity in Caenorhabditis Elegans by Promoting Mitochondrial Cristae Formation. Genes Cells 2009, 14, 319–327. [Google Scholar] [CrossRef]
  26. Han, S.M.; Lee, T.H.; Mun, J.Y.; Kim, M.J.; Kritikou, E.A.; Lee, S.J.; Han, S.S.; Hengartner, M.O.; Koo, H.S. Deleted in Cancer 1 (DICE1) Is an Essential Protein Controlling the Topology of the Inner Mitochondrial Membrane in C. Elegans. Development 2006, 133, 3597–3606. [Google Scholar] [CrossRef]
  27. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 Years of Image Analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef]
  28. Anders, S.; Huber, W. Differential Expression Analysis for Sequence Count Data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef]
  29. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef]
  30. Varet, H.; Brillet-Guéguen, L.; Coppée, J.Y.; Dillies, M.A. SARTools: A DESeq2- and EdgeR-Based R Pipeline for Comprehensive Differential Analysis of RNA-Seq Data. PLoS ONE 2016, 11, e0157022. [Google Scholar] [CrossRef] [PubMed]
  31. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. ClusterProfiler 4.0: A Universal Enrichment Tool for Interpreting Omics Data. Innovation 2021, 2, 100141. [Google Scholar] [CrossRef] [PubMed]
  32. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. ClusterProfiler: An R Package for Comparing Biological Themes among Gene Clusters. Omi. A J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  33. RStudio Team. RStudio: Integrated Development for R; RStudio, PBC: Boston, MA, USA, 2020; Available online: http://www.rstudio.com/ (accessed on 3 April 2024).
  34. Sievert, C. Interactive Web-Based Data Visualization with R, Plotly, and Shiny; Chapman and Hall/CRC: Boca Raton, FL, USA, 2020. [Google Scholar]
  35. Galili, T.; O’Callaghan, A.; Sidi, J.; Sievert, C. Heatmaply: An R Package for Creating Interactive Cluster Heatmaps for Online Publishing. Bioinformatics 2018, 34, 1600. [Google Scholar] [CrossRef] [PubMed]
  36. Alboukadel Kassambara Ggcorrplot: Visualization of a Correlation Matrix Using “Ggplot2”. R Package Version 0.1.3. 2019. Available online: https://CRAN.R-Project.org/package=ggcorrplot (accessed on 3 April 2024).
  37. Therneau, T. A Package for Survival Analysis in R_. R Package Version 3.2-11. 2021. Available online: https://CRAN.R-project.org/package=survival (accessed on 3 April 2024).
  38. Kassambara, A.; Kosinski, M.; Biecek, P.; Fabian, S. survminer: Drawing Survival Curves Using ‘ggplot2’. R Package Version 0.4.9. 2021. Available online: https://CRAN.R-project.org/package=survminer (accessed on 3 April 2024).
  39. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  40. Kassambara, A. ggpubr: ‘ggplot2’ Based Publication Ready Plots. R Package Version 0.4.0. 2020. Available online: https://CRAN.R-project.org/package=ggpubr (accessed on 3 April 2024).
  41. Liu, J.; Bai, X.; Zhang, M.; Wu, S.; Xiao, J.; Zeng, X.; Li, Y.; Zhang, Z. Energy Metabolism: A New Target for Gastric Cancer Treatment. Clin. Transl. Oncol. 2024, 26, 338–351. [Google Scholar] [CrossRef] [PubMed]
  42. Schumacher, H.R.; Szekely, I.E.; Patel, S.B.; Fisher, D.R. Leukemic Mitochondria. I. Acute Myeloblastic Leukemia. Am. J. Pathol. 1974, 74, 71–82. [Google Scholar] [PubMed]
  43. Pascale, R.M.; Calvisi, D.F.; Simile, M.M.; Feo, C.F.; Feo, F. The Warburg Effect 97 Years after Its Discovery. Cancers 2020, 12, 2819. [Google Scholar] [CrossRef] [PubMed]
  44. Schober, P.; Schwarte, L.A. Correlation Coefficients: Appropriate Use and Interpretation. Anesth. Analg. 2018, 126, 1763–1768. [Google Scholar] [CrossRef] [PubMed]
  45. Agresti, A. Categorical Data Analysis; Wiley: Hoboken, NJ, USA, 2002. [Google Scholar] [CrossRef]
  46. Gillis, L.D.; Lewis, S.M. Decreased EIF3e/Int6 Expression Causes Epithelial-to-Mesenchymal Transition in Breast Epithelial Cells. Oncogene 2013, 32, 3598–3605. [Google Scholar] [CrossRef] [PubMed]
  47. Expression of INTS6 in Cancer—Summary—The Human Protein Atlas. Available online: https://www.proteinatlas.org/ENSG00000102786-INTS6/pathology (accessed on 5 June 2024).
  48. Ding, X.; Chen, T.; Shi, Q.; Nan, P.; Wang, X.; Xie, D.; Li, J. INTS6 Promotes Colorectal Cancer Progression by Activating of AKT and ERK Signaling. Exp. Cell Res. 2021, 407, 112826. [Google Scholar] [CrossRef]
  49. INTS6 Protein Expression Summary—The Human Protein Atlas. Available online: https://www.proteinatlas.org/ENSG00000102786-INTS6 (accessed on 5 June 2024).
Figure 1. Analysis of INTS6 mRNA levels and U2-69 snRNA processing in colorectal biopsy samples from healthy individuals who served as controls (blue box) vs. colorectal-tumor patients. Numbers correspond to biopsy sample identifiers. Semiquantitative RT-PCR of INTS6 mRNA is shown in the top row. RT-PCR of the region immediately downstream of the U2-69 snRNA, as a measure of snRNA processing, is shown in the middle row. GAPDH mRNA levels (bottom row) are shown as normalization of the loading control. Red asterisks mark levels lower than the mean minus the standard deviation of the corresponding control samples and blue asterisks mark levels greater than the mean minus the standard deviation of the corresponding controls.
Figure 1. Analysis of INTS6 mRNA levels and U2-69 snRNA processing in colorectal biopsy samples from healthy individuals who served as controls (blue box) vs. colorectal-tumor patients. Numbers correspond to biopsy sample identifiers. Semiquantitative RT-PCR of INTS6 mRNA is shown in the top row. RT-PCR of the region immediately downstream of the U2-69 snRNA, as a measure of snRNA processing, is shown in the middle row. GAPDH mRNA levels (bottom row) are shown as normalization of the loading control. Red asterisks mark levels lower than the mean minus the standard deviation of the corresponding control samples and blue asterisks mark levels greater than the mean minus the standard deviation of the corresponding controls.
Cancers 16 02340 g001
Figure 2. Transcriptomic analysis of the tumor categories based on U2 snRNA processing and INTS6 mRNA level. (A) Principal Component Analysis (PCA) plot based on DESeq2 regularized log2 transformation (rlog) data shows clusters according to the similarity of the transcriptional profile. Three samples from each tumor category are shown in the 2D plane, spanned by their first two principal components. Controls corresponding to biopsies from healthy individuals (yellow) grouped very closely, indicating very similar transcriptional profiles. Biopsy samples containing high levels of unprocessed U2 snRNAs (brown and green) show a more similar transcriptional profile to those of controls than samples with normal levels of U2 snRNA processing (orange, purple, and blue). Samples containing high levels of INTS6 mRNA separate from those with low or normal levels of INTS6. Among all analyzed samples, those corresponding to normal processing of U2 and high levels of INTS6 (blue) show the most different transcriptional profile from those of healthy controls. (B) Gene ontology (GO) analysis of up- and downregulated gene biological processes (BP). GO of BP vs. number of genes within each category are shown as color bars, one bar per GO term. Bar length indicates the number of genes belonging to the different GO categories and color illustrates the statistical significance. Only the five most significant GO terms (adjusted p < 0.002, bars color-coded in red) are shown. All defined groups share common features such as downregulation of fatty acid metabolism or mitochondrial respiration and upregulation of genes involved in the organization of the extracellular matrix. However, those groups showing a more distant transcriptional profile from that of the control upregulate genes involved in ribosome biogenesis (highlighted with font color yellow), muscle contraction (highlighted with font color orange) or cell division (highlighted with font color cyan).
Figure 2. Transcriptomic analysis of the tumor categories based on U2 snRNA processing and INTS6 mRNA level. (A) Principal Component Analysis (PCA) plot based on DESeq2 regularized log2 transformation (rlog) data shows clusters according to the similarity of the transcriptional profile. Three samples from each tumor category are shown in the 2D plane, spanned by their first two principal components. Controls corresponding to biopsies from healthy individuals (yellow) grouped very closely, indicating very similar transcriptional profiles. Biopsy samples containing high levels of unprocessed U2 snRNAs (brown and green) show a more similar transcriptional profile to those of controls than samples with normal levels of U2 snRNA processing (orange, purple, and blue). Samples containing high levels of INTS6 mRNA separate from those with low or normal levels of INTS6. Among all analyzed samples, those corresponding to normal processing of U2 and high levels of INTS6 (blue) show the most different transcriptional profile from those of healthy controls. (B) Gene ontology (GO) analysis of up- and downregulated gene biological processes (BP). GO of BP vs. number of genes within each category are shown as color bars, one bar per GO term. Bar length indicates the number of genes belonging to the different GO categories and color illustrates the statistical significance. Only the five most significant GO terms (adjusted p < 0.002, bars color-coded in red) are shown. All defined groups share common features such as downregulation of fatty acid metabolism or mitochondrial respiration and upregulation of genes involved in the organization of the extracellular matrix. However, those groups showing a more distant transcriptional profile from that of the control upregulate genes involved in ribosome biogenesis (highlighted with font color yellow), muscle contraction (highlighted with font color orange) or cell division (highlighted with font color cyan).
Cancers 16 02340 g002
Figure 3. Spearman correlation matrix, highlighting rho significant coefficients. Coefficients with a p < 0.05 are shown sharply highlighted to indicate significance, while those with a p > 0.05 are blurred to indicate insignificance. The coefficients range from −1 to 1, indicating the strength and direction of the monotonic relationship between pairs of variables. Positive values signify a direct monotonic relationship, while negative values indicate an inverse monotonic relationship. Coefficients close to 1 or −1 denote a strong monotonic association, whereas coefficients near 0 suggest a weak or non-monotonic relationship. Each cell in the matrix represents the Spearman correlation coefficient between the corresponding pair of variables, providing insights into their mutual relationships. Thus, there is a negative correlation between levels of unprocessed U2 and patient survival (p = 0.029, rho = −0.229). This indicates that lack of U2 processing in the tumor correlates with lower patient survival. In addition, there is a positive correlation between levels of unprocessed U2 and levels of INTS6 mRNA (p = 0.0, rho = 0.436). This indicates that lack of U2 processing in the tumor correlates with high levels of INTS6 mRNA.
Figure 3. Spearman correlation matrix, highlighting rho significant coefficients. Coefficients with a p < 0.05 are shown sharply highlighted to indicate significance, while those with a p > 0.05 are blurred to indicate insignificance. The coefficients range from −1 to 1, indicating the strength and direction of the monotonic relationship between pairs of variables. Positive values signify a direct monotonic relationship, while negative values indicate an inverse monotonic relationship. Coefficients close to 1 or −1 denote a strong monotonic association, whereas coefficients near 0 suggest a weak or non-monotonic relationship. Each cell in the matrix represents the Spearman correlation coefficient between the corresponding pair of variables, providing insights into their mutual relationships. Thus, there is a negative correlation between levels of unprocessed U2 and patient survival (p = 0.029, rho = −0.229). This indicates that lack of U2 processing in the tumor correlates with lower patient survival. In addition, there is a positive correlation between levels of unprocessed U2 and levels of INTS6 mRNA (p = 0.0, rho = 0.436). This indicates that lack of U2 processing in the tumor correlates with high levels of INTS6 mRNA.
Cancers 16 02340 g003
Figure 4. Kaplan–Meier survival analysis according to U2 snRNA processing: Horizontal axis (x-axis) shows time in months and vertical axis (y-axis) shows the survival probability. Lines represent the survival curves of the two groups into which we have separated the samples: normal unprocessed U2 snRNA levels (blue line) and high unprocessed U2 snRNA levels (red line). At the start of the study (time zero), all the patients studied were alive and, therefore, the probability of survival was 1. After 6 months and onwards, the survival probability was approximately 20% lower for patients with high levels of unprocessed U2 snRNAs compared to those with normal processing of U2 snRNAs.
Figure 4. Kaplan–Meier survival analysis according to U2 snRNA processing: Horizontal axis (x-axis) shows time in months and vertical axis (y-axis) shows the survival probability. Lines represent the survival curves of the two groups into which we have separated the samples: normal unprocessed U2 snRNA levels (blue line) and high unprocessed U2 snRNA levels (red line). At the start of the study (time zero), all the patients studied were alive and, therefore, the probability of survival was 1. After 6 months and onwards, the survival probability was approximately 20% lower for patients with high levels of unprocessed U2 snRNAs compared to those with normal processing of U2 snRNAs.
Cancers 16 02340 g004
Figure 5. Cross-study of the lack of processing of U2 snRNA and population screening for colorectal cancer using fecal occult blood tests carried out in La Rioja (Spain) from 2017 to 2019. During the time period evaluated in the study, 83,111 individuals were selected, of which 36,526 completed the fecal occult blood test screening. A total of 1586 had a positive result. Of the 1405 colonoscopies performed, 96 tumors were diagnosed. From 24 colorectal cancer samples, a quantitatively and qualitatively adequate RNA sample was obtained, and the levels of unprocessed U2 were evaluated. In 22/24 cases, the levels of unprocessed U2 snRNA were normal and 2/24 cases (8.3%) had high levels of unprocessed U2 snRNA.
Figure 5. Cross-study of the lack of processing of U2 snRNA and population screening for colorectal cancer using fecal occult blood tests carried out in La Rioja (Spain) from 2017 to 2019. During the time period evaluated in the study, 83,111 individuals were selected, of which 36,526 completed the fecal occult blood test screening. A total of 1586 had a positive result. Of the 1405 colonoscopies performed, 96 tumors were diagnosed. From 24 colorectal cancer samples, a quantitatively and qualitatively adequate RNA sample was obtained, and the levels of unprocessed U2 were evaluated. In 22/24 cases, the levels of unprocessed U2 snRNA were normal and 2/24 cases (8.3%) had high levels of unprocessed U2 snRNA.
Cancers 16 02340 g005
Table 1. Pathological characteristics associated with U2 expression levels. p-Values and correlation coefficients.
Table 1. Pathological characteristics associated with U2 expression levels. p-Values and correlation coefficients.
VariableUnprocessed U2 snRNATotalp Valuerho Coefficient
NormalHigh
All the cases731790
Localization within colon7317900.951−0.059
024529
128735
221526
Reason for colonoscopy7317900.2460.124
022224
1341246
215318
3202
IHC (Immunohistochemistry)7317900.259−0.135
024630
1391049
2404
3415
4202
Adjuvant therapy7216880.628−0.050
09312
130434
225631
38311
Findings in screening colonoscopy7317900.130−0.132
036844
120727
213215
3404
Age7317900.5910.050
<7023831
>7050959
<7534943
≥7539847
Colonic structures7317900.682−0.012
040848
133942
Anatomical Pathology7317900.0540.213 (*)
0303
1661581
2426
CEA (Carcinoembryonic antigen)7317900.855−0.002
Progression7216880.6530.044
047856
125833
Concomitant diagnosis of other tumors7216880.4840.078
0601373
112315
Survivorship7316890.029−0.229 (*)
Tumor stage7317900.5940.050
015217
116319
226834
316420
INTS6 mRNA7317900.0000.436 (**)
012416
129130
2321244
*. The correlation is significant at the 0.05 level (two-tailed). **. The correlation is significant at the 0.01 level (two-tailed).
Table 2. Classification of the analyzed colorectal tumor samples based on U2 snRNA processing and INTS6 mRNA levels. Numbers correspond to biopsy samples. Samples were grouped in the different categories depending on whether the optical density values of the semi-quantitative RT-PCRs gel bands were high (greater than the mean + SD of the controls), medium (between the mean ± SD of the controls) or low (less than the mean − SD of the controls).
Table 2. Classification of the analyzed colorectal tumor samples based on U2 snRNA processing and INTS6 mRNA levels. Numbers correspond to biopsy samples. Samples were grouped in the different categories depending on whether the optical density values of the semi-quantitative RT-PCRs gel bands were high (greater than the mean + SD of the controls), medium (between the mean ± SD of the controls) or low (less than the mean − SD of the controls).
LOW Levels of INTS6NORMAL Levels of INTS6HIGH Levels of INTS6
NORMAL levels of non-processed U21, 15, 18, 19, 20, 25, 28, 29, 48, 50, 175, 17672, 107, 111, 131, 137, 138, 140, 146, 160, 161, 166, 167, 168, 173, 174, 180, 182, 183, 184, 187, 188, 189, 190, 194, 195, 197, 201, 202, 20323, 40, 45, 54, 55, 58, 60, 62, 67, 68, 71, 73, 79, 80, 101, 105, 115, 118, 121, 122, 124, 129, 141, 164, 169, 171, 172, 179, 191, 193, 199, 200
HIGH levels of non-processed U217, 21, 32, 3311922, 35, 37, 38, 41, 43, 56, 75, 90, 94, 120, 185
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Escrich, V.; Romero-Aranda, C.; López, R.; de Toro, M.; Metola, Á.; Ezcurra, B.; Gómez-Orte, E.; Cabello, J. Unprocessed snRNAs Are a Prognostic Biomarker and Correlate with a Poorer Prognosis in Colorectal Cancer. Cancers 2024, 16, 2340. https://doi.org/10.3390/cancers16132340

AMA Style

Escrich V, Romero-Aranda C, López R, de Toro M, Metola Á, Ezcurra B, Gómez-Orte E, Cabello J. Unprocessed snRNAs Are a Prognostic Biomarker and Correlate with a Poorer Prognosis in Colorectal Cancer. Cancers. 2024; 16(13):2340. https://doi.org/10.3390/cancers16132340

Chicago/Turabian Style

Escrich, Víctor, Cristina Romero-Aranda, Rosario López, María de Toro, Ángela Metola, Begoña Ezcurra, Eva Gómez-Orte, and Juan Cabello. 2024. "Unprocessed snRNAs Are a Prognostic Biomarker and Correlate with a Poorer Prognosis in Colorectal Cancer" Cancers 16, no. 13: 2340. https://doi.org/10.3390/cancers16132340

APA Style

Escrich, V., Romero-Aranda, C., López, R., de Toro, M., Metola, Á., Ezcurra, B., Gómez-Orte, E., & Cabello, J. (2024). Unprocessed snRNAs Are a Prognostic Biomarker and Correlate with a Poorer Prognosis in Colorectal Cancer. Cancers, 16(13), 2340. https://doi.org/10.3390/cancers16132340

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop