Aberrant Splicing Research Articles

Use of RNA‑sequencing to detect abnormal transcription of thecollagen α‑2 (VI) chain gene that can lead to Bethlem myopathy.

Bethlem myopathy (BM) is an autosomal dominant or autosomal recessive disorder and is usually associated with mutations in the collagen VI genes. In the present study, the pathogenicity of a novel splice-site mutation was explored using RNA-sequencing in a family with suspected BM, and a myopathy panel was performed in the proband. The genetic status of all family members was confirmed using Sanger sequencing. Clinical data and magnetic resonance imaging (MRI) features were also documented. In silico analysis was performed to predict the effects of the splice mutation. RNA-sequencing and reverse transcription (RT)-PCR were used to assess aberrant splicing. Immunocytochemistry was conducted to measure collagen VI protein levels within the gastrocnemius and in cultured skin fibroblasts. The results revealed that three patients in the family shared a similar classic BM presentation. MRI revealed distinct patterns of fatty infiltration in the lower extremities. A novel splicing mutation c.736-1G>C in the collagen α-2 (VI) chain (COL6A2) gene was found in all three patients. In silico analysis predicted that the mutation would destroy the normal splice acceptor site. RNA-sequencing detected two abnormal splicing variants adjacent to the mutation site, and RT-PCR confirmed the RNA-sequencing findings. Furthermore, a defect in the collagen protein within cultured fibroblasts was detected using immunocytochemistry. The mutation c.736-1G>C in the COL6A2 gene caused aberrant splicing and led to premature termination of protein translation. In conclusion, these findings may improve our knowledge of mutations of the COL6A2 gene associated with BM and demonstrated that RNA-sequencing can be a powerful tool for finding the underlying mechanism of a disease-causing mutations at a splice site.

A biallelic variant in CLRN2 causes non-syndromic hearing loss in humans

Deafness, the most frequent sensory deficit in humans, is extremely heterogeneous with hundreds of genes involved. Clinical and genetic analyses of an extended consanguineous family with pre-lingual, moderate-to-profound autosomal recessive sensorineural hearing loss, allowed us to identify CLRN2, encoding a tetraspan protein, as a new deafness gene. Homozygosity mapping followed by exome sequencing identified a 14.96 Mb locus on chromosome 4p15.32p15.1 containing a likely pathogenic missense variant in CLRN2 (c.494C > A, NM_001079827.2) segregating with the disease. Using in vitro RNA splicing analysis, we show that the CLRN2 c.494C > A variant leads to two events: (1) the substitution of a highly conserved threonine (uncharged amino acid) to lysine (charged amino acid) at position 165, p.(Thr165Lys), and (2) aberrant splicing, with the retention of intron 2 resulting in a stop codon after 26 additional amino acids, p.(Gly146Lysfs*26). Expression studies and phenotyping of newly produced zebrafish and mouse models deficient for clarin 2 further confirm that clarin 2, expressed in the inner ear hair cells, is essential for normal organization and maintenance of the auditory hair bundles, and for hearing function. Together, our findings identify CLRN2 as a new deafness gene, which will impact future diagnosis and treatment for deaf patients.

Aberrant Splicing Events and Epigenetics in Viral Oncogenomics: Current Therapeutic Strategies.

Global cancer incidence and mortality are on the rise. Although cancer is fundamentally a non-communicable disease, a large number of cancers are known to have a viral aetiology. A high burden of infectious agents (Human immunodeficiency virus (HIV), human papillomavirus (HPV), hepatitis B virus (HBV)) in certain Sub-Saharan African countries drives the rates of certain cancers. About one-third of all cancers in Africa are attributed to infection. Seven viruses have been identified with carcinogenic characteristics, namely the HPV, HBV, Hepatitis C virus (HCV), Epstein–Barr virus (EBV), Human T cell leukaemia virus 1 (HTLV-1), Kaposi’s Sarcoma Herpesvirus (KSHV), and HIV-1. The cellular splicing machinery is compromised upon infection, and the virus generates splicing variants that promote cell proliferation, suppress signalling pathways, inhibition of tumour suppressors, alter gene expression through epigenetic modification, and mechanisms to evade an immune response, promoting carcinogenesis. A number of these splice variants are specific to virally-induced cancers. Elucidating mechanisms underlying how the virus utilises these splice variants to maintain its latent and lytic phase will provide insights into novel targets for drug discovery. This review will focus on the splicing genomics, epigenetic modifications induced by and current therapeutic strategies against HPV, HBV, HCV, EBV, HTLV-1, KSHV and HIV-1.

Detection of aberrant splicing events in RNA-seq data using FRASER

Aberrant splicing is a major cause of rare diseases. However, its prediction from genome sequence alone remains in most cases inconclusive. Recently, RNA sequencing has proven to be an effective complementary avenue to detect aberrant splicing. Here, we develop FRASER, an algorithm to detect aberrant splicing from RNA sequencing data. Unlike existing methods, FRASER captures not only alternative splicing but also intron retention events. This typically doubles the number of detected aberrant events and identified a pathogenic intron retention in MCOLN1 causing mucolipidosis. FRASER automatically controls for latent confounders, which are widespread and affect sensitivity substantially. Moreover, FRASER is based on a count distribution and multiple testing correction, thus reducing the number of calls by two orders of magnitude over commonly applied z score cutoffs, with a minor loss of sensitivity. Applying FRASER to rare disease diagnostics is demonstrated by reprioritizing a pathogenic aberrant exon truncation in TAZ from a published dataset. FRASER is easy to use and freely available.

The RNA-binding protein SFPQ preserves long-intron splicing and regulates circRNA biogenesis in mammals.

Circular RNAs (circRNAs) represent an abundant and conserved entity of non-coding RNAs; however, the principles of biogenesis are currently not fully understood. Here, we identify two factors, splicing factor proline/glutamine rich (SFPQ) and non-POU domain-containing octamer-binding protein (NONO), to be enriched around circRNA loci. We observe a subclass of circRNAs, coined DALI circRNAs, with distal inverted Alu elements and long flanking introns to be highly deregulated upon SFPQ knockdown. Moreover, SFPQ depletion leads to increased intron retention with concomitant induction of cryptic splicing, premature transcription termination, and polyadenylation, particularly prevalent for long introns. Aberrant splicing in the upstream and downstream regions of circRNA producing exons are critical for shaping the circRNAome, and specifically, we identify missplicing in the immediate upstream region to be a conserved driver of circRNA biogenesis. Collectively, our data show that SFPQ plays an important role in maintaining intron integrity by ensuring accurate splicing of long introns, and disclose novel features governing Alu-independent circRNA production.

Therapeutic approaches targeting splicing factor mutations in myelodysplastic syndromes and acute myeloid leukemia.

Mutations in components of the spliceosome are the most common acquired lesions in myelodysplastic syndromes (MDS) and are frequently identified in other myeloid malignancies with a high rate of progression to acute myeloid leukemia (AML) including chronic myelomonocytic leukemia and primary myelofibrosis. The only curative option for these disorders remains allogeneic stem-cell transplantation, which is associated with high morbidity and mortality in these patients. The purpose of this review is to highlight the recent therapeutic developments and strategies being pursued for clinical benefit in splicing factor mutant myeloid malignancies. Cells harboring splicing factor mutations have increased aberrant splicing leading to R-loop formation and cell cycle stalling that create dependencies on Checkpoint kinase 1 (CHK1) activation and canonical splicing maintained by protein arginine methyltransferase activity. Both targeting of the spliceosome and targeting of the downstream consequences of splicing factor mutation expression show promise as selective strategies for the treatment of splicing factor-mutant myeloid malignancies. An improved understanding of the therapeutic vulnerabilities in splicing factor-mutant MDS and AML has led to the development of clinical trials of small molecule inhibitors that target the spliceosome, ataxia telangectasia and Rad3 related (ATR)-CHK1 pathway, and methylation of splicing components.

Interpretation of BRCA2 Splicing Variants: A Case Series of Challenging Variant Interpretations and the Importance of Functional RNA Analysis

A substantial proportion of pathogenic variants associated with an increased risk of hereditary cancer are sequence variants affecting RNA splicing. The classification of these variants can be complex when both non-functional and functional transcripts are produced from the variant allele. We present four BRCA2 splice site variants with complex variant interpretations (BRCA2 c.68-3T>G, c.68-2A>G, c.425G>T, c.8331+2T>C). Evidence supporting a pathogenic classification is available for each variant, including in silico models, absence in population databases, and published functional data. However, comprehensive RNA analysis showed that some functional transcript may be produced by each variant. BRCA2 c.68-3T>G results in a partial splice defect. For BRCA2 c.68-2A>G and c.425G>T, aberrant splicing was shown to produce a potentially functional, in-frame transcript. BRCA2 c.8331+2T>C may utilize a functional GC donor in place of the wild-type GT donor. The severity of cancer history for carriers of these variants was also assessed using a history weighting algorithm and was not consistent with pathogenic controls (carriers of known pathogenic variants in BRCA2). Due to the conflicting evidence, our laboratory classifies these BRCA2 variants as variants of uncertain significance. This highlights the importance of evaluating new and existing evidence to ensure accurate variant classification and appropriate patient care.

Splicing Genomics Events in Cervical Cancer: Insights for Phenotypic Stratification and Biomarker Potency.

Gynaecological cancers are attributed to the second most diagnosed cancers in women after breast cancer. On a global scale, cervical cancer is the fourth most common cancer and the most common cancer in developing countries with rapidly increasing mortality rates. Human papillomavirus (HPV) infection is a major contributor to the disease. HPV infections cause prominent cellular changes including alternative splicing to drive malignant transformation. A fundamental characteristic attributed to cancer is the dysregulation of cellular transcription. Alternative splicing is regulated by several splicing factors and molecular changes in these factors lead to cancer mechanisms such as tumour development and progression and drug resistance. The serine/arginine-rich (SR) proteins and heterogeneous ribonucleoproteins (hnRNPs) have prominent roles in modulating alternative splicing. Evidence shows molecular alteration and expression levels in these splicing factors in cervical cancer. Furthermore, aberrant splicing events in cancer-related genes lead to chemo- and radioresistance. Identifying clinically relevant modifications in alternative splicing events and splicing variants, in cervical cancer, as potential biomarkers for their role in cancer progression and therapy resistance is scrutinised. This review will focus on the molecular mechanisms underlying the aberrant splicing events in cervical cancer that may serve as potential biomarkers for diagnosis, prognosis, and novel drug targets.

Detection of aberrant gene expression events in RNA sequencing data

RNA sequencing (RNA-seq) has emerged as a powerful approach to discover disease-causing gene regulatory defects in individuals affected by genetically undiagnosed rare disorders. Pioneering studies have shown that RNA-seq could increase the diagnosis rates over DNA sequencing alone by 8-36%, depending on the disease entity and tissue probed. To accelerate adoption of RNA-seq by human genetics centers, detailed analysis protocols are now needed. We present a step-by-step protocol that details how to robustly detect aberrant expression levels, aberrant splicing and mono-allelic expression in RNA-seq data using dedicated statistical methods. We describe how to generate and assess quality control plots and interpret the analysis results. The protocol is based on the detection of RNA outliers pipeline (DROP), a modular computational workflow that integrates all the analysis steps, can leverage parallel computing infrastructures and generates browsable web page reports.

X-Linked Retinitis Pigmentosa Caused by Non-Canonical Splice Site Variants in RPGR.

We aimed to validate the effect of non-canonical splice site variants in the RPGR gene in five patients from four families diagnosed with retinitis pigmentosa. Four variants located in intron 2 (c.154 + 3_154 + 6del), intron 3 (c.247 + 5G>A), intron 7 (c.779-5T>G), and intron 13 (c.1573-12A>G), respectively, were analyzed by means of in vitro splice assays. Splicing analysis revealed different aberrant splicing events, including exon skipping and intronic nucleotide addition, which are predicted to lead either to an in-frame deletion affecting relevant protein domains or to a frameshift of the open reading frame. Our data expand the landscape of pathogenic variants in RPGR, thereby increasing the genetic diagnostic rate in retinitis pigmentosa and allowing patients harboring the analyzed variants to be enrolled in clinical trials.

Rare missense variant in MSH4 associated with primary gonadal failure in both 46, XX and 46, XY individuals.

Can whole-exome sequencing (WES) reveal a shared pathogenic variant responsible for primary gonadal failure in both male and female patients from a consanguineous family? Patients with primary ovarian insufficiency (POI) and non-obstructive azoospermia (NOA) were homozygous for the rare missense variant p. S754L located in the highly conserved MSH4 MutS signature motif of the ATPase domain. An oligozoospermic patient was heterozygous for the variant. MSH4 is a meiosis-specific protein expressed at a certain level in the testes and ovaries. Along with its heterodimer partner MSH5, it is responsible for double-strand Holliday junction recognition and stabilization, to ensure accurate chromosome segregation during meiosis. Knockout male and female mice for Msh4 and Msh5 are reportedly infertile due to meiotic arrest. In humans, MSH4 is associated with male and female gonadal failure, with distinct variations in the MutS domain V. This was a retrospective genetics study of a consanguineous family with multiple cases of gonadal failure in both genders. The subject family was recruited in Iran, in 2018. The proband who is affected by POI, an NOA brother, a fertile sister and their parents were subjected to WES. The discovered variant was validated in these individuals, and the rest of the family was also genotyped by Sanger sequencing. The variant was not detected in 800 healthy Iranian individuals from the Iranome database nor in 30 sporadic NOA and 30 sporadic POI patients. Suggested effect in aberrant splicing was studied by RT-PCR. Moreover, protein homology modeling was used to further investigate the amino acid substitution in silico. The discovered variant is very rare and has never been reported in the homozygous state. It occurs in the ATPase domain at Serine 754, the first residue within the highly conserved MutS signature motif, substituting it with a Leucine. All variant effect prediction tools indicated this variant as deleterious. Since the substitution occurs immediately before the Walker B motif at position 755, further investigations based on protein homology were conducted. Considering the modeling results, the nature of the substituted amino acid residue and the distances between p. S754L variation and the residues of the Walker B motif suggested the possibility of conformational changes affecting the ATPase activity of the protein. We have submitted dbSNP entry rs377712900 to ClinVar under SCV001169709, SCV001169708 and SCV001142647 for oligozoospermia, NOA and POI, respectively. Studies in model organisms can shed more light on the role of this variant as our results were obtained by variant effect prediction tools and protein homology modeling. Identification of variants in meiotic genes should improve genetic counseling for both male and female infertility. Also, as two of our NOA patients underwent testicular sperm extraction (TESE) with no success, ruling out the existence of pathogenic variants in meiotic genes in such patients prior to TESE could prove useful. This study was financially supported by Royan Institute in Tehran, Iran, and Institut Pasteur in Paris, France. The authors declare no competing interests. N/A.

Mutations causing Lopes-Maciel-Rodan syndrome are huntingtin hypomorphs.

Huntington's disease pathogenesis involves a genetic gain-of-function toxicity mechanism triggered by the expanded HTT CAG repeat. Current therapeutic efforts aim to suppress expression of total or mutant huntingtin, though the relationship of huntingtin's normal activities to the gain-of-function mechanism and what the effects of huntingtin-lowering might be are unclear. Here, we have re-investigated a rare family segregating two presumed HTT loss-of-function (LoF) variants associated with the developmental disorder, Lopes-Maciel-Rodan syndrome (LOMARS), using whole-genome sequencing of DNA from cell lines, in conjunction with analysis of mRNA and protein expression. Our findings correct the muddled annotation of these HTT variants, reaffirm they are the genetic cause of the LOMARS phenotype and demonstrate that each variant is a huntingtin hypomorphic mutation. The NM_002111.8: c.4469+1G>A splice donor variant results in aberrant (exon 34) splicing and severely reduced mRNA, whereas, surprisingly, the NM_002111.8: c.8157T>A NP_002102.4: Phe2719Leu missense variant results in abnormally rapid turnover of the Leu2719 huntingtin protein. Thus, although rare and subject to an as yet unknown LoF intolerance at the population level, bona fide HTT LoF variants can be transmitted by normal individuals leading to severe consequences in compound heterozygotes due to huntingtin deficiency.

Cytoprotective activities of kinetin purine isosteres

Kinetin (N6-furfuryladenine), a plant growth substance of the cytokinin family, has been shown to modulate aging and various age-related conditions in animal models. Here we report the synthesis of kinetin isosteres with the purine ring replaced by other bicyclic heterocycles, and the biological evaluation of their activity in several in vitro models related to neurodegenerative diseases. Our findings indicate that kinetin isosteres protect Friedreich́s ataxia patient-derived fibroblasts against glutathione depletion, protect neuron-like SH-SY5Y cells from glutamate-induced oxidative damage, and correct aberrant splicing of the ELP1 gene in fibroblasts derived from a familial dysautonomia patient. Although the mechanism of action of kinetin derivatives remains unclear, our data suggest that the cytoprotective activity of some purine isosteres is mediated by their ability to reduce oxidative stress. Further, the studies of permeation across artificial membrane and model gut and blood-brain barriers indicate that the compounds are orally available and can reach central nervous system. Overall, our data demonstrate that isosteric replacement of the kinetin purine scaffold is a fruitful strategy for improving known biological activities of kinetin and discovering novel therapeutic opportunities.

Mutant SF3B1 promotes AKT- and NF-κB-driven mammary tumorigenesis.

Mutations in the core RNA splicing factor SF3B1 are prevalent in leukemias and uveal melanoma, but hotspot SF3B1 mutations are also seen in epithelial malignancies such as breast cancer. Although hotspot mutations in SF3B1 alter hematopoietic differentiation, whether SF3B1 mutations contribute to epithelial cancer development and progression is unknown. Here, we identify that SF3B1 mutations in mammary epithelial and breast cancer cells induce a recurrent pattern of aberrant splicing leading to activation of AKT and NF-κB, enhanced cell migration, and accelerated tumorigenesis. Transcriptomic analysis of human cancer specimens, MMTV-cre Sf3b1K700E/WT mice, and isogenic mutant cell lines identified hundreds of aberrant 3' splice sites (3'ss) induced by mutant SF3B1. Consistently between mouse and human tumors, mutant SF3B1 promoted aberrant splicing (dependent on aberrant branchpoints as well as pyrimidines downstream of the cryptic 3'ss) and consequent suppression of PPP2R5A and MAP3K7, critical negative regulators of AKT and NF-κB. Coordinate activation of NF-κB and AKT signaling was observed in the knockin models, leading to accelerated cell migration and tumor development in combination with mutant PIK3CA but also hypersensitizing cells to AKT kinase inhibitors. These data identify hotspot mutations in SF3B1 as an important contributor to breast tumorigenesis and reveal unique vulnerabilities in cancers harboring them.

Aberrant splicing isoforms detected by full-length transcriptome sequencing as transcripts of potential neoantigens in non-small cell lung cancer

BackgroundLong-read sequencing of full-length cDNAs enables the detection of structures of aberrant splicing isoforms in cancer cells. These isoforms are occasionally translated, presented by HLA molecules, and recognized as neoantigens. This study used a long-read sequencer (MinION) to construct a comprehensive catalog of aberrant splicing isoforms in non-small-cell lung cancers, by which novel isoforms and potential neoantigens are identified.ResultsFull-length cDNA sequencing is performed using 22 cell lines, and a total of 2021 novel splicing isoforms are identified. The protein expression of some of these isoforms is then validated by proteome analysis. Ablations of a nonsense-mediated mRNA decay (NMD) factor, UPF1, and a splicing factor, SF3B1, are found to increase the proportion of aberrant transcripts. NetMHC evaluation of the binding affinities to each type of HLA molecule reveals that some of the isoforms potentially generate neoantigen candidates. We also identify aberrant splicing isoforms in seven non-small-cell lung cancer specimens. An enzyme-linked immune absorbent spot assay indicates that approximately half the peptide candidates have the potential to activate T cell responses through their interaction with HLA molecules. Finally, we estimate the number of isoforms in The Cancer Genome Atlas (TCGA) datasets by referring to the constructed catalog and found that disruption of NMD factors is significantly correlated with the number of splicing isoforms found in the TCGA-Lung Adenocarcinoma data collection.ConclusionsOur results indicate that long-read sequencing of full-length cDNAs is essential for the precise identification of aberrant transcript structures in cancer cells.

RNA-binding motif protein 39 (RBM39): An emerging cancer target.

RNA-binding motif protein 39 (RBM39) is an RNA-binding protein involved in transcriptional co-regulation and alternative RNA splicing. Recent studies have revealed that RBM39 is the unexpected target of aryl sulphonamides, which act as molecular glues between RBM39 and the DCAF15-associated E3 ubiquitin ligase complex leading to selective degradation of the target. Loss of RBM39 leads to aberrant splicing events and differential gene expression, thereby inhibiting cell cycle progression and causing tumour regression in a number of preclinical models. Many clinical studies have shown that aryl sulphonamides were well tolerated, but their clinical performance was limited due to an insufficient understanding of the target, RBM39 biology and a lack of predictive biomarkers. This review summarises the current knowledge of RBM39 function and discusses the therapeutic potential of this spliceosome target in cancer therapy. LINKED ARTICLES: This article is part of a themed issue on New avenues in cancer prevention and treatment (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.12/issuetoc.

Novel insights into clear cell renal cell carcinoma prognosis by comprehensive characterization of aberrant alternative splicing signature: a study based on large-scale sequencing data

ABSTRACT Clear cell renal cell carcinoma (ccRCC) is the most common type with poor prognosis in kidney tumor. Growing evidence has indicated that aberrant alternative splicing (AS) events are efficacious signatures for tumor prognosis prediction and therapeutic targets. However, the detailed roles of AS events in ccRCC are largely unknown. In our study, level 3 RNA-seq data was acquired from The Cancer Genome Atlas dataset and corresponding AS profiles were detected with the assistance of SpliceSeq software. A total of 2100 aberrant survival-associated AS events were identified via differential expression and univariate cox regression analysis. The final prognostic panel formed by 17 specific events was developed by stepwise least absolute shrinkage and selection operator (LASSO) penalty, with the area under curve (AUC) values of receiver operator characteristic (ROC) curves keeping above 0.7 spanning 1 year to 5 years. And the results from functional enrichment analyses are unanimous that autophagy could be a potential mechanism of splicing regulation in ccRCC. Furthermore, splicing regulatory network was constructed via Spearman correlation between splicing factors and AS events. Finally, unsupervised clustering analysis revealed three clusters with distinct survival patterns, and associated with specific clinicopathological phenotypes. In overall, we developed a robust and individualized predictive model based on large-scale sequencing data. The identified AS events and splicing network may be valuable in deciphering the crucial posttranscriptional mechanisms on tumorigenesis of ccRCC.

The Impact of Variants at Branchpoint Splicing Elements in Cancer Genes

Purpose: Branchpoint elements are required for intron removal, and variations at these elements can result in aberrant splicing. We aimed to assess the value of branchpoint annotations generated from recent large-scale studies to select branchpoint-abrogating variants, using hereditary cancer genes as model.Methods: We identified branchpoint elements in 119 genes associated with hereditary cancer from three genome-wide experimentally-inferred and two predicted branchpoint datasets. We then identified variants from public databases that occur within branchpoint elements. We compared conservation, unique variant observations, and population frequencies at different nucleotides within branchpoint motifs. Finally, selected minigene assays were performed to assess the functional impact of variants in branchpoint elements within mismatch repair genes.Results: There was poor overlap between predicted and experimentally-inferred branchpoints. Our analysis of cancer genes suggested that variants at -2nt, -1nt, and branchpoint positions in experimentally-inferred canonical motifs are more likely to be clinically relevant. Minigene assay data showed the -2nt to be more important to branchpoint motif integrity, but also demonstrated fluidity in branchpoint usage.Conclusion: Data from cancer gene analysis suggests that there will be few high-risk alleles that severely impact function via branchpoint abrogation. Results inform a general scheme to prioritize branchpoint motif variants for further study.

Signature of gene aberrant alternative splicing events in pancreatic adenocarcinoma prognosis.

Alternative splicing (AS), as an effective and universal mechanism of transcriptional regulation, is involved in the development and progression of cancer. Therefore, systematic analysis of alternative splicing in pancreatic adenocarcinoma (PAAD) is warranted. The corresponding clinical information of the RNA-Seq data and PAAD cohort was downloaded from the TCGA data portal. Then, a java application, SpliceSeq, was used to evaluate the RNA splicing pattern and calculate the splicing percentage index (PSI). Differentially expressed AS events (DEAS) were identified based on PSI values between PAAD cancer samples and normal samples of adjacent tissues. Kaplan-Meier and Cox regression analyses were used to assess the association between DEAS and patient clinical characteristics. Unsupervised cluster analysis used to reveal four clusters with different survival patterns. At the same time, GEO and TCGA combined with GTEx to verify the differential expression of AS gene and splicing factor. After rigorous filtering, a total of 45,313 AS events were identified, 1,546 of which were differentially expressed AS events. Nineteen DEAS were found to be associated with OS with a five-year overall survival rate of 0.946. And the subtype clusters results indicate that there are differences in the nature of individual AS that affect clinical outcomes. Results also identified 15 splicing factors associated with the prognosis of PAAD. And the splicing factors ESRP1 and RBM5 played an important role in the PAAD-associated AS events. The PAAD-associated AS events, splicing networks, and clusters identified in this study are valuable for deciphering the underlying mechanisms of AS in PAAD and may facilitate the establishment of therapeutic goals for further validation.

ATP5A1 Participates in Transcriptional and Posttranscriptional Regulation of Cancer-Associated Genes by Modulating Their Expression and Alternative Splicing Profiles in HeLa Cells.

Background: Aberrant expression and alternative splicing of oncogenes are the driving events in tumor initiation and development. But how these events are regulated in cancer cells is largely unknown. Functions of ATP5A1, an important mitochondrial ATP synthase gene, in transcriptional and posttranscriptional regulation were explored in this study. Methods: ATP5A1 was overexpressed using plasmid-transformed HeLa cells, and its influence on cell apoptosis and proliferation is evaluated. Transcriptome sequencing was then performed using RNA-seq to study the changes in gene expression and regulation of alternative splicing events. Validation of the implicated genes was achieved using RT-qPCR analysis. Results: It was found that ATP5A1 could significantly promote cellular apoptosis, but it had no influence on cell proliferation. ATP5A1 overexpression significantly increased the expression levels of genes associated with the innate immune response, angiogenesis, and collagen catabolic processes. This included enrichment of MMP2 and MMP19. It was also found that ATP5A1 could interfere with the alternative splicing of hundreds of genes associated with glucose homeostasis, HIF-1 signaling activation, and several pathways associated with cancers. Eight ATP5A1-regulated differentially expressed genes and 3 genes altered by splicing were selected and validated using RT-qPCR analysis. Conclusions: In summary, we illustrate the regulatory functions of ATP5A1 on the transcriptome of HeLa cells by exploring its influence on gene expression and alternative splicing. The results suggest that ATP5A1 may play an important regulatory role in cervical cancer cells by regulating expression and alternative splicing of cancer-associated genes. This study provides novel insights into the current understanding of the mechanisms of ATP5A1 on carcinogenesis and cancer progression.