Systematic Profiling of Alternative Splicing for Sarcoma Patients Reveals Novel Prognostic Biomarkers Associated with Tumor Microenvironment and Immune Cells

Identification of Prognostic lncRNA Related to the Immune Microenvironment of Soft Tissue Sarcoma

Background. Soft tissue sarcoma is a malignant tumor with high degree of malignancy and poor prognosis, originating from mesenchymal tissue. Long noncoding RNAs (lncRNAs) are involved in various biological and pathological processes in the body. They perform preprocessing, splicing, transport, degradation, and translation of mRNA to achieve posttranscriptional level regulation, resulting in the occurrence, invasion, and metastasis of tumors. Therefore, they are highly relevant with regard to early diagnoses and as prognostic indicators. Objective. The objective of the present study was to identify immune microenvironment-related lncRNAs that can be used to predict soft tissue sarcomas. Methods. Clinical data and follow-up data were obtained from the cBioPortal database, and RNA sequencing data used for the model structure can be accessed from The Cancer Genome Atlas (TCGA) database. LncRNAs were screened by differential expression analysis and coexpression analysis. The Cox regression model and Kaplan–Meier analysis were used to study the association between lncRNAs and soft tissue sarcoma prognosis in the immune microenvironment. Unsupervised cluster analysis was then completed to discover the impact of screening lncRNAs on disease. We constructed an mRNA-lncRNA network by Cytoscape software. Finally, qRT-PCR was used to verify the difference in the expression of the lncRNAs in normal cells and sarcoma cells. Results. Unsupervised cluster analysis revealed that the 210 lncRNAs screened showed strong correlation with the tumor immune microenvironment. Two signatures containing seven and five lncRNAs related to the tumor microenvironment were constructed and used to predict overall survival (OS) and disease-free survival (DFS). The Kaplan–Meier (K-M) survival curve showed that the prognoses of patients in the high-risk and low-risk groups differed significantly, and the prognosis associated with the low-risk group was better than that associated with the high-risk group. Two nomograms with predictive capabilities were established. qRT-PCR results showed that the expression of AC108134.3 and AL031717.1 was significantly different in normal and sarcoma cells. Conclusion. In summary, the experimental results showed that lncrnA associated with immune microenvironment was related to tumor, which may provide a new idea for immunotherapy of STS.

Prognostic Signatures of Alternative Splicing Events in Esophageal Carcinoma Based on TCGA Splice-Seq Data

An alternative splicing (AS) event is a highly complex process that plays an essential role in post-transcriptional gene expression. Several studies have suggested that abnormal AS events were the primary element in the pathological process of cancer. However, few works are dedicated to the study of AS events in esophageal carcinoma (EC). In the present study, clinical information and RNA-seq data of EC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The percent spliced in (PSI) values of AS events were acquired from the TCGA Splice-seq. A total of 183 EC patients were enrolled in this study, and 2,212 AS events were found significantly associated with the overall survival of these patients by univariate Cox regression analysis. The prognostic signatures based on AS events were built by multivariate Cox analysis. Receiver operating characteristic (ROC) curves displayed that the area under the curve (AUC) of the following prognostic signatures, including exon skip (ES), alternate terminator (AT), alternate acceptor site (AA), alternate promoter (AP), alternate donor site (AD), retained intron (RI), and total events, was greater than 0.8, suggesting that these seven signatures had valuable prognosis prediction capacity. Finally, the risk score of prognostic signatures was indicated as an independent risk factor of survival. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to explore the function of splicing factors (SFs) that were associated with AS events. Also, the interactive network between AS events and SFs identified several hub genes and AS events which need further study. This was a comprehensive study that explored prognosis-related AS events and established valuable prognosis signatures in EC patients. The network of interactions between AS events and SFs might offer novel insights into the fundamental mechanisms of tumorigenesis and progression of EC.

Interplay within tumor microenvironment orchestrates neoplastic RNA metabolism and transcriptome diversity

The heterogeneous population of cancer cells within a tumor mass interacts intricately with the multifaceted aspects of the surrounding microenvironment. The reciprocal crosstalk between cancer cells and the tumor microenvironment (TME) shapes the cancer pathophysiome in a way that renders it uniquely suited for immune tolerance, angiogenesis, metastasis, and therapy resistance. This dynamic interaction involves a dramatic reconstruction of the transcriptomic landscape of tumors by altering the synthesis, modifications, stability, and processing of gene readouts. In this review, we categorically evaluate the influence of TME components, encompassing a myriad of resident and infiltrating cells, signaling molecules, extracellular vesicles, extracellular matrix, and blood vessels, in orchestrating the cancer-specific metabolism and diversity of both mRNA and noncoding RNA, including micro RNA, long noncoding RNA, circular RNA among others. We also highlight the transcriptomic adaptations in response to the physicochemical idiosyncrasies of TME, which include tumor hypoxia, extracellular acidosis, and osmotic stress. Finally, we provide a nuanced analysis of existing and prospective therapeutics targeting TME to ameliorate cancer-associated RNA metabolism, consequently thwarting the cancer progression. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Turnover and Surveillance > Regulation of RNA Stability RNA in Disease and Development > RNA in Disease.

TGFβ1: An Indicator for Tumor Immune Microenvironment of Colon Cancer From a Comprehensive Analysis of TCGA

Background

Tumor microenvironment (TME) and tumor-infiltrating immune cells (TICs) greatly participate in the genesis and development of colon cancer (CC). However, there is little research exploring the dynamic modulation of TME.

Methods

We analyzed the proportion of immune/stromal component and TICs in the TME of 473 CC samples and 41 normal samples from The Cancer Genome Atlas (TCGA) database through ESTIMATE and CIBERSORT algorithms. Correlation analysis was conducted to evaluate the association between immune/stromal component in the TME and clinicopathological characteristics of CC patients. The difference analysis was performed to obtain the differentially expressed genes (DEGs). These DEGs were further analyzed by GO and KEGG enrichment analyses, PPI network, and COX regression analysis. Transforming growth factor β1 (TGFβ1) was finally overlapped from the above analysis. Paired analysis and GSEA were carried out to understand the role of TGFβ1 in colon cancer. The intersection between the difference analysis and correlation analysis was conducted to learn the association between TGFβ1 and TICs.

Results

Our results showed that the immune component in the TME was negatively related with the stages of CC. GO and KEGG enrichment analysis revealed that 1,110 DEGs obtained from the difference analysis were mainly enriched in immune-related activities. The intersection analysis between PPI network and COX regression analysis indicated that TGFβ1 was significantly associated with the communication of genes in the PPI network and the survival of CC patients. In addition, TGFβ1 was up-regulated in the tumor samples and significantly related with poor prognosis of CC patients. Further GSEA suggested that genes in the TGFβ1 up-regulated group were enriched in immune-related activities and the function of TGFβ1 might depend on the communications with TICs, including T cells CD4 naïve and T cells regulatory.

Conclusion

The expression of TGFβ1 might be an indicator for the tumor immune microenvironment of CC and serve as a prognostic factor. Drugs targeting TGFβ1 might be a potential immunotherapy for CC patients in the future.

Development and validation of an immune gene-set based prognostic signature for soft tissue sarcoma

Background

Sarcomas is a group of heterogeneous malignant tumors originated from mesenchymal tissue and different types of sarcomas have disparate outcomes. The present study aims to identify the prognostic value of immune-related genes (IRGs) in sarcoma and establish a prognostic signature based on IRGs.

Methods

We collected the expression profile and clinical information of 255 soft tissue sarcoma samples from The Cancer Genome Atlas (TCGA) database and 2498 IRGs from the ImmPort database. The LASSO algorithm and Cox regression analysis were used to identify the best candidate genes and construct a signature. The prognostic ability of the signature was evaluated by ROC curves and Kaplan-Meier survival curves and validated in an independent cohort. Besides, a nomogram based on the IRGs and independent prognostic clinical variables was developed.

Results

A total of 19 IRGs were incorporated into the signature. In the training cohort, the AUC values of signature at 1-, 2-, and 3-years were 0.938, 0.937 and 0.935, respectively. The Kaplan-Meier survival curve indicated that high-risk patients were significantly worse prognosis (P < 0.001). In the validation cohort, the AUC values of signature at 1-, 2-, and 3-years were 0.730, 0.717 and 0.647, respectively. The Kaplan-Meier survival curve also showed significant distinct survival outcome between two risk groups. Furthermore, a nomogram based on the signature and four prognostic variables showed great accuracy in whole sarcoma patients and subgroup analyses. More importantly, the results of the TF regulatory network and immune infiltration analysis revealed the potential molecular mechanism of IRGs.

Conclusions

In general, we identified and validated an IRG-based signature, which can be used as an independent prognostic signature in evaluating the prognosis of sarcoma patients and provide potential novel immunotherapy targets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-07852-2.

Alternative splicing: Human disease and quantitative analysis from high-throughput sequencing

Alternative splicing contributes to the majority of protein diversity in higher eukaryotes by allowing one gene to generate multiple distinct protein isoforms. It adds another regulation layer of gene expression. Up to 95% of human multi-exon genes undergo alternative splicing to encode proteins with different functions. Moreover, around 15% of human hereditary diseases and cancers are associated with alternative splicing. Regulation of alternative splicing is attributed to a set of delicate machineries interacting with each other in aid of important biological processes such as cell development and differentiation. Given the importance of alternative splicing events, their accurate mapping and quantification are paramount for downstream analysis, especially for associating disease with alternative splicing. However, deriving accurate isoform expression from high-throughput RNA-seq data remains a challenging task. In this mini-review, we aim to illustrate I) mechanisms and regulation of alternative splicing, II) alternative splicing associated human disease, III) computational tools for the quantification of isoforms and alternative splicing from RNA-seq.

Clusterization in acute myeloid leukemia based on prognostic alternative splicing signature to reveal the clinical characteristics in the bone marrow microenvironment

Background

Alternative splicing (AS), a crucial post-transcriptional regulatory mechanism in expanding the coding capacities of genomes and increasing the diversity of proteins, still faces various challenges in the splicing regulation mechanism of acute myeloid leukemia (AML) and microenvironmental changes.

Results

A total of 27,833 AS events were detected in 8337 genes in 178 AML patients, with exon skip being the predominant type. Approximately 11% of the AS events were significantly related to prognosis, and the prediction models based on various events demonstrated high classification efficiencies. Splicing factors correlation networks further altered the diversity of AS events through epigenetic regulation and clarified the potential mechanism of the splicing pathway. Unsupervised cluster analysis revealed significant correlations between AS and immune features, molecular mutations, immune checkpoints and clinical outcome. The results suggested that AS clusters could be used to identify patient subgroups with different survival outcomes in AML, among which C1 was both associated with good outcome in overall survival. Interestingly, C1 was associated with lower immune scores compared with C2 and C3, and favorable-risk cytogenetics was rarely distributed in C2, but much more common in C1.

Conclusions

This study revealed a comprehensive landscape of AS events, and provides new insight into molecular targeted therapy and immunotherapy strategy for AML.