SUVA: splicing site usage variation analysis from RNA-seq data reveals highly conserved complex splicing biomarkers in liver cancer

Identification of Regulatory RNA-Binding Proteins Associated with Immune Infiltration in Laryngeal Squamous Cell Carcinoma

We analyzed bulk RNA sequencing and single-cell RNA sequencing (scRNA-seq) data to identify alternative splicing (AS) events and regulatory RNA-binding proteins (RBPs) associated with immune infiltration in human laryngeal squamous cell carcinoma (LSCC). Whole-transcriptome sequencing data of 20 human laryngeal cancer and paracancerous tissues were downloaded from the Gene Expression Omnibus public database, using newly published splicing-site usage variation analysis software to obtain highly conserved regulated AS (RAS) events, and scientific reverse convolution algorithm analysis was used to identify significantly different immune cells and perform a correlation analysis between the two. The software package edgeR was used to identify differentially expressed RBPs and the immune infiltration-related LSCC-RAS they may regulate. Finally, we present the expression profiles and survival curves of 117 human laryngeal cancer samples from The Cancer Genome Atlas dataset for the identified RBPs and LSCC-RAS. We also downloaded the gene set enrichment 150321 scRNA-seq data for two human LSCC tissue samples. The RBP expression pattern and the expression of prophase RBP genes were analyzed in different LSCC cell populations. RNA-binding motif protein 47 (RBM47) and filamin A, as well as the RBP-RAS events that were screened in both the fibulin 2 and fibronectin 1 genes, were all significantly associated with the prognosis, and the RBM47 gene was upregulated in myeloid cells. Because the prognosis was significantly associated with two RBP regulators and two LSCC-RAS events, they may be critical regulators of immune cell survival during laryngeal cancer progression, and RBM47 may regulate macrophage-associated AS and affect immunity.

RNA Sequencing Analyses Reveal the Potential Anti-Inflammatory Mechanisms of Acacetin Against ODG/R Injuries in Microglia

Background

Acacetin is a natural flavonoid known for its anti-tumor, antioxidant, and anti-inflammatory properties. Our previous studies have shown its protective effects against cerebral ischemia-reperfusion injury (IRI), but the underlying molecular mechanisms remain unclear.

Purpose

The study delves into acacetin's mechanism in mitigating cerebral IRI, with a focus on transcriptomic insights.

Methods

We established the oxygen-glucose deprivation/re-oxygenation (OGD/R) model in BV2 microglia, treating them with 10μM acacetin. Then we assessed cell proliferation using CCK-8 and measured Lactate Dehydrogenase (LDH) release. High-throughput RNA sequencing (RNA-seq) underpinned the analysis of differentially expressed genes (DEGs) and long non-coding RNAs (lncRNAs), functional enrichment, and alternative splicing events (ASEs), validated by reverse transcription-quantitative polymerase chain reaction (RT-qPCR).

Results

OGD/R injury significantly impaired cell proliferation and increased LDH release, effects mitigated by acacetin. RNA-seq identified 2148 upregulated and 2135 downregulated DEGs post-OGD/R. In contrast, the acacetin-treated group showed 248 upregulated and 240 downregulated DEGs compared to the OGD/R group. All DEGs were enriched in both Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Overlapping analysis indicated that acacetin treatment reversed the expression of 203 genes affected by OGD/R, including inflammation-related genes such as Isg15, Fcgr1, Il1b, and Parp12. Moreover, the oxidative stress-related gene, Mt2, was downregulated post-OGD/R but upregulated following acacetin treatment. We further found that OGD/R and acacetin treatment could modulate gene splicing events, impacting cell apoptosis or inflammatory responses, such as the A3SS splicing event in the Trim47 gene. RNA-seq also highlighted differential expression of numerous lncRNAs, particularly the upregulation of lncRNA Rmrp and Terc post-OGD/R and their subsequent downregulation post-acacetin treatment. These lncRNAs might regulate cell proliferation through mediating target gene expressions. RT-qPCR validation confirmed these findings.

Conclusion

Significant upregulation of genes and ASEs linked to oxidative stress and inflammatory response is observed in cerebral IRI. Acacetin intervention reverses these effects, highlighting its mechanism in alleviating the injury by modulating gene expression and splicing events.

Genome-wide identification and analysis of epithelial-mesenchymal transition-related RNA-binding proteins and alternative splicing in a human breast cancer cell line

Exploring the mechanism of breast cancer metastasis and searching for new drug therapeutic targets are still the focuses of current research. RNA-binding proteins (RBPs) may affect breast cancer metastasis by regulating alternative splicing (AS) during epithelial-mesenchymal transition (EMT). We hypothesised that during EMT development in breast cancer cells, the expression level of RBPs and the gene AS pattern in the cell were significantly changed on a genome-wide scale. Using GEO database, this study identified differentially expressed RBPs and differential AS events at different stages of EMT in breast cancer cells. By establishing the correlation network of differential RBPs and differential AS events, we found that RBM47, PCBP3, FRG1, SRP72, RBMS3 and other RBPs may regulate the AS of ITGA6, ADGRE5, TNC, COL6A3 and other cell adhesion genes. By further analysing above EMT-related RBPs and AS in breast cancer tissues in TCGA, it was found that the expression levels of ADAT2, C2orf15, SRP72, PAICS, RBMS3, APOBEC3G, NOA1, ACO1 and the AS of TNC and COL6A3 were significantly correlated with the prognosis of breast cancer patients. The expression levels of all 8 RBPs were significantly different in breast cancer tissues without metastasis compared with normal breast tissues. Conclusively, eight RBPs such as RBMS3 and AS of TNC and COL6A3 could be used as predictors of breast cancer prognosis. These findings need to be further explored as possible targets for breast cancer treatment.

RNA-binding proteins potentially regulate the alternative splicing of apoptotic genes during knee osteoarthritis progression

BACKGROUND

Alternative splicing (AS) is a principal mode of genetic regulation and one of the most widely used mechanisms to generate structurally and functionally distinct mRNA and protein variants. Dysregulation of AS may result in aberrant transcription and protein products, leading to the emergence of human diseases. Although considered important for regulating gene expression, genome-wide AS dysregulation, underlying mechanisms, and clinical relevance in knee osteoarthritis (OA) remain unelucidated. Therefore, in this study, we elucidated and validated AS events and their regulatory mechanisms during OA progression.

RESULTS

In this study, we identified differentially expressed genes between human OA and healthy meniscus samples. Among them, the OA-associated genes were primarily enriched in biological pathways such as extracellular matrix organization and ossification. The predominant OA-associated regulated AS (RAS) events were found to be involved in apoptosis during OA development. The expression of the apoptosis-related gene BCL2L13, XAF1, and NF2 were significantly different between OA and healthy meniscus samples. The construction of a covariation network of RNA-binding proteins (RBPs) and RAS genes revealed that differentially expressed RBP genes LAMA2 and CUL4B may regulate the apoptotic genes XAF1 and BCL2L13 to undergo AS events during OA progression. Finally, RT-qPCR revealed that CUL4B expression was significantly higher in OA meniscus samples than in normal controls and that the AS ratio of XAF1 was significantly different between control and OA samples; these findings were consistent with their expected expression and regulatory relationships.

CONCLUSIONS

Differentially expressed RBPs may regulate the AS of apoptotic genes during knee OA progression. XAF1 and its regulator, CUL4B, may serve as novel biomarkers and potential therapeutic targets for this disease.

Genome-wide identification of dysregulated alternative splicing and RNA-binding proteins involved in atopic dermatitis

Objectives: We explored the role and molecular mechanisms of RNA-binding proteins (RBPs) and their regulated alternative splicing events (RASEs) in the pathogenesis of atopic dermatitis (AD). Methods: We downloaded RNA-seq data (GSE121212) from 10 healthy control skin samples (healthy, Ctrl), 10 non-lesional skin samples with AD damage (non-lesional, NL), and 10 lesional skin samples with AD damage (lesional, LS). We performed the analysis of differentially expressed genes (DEGs), differentially expressed RBPs (DE-RBPs), alternative splicing (AS), functional enrichment, the co-expression of RBPs and RASEs, and quantitative polymerase chain reaction (qPCR). Results: We identified 60 DE-RBP genes by intersecting 2141 RBP genes from existing reports with overall 2697 DEGs. Most of the DE-RBP genes were found to be upregulated in the AD LS group and related to immune and apoptosis pathways. We observed different ASEs and RASEs among the healthy, AD NL, and AD LS groups. In particular, alt3p and alt5p were the main ASEs and RASEs in AD NL and AD LS groups, compared to the healthy group. Furthermore, we constructed co-expression networks of DE-RBPs and RAS, with particular enrichment in biological pathways including cytoskeleton organization, inflammation, and immunity. Subsequently, we selected seven genes that are commonly present in these three pathways to assess their expression levels in the peripheral blood mononuclear cells (PBMCs) from both healthy individuals and AD patients. The results demonstrated the upregulation of four genes (IFI16, S100A9, PKM, and ENO1) in the PBMCs of AD patients, which is highly consistent with DE-RBP genes analysis. Finally, we selected four RAS genes regulated by RBPs that were related to immune pathways and examined their RASEs in PBMCs from both AD patients and healthy controls. The results revealed an increased percentage of RASEs in the DDX60 gene in AD, which is highly consistent with AS analysis. Conclusion: Dysregulated RBPs and their associated RASEs may have a significant regulatory role in the development of AD and could be potential therapeutic targets in the future.

Identification of alternative splicing and RNA-binding proteins involved in myocardial ischemia-reperfusion injury

Alternative splicing (AS) and RNA-binding proteins (RBPs) have been implicated in various cardiovascular diseases. Yet, a comprehensive understanding of their role in myocardial ischemia-reperfusion injury (MIRI) remains elusive. We aimed to identify potential therapeutic targets for MIRI by studying genome-wide changes in AS events and RBPs. We analyzed RNA-seq data from ischemia-reperfusion mouse models and the control group from the GSE130217 data set using Splicing Site Usage Variation Analysis software. We identified 28 regulated alternative splicing events (RASEs) and 47 differentially expressed RBP (DE-RBP) genes in MIRI. Most variable splicing events were involved in cassette exon, alternative 5' splice, alternative 3' splice, and retained intron types. Gene Ontology and Kyoto Encyclopedia of Genes (KOBAS 2.0 server) and Genomes pathway enrichment analyses showed that the differentially expressed variable splicing and RBP genes were mainly enriched in pathways related to myocardial function. The RBP-RASE network demonstrated a common variance relationship between DE-RBPs and RASEs, indicating that RBPs regulate variable shear events in MIRI. This study systematically identified important alterations in RASEs and RBPs in MIRI, expanding our understanding of the underlying pathogenesis of MIRI.

Genome-wide identification and functional analysis of dysregulated alternative splicing profiles in sepsis

BACKGROUND

An increasing body of evidence now shows that the long-term mortality of patients with sepsis are associated with various sepsis-related immune cell defects. Alternative splicing (AS), as a sepsis-related immune cell defect, is considered as a potential immunomodulatory therapy target to improve patient outcomes. However, our understanding of the role AS plays in sepsis is currently insufficient.

AIM

This study investigated possible associations between AS and the gene regulatory networks affecting immune cells. We also investigated apoptosis and AS functionality in sepsis pathophysiology.

METHODS

In this study, we assessed publicly available mRNA-seq data that was obtained from the NCBI GEO dataset (GSE154918), which included a healthy group (HLTY), a mild infection group (INF1), asepsis group (Seps), and a septic shock group (Shock). A total of 79 samples (excluding significant outliers) were identified by a poly-A capture method to generate RNA-seq data. The variable splicing events and highly correlated RNA binding protein (RBP) genes in each group were then systematically analyzed.

RESULTS

For the first time, we used systematic RNA-seq analysis of sepsis-related AS and identified 1505 variable AS events that differed significantly (p <= 0.01) across the four groups. In the sepsis group, the genes related to significant AS events, such as, SHISA5 and IFI27, were mostly enriched in the cell apoptosis pathway. Furthermore, we identified differential splicing patterns within each of the four groups. Significant differences in the expression of RNA Binding Protein(RBP) genes were observed between the control group and the sepsis group. RBP gene expression was highly correlated with variant splicing events in sepsis, as determined by co-expression analysis; The expression of DDX24, CBFA2T2, NOP, ILF3, DNMT1, FTO, PPRC1, NOLC1 RBPs were significant reduced in sepsis compared to the healthy group. Finally, we constructed an RBP-AS functional network.

CONCLUSION

Analysis indicated that the RBP-AS functional network serves as a critical post-transcriptional mechanism that regulates the development of sepsis. AS dysregulation is associated with alterations in the regulatory gene expression network that is involved in sepsis. Therefore, the RBP-AS expression network could be useful in refining biomarker predictions in the development of new therapeutic targets for the pathogenesis of sepsis.

Genome-wide identification of aberrant alternative splicing and RNA-binding protein regulators in acute myeloid leukaemia which may contribute to immune microenvironment remodelling

Acute myeloid leukaemia (AML) is one of the most lethal cancers of the haematopoietic system with a poorly understood aetiology. Recent studies have shown that aberrant alternative splicing (AS) and a (RBP) regulators are highly associated with the pathogenesis of AML. This study presents an overview of the abnormal AS and differential expression of RNA-binding proteins (RBPs) in AML and further highlights their close relation to the remodelling of the immune microenvironment in AML patients. An in-depth understanding of the regulatory mechanism underlying AML will contribute to the future development of strategies for the prevention, diagnosis and therapy of AML and thus improve the overall survival of patients with AML.

The regulation of expression and splicing of transcription factors are related to the muscle damage caused by radiation in tree shrews

Radiotherapy-induced muscle injury (RIMI) is a major complication of radiotherapy for nasopharyngeal carcinoma. Transcription factor (TF) expression and alternative splicing are crucial events in transcriptional and posttranscriptional regulation, respectively, and are known to be involved in key signaling pathways contributing to a variety of human disorders, including radiation injury. To investigate the TFs and alternative splicing events involved in RIMI, we constructed a tree shrew model as described previously in which the RIMI group received 20 Gy of irradiation on the tensor veli palatini (TVP) muscles. The irradiated muscles were evaluated by RNA sequencing (RNA-seq) 6 months later, and the results compared with those for normal TVP muscles. The alt5p and alt3p events were the two main types of differentially regulated alternative splicing events (RASEs) identified via the Splice sites Usage Variation Analysis (SUVA) software, and these RASEs were highly conserved in RIMI. According to functional enrichment analysis, the differentially RASEs were primarily enriched in pathways related to transcriptional regulation. Furthermore, we identified 16 alternative splicing TFs (ASTFs) in ASTF-differentially expressed gene (DEG) networks based on co-expression analysis, and the regulatory networks were chiefly enriched in pathways linked to cell proliferation and differentiation. This study revealed that RASEs and ASTF-DEG networks may both play important regulatory roles in gene expression network alteration in RIMI. Future studies on the targeting mechanisms and early interventions directed at RASEs and ASTF-DEG networks may aid in the treatment of RIMI.

Genome-wide analysis of dysregulated RNA-binding proteins and alternative splicing genes in keloid

Introduction: The pathogenesis of keloids remains unclear. Methods: In this study, we analyzed RNA-Seq data (GSE113619) of the local skin tissue of 8 keloid-prone individuals (KPI) and 6 healthy controls (HC) before and 42 days after trauma from the gene expression omnibus (GEO) database. The differential alternative splicing (AS) events associated with trauma healing between KPIs and HCs were identifified, and their functional differences were analyzed by gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathways. The co-expression relationship of differentially alternative splicing genes and differentially expressed RNA binding proteins (RBPs) was established subsequently. Results: A total of 674 differential AS events between the KD42 and the KD0 and 378 differential AS events between the HD42 and the HD0 were discovered. Notably, most of the differential genes related to keloids are enriched in actin, microtubule cells, and cortical actin cytoskeletal tissue pathway. We observed a signifificant association between AS genes (EPB41, TPM1, NF2, PARD3) and trauma healing in KPIs and HCs. We also found that the differential expression of healthy controls-specifific trauma healing-related RBPs (TKT, FDPS, SAMHD1) may affect the response of HCs to trauma healing by regulating the AS of downstream trauma healing-related genes such as DCN and DST. In contrast, KPIs also has specifific differential expression of trauma healing related RBPs (S100A9, HspB1, LIMA1, FBL), which may affect the healing response of KPIs to trauma by regulating the AS of downstream trauma healing-related genes such as FN1 and TPM1. Discussion: Our results were innovative in revealing early wound healing-related genes (EPB41, TPM1, NF2, PARD3) in KPI from the perspective of AS regulated by RBPs.

Genome-wide identification of alternative splicing and splicing regulated in immune infiltration in osteosarcoma patients

Background: Osteosarcoma typically occurs in adolescents, and the survival rate of patients with metastatic and recurrent osteosarcoma remains low. Abnormal regulation of alternative splicing is associated with the development of osteosarcoma. However, there is no genome-wide analysis of the function and regulatory mechanisms of aberrant alternative splicing associated with osteosarcoma. Methods: Published transcriptome data on osteosarcoma (GSE126209) derived from osteosarcoma patient tissue were downloaded. Gene expression profiling by high-throughput sequencing was performed on 9 normal samples and 10 tumor samples for genome-wide identification of osteosarcoma-related alternative splicing events. The potential function of osteosarcoma-associated alternative splicing events was examined by immune infiltration and correlation analysis. Regulation of aberrantly expressed RNA-binding proteins (RBPs) related to alternative splicing in osteosarcoma was clarified by co-expression analysis. Results: A total of 63 alternative splicing events, which are highly credible and dominant, were identified. GO enrichment analysis indicated that alternative splicing may be closely related to the immune response process. Immune infiltration analysis showed significant changes in the percentages of CD8 T cells, resting memory CD4 T cells, activated memory CD4 T cells, monocytes, resting dendritic cells, and activated mast cells in tumors compared to normal tissues, indicating the involvement of these immune cell types in the occurrence of osteosarcoma. Moreover, the analysis identified alternative splicing events that were co-altered with resting memory CD4 T cells, resting dendritic cells, and activated mast cells, events that may be associated with regulation of the osteosarcoma immune microenvironment. In addition, a co-regulatory network (RBP-RAS-immune) of osteosarcoma-associated RBPs with aberrant alternative splicing and altered immune cells was established. These RBPs include NOP58, FAM120C, DYNC1H1, TRAP1, and LMNA, which may serve as molecular targets for osteosarcoma immune regulation. Conclusion: These findings allow us to further understand the causes of osteosarcoma development and provide a new research direction for osteosarcoma immunotherapy or targeted therapy.

LncRNA TP53TG1 plays an anti-oncogenic role in cervical cancer by synthetically regulating transcriptome profile in HeLa cells

Long non-coding RNAs (lncRNAs) have been extensively studied as important regulators of tumor development in various cancers. Tumor protein 53 target gene 1 (TP53TG1) is a newly identified lncRNA in recent years, and several studies have shown that TP53TG1 may play oncogenic or anti-oncogenic roles in different cancers. Nevertheless, the role of TP53TG1 in the development of cervical cancer is unclear. In our study, pan-cancer analysis showed that high expression of TP53TG1 was significantly associated with a better prognosis. We then constructed a TP53TG1 overexpression model in HeLa cell line to explore its functions and molecular targets. We found that TP53TG1 overexpression significantly inhibited cell proliferation and induced apoptosis, demonstrating that TP53TG1 may be a novel anti-oncogenic factor in cervical cancer. Furthermore, overexpression of TP53TG1 could activate type I interferon signaling pathways and inhibit the expression of genes involved in DNA damage responses. Meanwhile, TP53TG1 could affect alternative splicing of genes involved in cell proliferation or apoptosis by regulating the expression of many RNA-binding protein genes. Competing endogenous RNA (ceRNA) network analysis demonstrated that TP53TG1 could act as the sponge of several miRNAs to regulate the expression level of target genes. In conclusion, our study highlights the essential role of lncRNA TP53TG1 in the development of cervical cancer and suggests the potential regulatory mechanisms.

Genome-Wide Identification of Immune-Related Alternative Splicing and Splicing Regulators Involved in Abdominal Aortic Aneurysm

The molecular mechanism of AAA formation is still poorly understood and has not been fully elucidated. The study was designed to identify the immune-related genes, immune-RAS in AAA using bioinformatics methods. The GSE175683 datasets were downloaded from the GEO database. The DEseq2 software was used to identify differentially expressed genes (DEGs). SUVA pipeline was used to quantify AS events and RAS events. KOBAS 2.0 server was used to identify GO terms and KEGG pathways to sort out functional categories of DEGs. The CIBERSORT algorithm was used with the default parameter for estimating immune cell fractions. Nine samples from GSE175683 were used to construct the co-disturbed network between expression of SFs and splicing ratio of RAS events. PCA analysis was performed by R package factoextra to show the clustering of samples, and the pheatmap package in R was used to perform the clustering based on Euclidean distance. The results showed that there were 3,541 genes significantly differentially expressed, of which 177 immune-related genes were upregulated and 48 immune-related genes were downregulated between the WT and WTA group. Immune-RAS events were mainly alt5P and IR events, and about 60% of it was complex splicing events in AAA. The WT group and the WTA group can be clearly distinguished in the first principal component by using the splicing ratio of immune-RAS events. Two downregulated genes, Nr4a1 and Nr4a2, and eight upregulated genes, Adipor2, Akt2, Bcl3, Dhx58, Pparg, Ptgds, Sytl1, and Vegfa were identified among the immune-related genes with RAS and DEGs. Eighteen differentially expressed SFs were identified and displayed by heatmap. The proportion of different types of cells and ratio of the average ratio of different cells were quite different. Both M1 and M2 types of macrophages and plasma cells were upregulated, while M0 type was downregulated in AAA. The proportion of plasma cells in the WTA group had sharply increased. There is a correlation between SF expression and immune cells/immune-RAS. Sf3b1, a splicing factor with significantly different expression, was selected to bind on a mass of immune-related genes. In conclusion, our results showed that immune-related genes, immune-RAS, and SFs by genome-wide identification were involved in AAA.