Intron retention is a widespread mechanism of tumor-suppressor inactivation

An IGF-1R-mTORC1-SRPK2 signaling Axis contributes to FASN regulation in breast cancer

Background

Fatty acid synthase (FASN) expression is associated with a more aggressive breast cancer phenotype and is regulated downstream of receptor tyrosine kinase (RTK) signaling pathways. Recently, post transcriptional regulation of lipogenic transcripts have been demonstrated as being mediated downstream of serine-arginine rich protein kinase 2 (SRPK2), which acts to phosphorylate serine-arginine rich splicing factors (SRSFs), resulting in RNA binding and various RNA regulatory processes. Though post-transcriptional regulation of FASN has been studied previously, the upstream mediators of these pathways have not been elucidated.

Methods

Western blotting and RT-qPCR were utilized to demonstrate alterations in FASN and mRNA expression upon modulation of the IGF-1-mTORC1-SRPK2 pathway by small molecule inhibitors or RNAi mediated silencing. RNA stability was accessed by using the transcriptional inhibitor actinomycin-D followed by RT-qPCR. Further, we employed RNA-immunoprecipitation to demonstrate the direct binding of SRSF-1 to FASN transcripts.

Results

In the current study, we demonstrated an IGF-1 induced increase in FASN mRNA and protein expression that was attenuated by mTORC1 inhibition. This mTORC1 inhibition also resulted in decreases in total and nuclear p-SRPK2 in response to IGF-1 exposure. Upon SRPK2 knockdown and inhibition, we observed a decrease in FASN protein and mRNA stability, respectively, in response to IGF-1 exposure that was specific to triple negative and HER2+ breast cancer cell lines. As we explored further, IGF-1 exposure resulted in an altered localization of eGFP expressed SRSF-1, pEGFP-SRSF-1 that was rescued upon both SRPK2 knockdown and mTORC1 inhibition. Further, we observed an increase binding of SRSF-1 to FASN RNA upon IGF-1 exposure, which was abrogated by SRPK2 knockdown.

Conclusion

These current findings establish a potential IGF-1-mTORC1-SRPK2-FASN axis in breast cancer, which could be a potential therapeutic target for cancers that overexpress FASN and components of the IGF-1R pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10062-z.

Fyn and Lyn gene polymorphisms impact the risk of thyroid cancer

Thyroid cancer is the most common malignancy of the endocrine glands, and during last couple of decades, its incidence has risen alarmingly, across the globe. Etiology of thyroid cancer is still debatable. There are a few worth mentioning risk factors which contribute to initiation of abnormalities in thyroid gland leading to cancer. Genetic instability is major risk factors in thyroid carcinogenesis. Among the genetic factors, the Src family of genes (Src, Yes1, Fyn and Lyn) have been implicated in many cancers but there is little data regarding the association of these (Src, Yes1, Fyn and Lyn) genes with thyroid carcinogenesis. Fyn and Lyn genes of Src family found engaged in proliferation, migration, invasion, angiogenesis, and metastasis in different cancers. This study was planned to examine the effect of Fyn and Lyn SNPs on thyroid cancer risk in Pakistani population in 500 patients and 500 controls. Three polymorphisms of Fyn gene (rs6916861, rs2182644 and rs12910) and three polymorphisms of Lyn gene (rs2668011, rs45587541 and rs45489500) were analyzed using Tetra-primer ARMS-PCR followed by DNA sequencing. SNP rs6916861 of Fyn gene mutant genotype (CC) showed statistically significant threefold increased risk of thyroid cancer (P < 0.0001). In case of rs2182644 of Fyn gene, mutant genotype (AA) indicated statistically significant 17-fold increased risk of thyroid cancer (P < 0.0001). Statistically significant threefold increased risk of thyroid cancer was observed in genotype AC (P < 0.0001) of Fyn gene polymorphism rs12910. In SNP rs2668011 of Lyn gene, TT genotype showed statistically significant threefold increased risk of thyroid cancer (P < 0.0001). In case of rs45587541 of Lyn gene, GA genotypes showed statistically significant 11-fold increased risk in thyroid cancer (P < 0.0001). Haplotype analysis revealed that AAATAG*, AGACAG*, AGCCAA*, AGCCAG*, CAATAG*, CGCCAG* and CGCCGA* haplotypes of Fyn and Lyn polymorphisms are associated with increased thyroid cancer risk. These results showed that genotypes and allele distribution of Fyn and Lyn are significantly linked with increased thyroid cancer risk and could be genetic adjuster for said disease.

Nonsense-mediated RNA decay: an emerging modulator of malignancy

Nonsense-mediated RNA decay (NMD) is a highly conserved RNA turnover pathway that selectively degrades RNAs harbouring truncating mutations that prematurely terminate translation, including nonsense, frameshift and some splice-site mutations. Recent studies show that NMD shapes the mutational landscape of tumours by selecting for mutations that tend to downregulate the expression of tumour suppressor genes but not oncogenes. This suggests that NMD can benefit tumours, a notion further supported by the finding that mRNAs encoding immunogenic neoantigen peptides are typically targeted for decay by NMD. Together, this raises the possibility that NMD-inhibitory therapy could be of therapeutic benefit against many tumour types, including those with a high load of neoantigen-generating mutations. Complicating this scenario is the evidence that NMD can also be detrimental for many tumour types, and consequently tumours often have perturbed NMD. NMD may suppress tumour generation and progression by degrading subsets of specific normal mRNAs, including those encoding stress-response proteins, signalling factors and other proteins beneficial for tumours, as well as pro-tumour non-coding RNAs. Together, these findings suggest that NMD-modulatory therapy has the potential to provide widespread therapeutic benefit against diverse tumour types. However, whether NMD should be stimulated or repressed requires careful analysis of the tumour to be treated.

The splicing factor SF3B4 drives proliferation and invasion in cervical cancer by regulating SPAG5

Regulation of alternative splicing (AS) by the splicing factor 3b (SF3B) family plays an essential role in cancer. However, the biological function of SF3B family members in cervical cancer (CC) needs to be further elucidated. In this study, we found that splicing factor 3b subunit 4 (SF3B4) was highly expressed in CC by bioinformatics analysis using cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) data from The Cancer Genome Atlas (TCGA). Then, we demonstrated that high expression of SF3B4 promoted proliferation and invasion abilities of CC cells in vitro and in vivo and that reduced expression of SF3B4 performed the opposite effect. Further RNA-seq and AS analysis showed that sperm-associated antigen 5 (SPAG5) was a downstream target gene of SF3B4. Interestingly, SPAG5 expression was decreased after SF3B4 knockdown because of retained introns (RIs) and reduced maturation of SPAG5 pre-mRNA. Importantly, SPAG5 deficiency impaired the oncogenic effects of SF3B4 overexpression on CC cells. In conclusion, SF3B4 promotes CC progression by regulating the effective splicing of SPAG5. SF3B4 could be a promising target for CC.

The broader sense of nonsense

The term 'nonsense-mediated mRNA decay' (NMD) was initially coined to describe the translation-dependent degradation of mRNAs harboring premature termination codons (PTCs), but it is meanwhile known that NMD also targets many canonical mRNAs with numerous biological implications. The molecular mechanisms determining on which RNAs NMD ensues are only partially understood. Considering the broad range of NMD-sensitive RNAs and the variable degrees of their degradation, we highlight here the hallmarks of mammalian NMD and point out open questions. We review the links between NMD and disease by summarizing the role of NMD in cancer, neurodegeneration, and viral infections. Finally, we describe strategies to modulate NMD activity and specificity as potential therapeutic approaches for various diseases.

Dysregulation and therapeutic targeting of RNA splicing in cancer

High-throughput sequencing and functional characterization of the cancer transcriptome have uncovered cancer-specific dysregulation of RNA splicing across a variety of cancers. Alterations in the cancer genome and dysregulation of RNA splicing factors lead to missplicing, splicing alteration-dependent gene expression and, in some cases, generation of novel splicing-derived proteins. Here, we review recent advances in our understanding of aberrant splicing in cancer pathogenesis and present strategies to harness cancer-specific aberrant splicing for therapeutic intent.

Therapeutic Targeting of Alternative Splicing: A New Frontier in Cancer Treatment

The ability for cells to harness alternative splicing enables them to diversify their proteome in order to carry out complex biological functions and adapt to external and internal stimuli. The spliceosome is the multiprotein-RNA complex charged with the intricate task of alternative splicing. Aberrant splicing can arise from abnormal spliceosomes or splicing factors and drive cancer development and progression. This review will provide an overview of the alternative splicing process and aberrant splicing in cancer, with a focus on serine/arginine-rich (SR) proteins and their recently reported roles in cancer development and progression and beyond. Recent mapping of the spliceosome, its associated splicing factors, and their relationship to cancer have opened the door to novel therapeutic approaches that capitalize on the widespread influence of alternative splicing. We conclude by discussing small molecule inhibitors of the spliceosome that have been identified in an evolving era of cancer treatment.

miRNA-Mediated Low Expression of EPHX3 Is Associated with Poor Prognosis and Tumor Immune Infiltration in Head and Neck Squamous Cell Carcinomas

The aim of this study was to explore the regulatory role of epoxide hydrolase 3 (EPHX3) in head and neck squamous cell carcinoma (HNSCC) and to analyze its bioinformatic function, as well as, to screen and predict the miRNAs that can regulate EPHX3 expression in HNSCC. We examined the expression profile and prognostic potential of EPHX3 in TCGA and GTEX databases and performed functional enrichment analysis of EPHX3 using string database. Subsequently, we analyzed the regulatory role of miRNAs on EPHX3, including expression analysis, correlation analysis, and survival analysis. In addition, we also used TIMER to investigate the relationship among EPHX3 expression level, immune checkpoints, and immune infiltration in HNSCC. The results of data analysis after TGCA showed that EPHX3 is a key regulator of tumorigenesis in 13 cancers and can be used as a marker of poor prognosis in HNSCC patients. Bioinformatics analysis revealed that miR-4713-3p is a key miRNA of EPHX3 in HNSCC. Together, our findings indicate that EPHX3 exerts its anticancer effects by suppressing tumor immune checkpoint expression and immune cell infiltration. Overall, our data uncovered miRNA-mediated EPHX3 downregulation as a contributor to poor HNSCC prognosis and reduced tumor immune infiltration.

CRLF2 overexpression defines an immature-like subgroup which is rescued through restoration of the PRC2 function in T-cell precursor acute lymphoblastic leukaemia

CRLF2 overexpression has been described as a biomarker of poor prognosis in T-cell acute lymphoblastic leukaemia (T-ALL). In the present study, we aimed to unravel the genomic profile underlying CRLF2 overexpression (CRLF2-high) by analysing RNA-seq, WES and SNP-array data from 264 T-ALL patients and five cell lines deposited on the TARGET initiative, Cancer Cell Line Encyclopedia and Gene Expression Omnibus. These data allowed us to delineate the genomic landscape of CRLF2-high in T-ALL, which was associated with PTEN, JAK3, PHF6, EZH2 and RUNX1 mutations. We also observed an enrichment of CRLF2-high in early T-precursor (ETP)-ALL (23.08% vs 4.02%, P = 7.579e^-06 ) and a very similar gene upregulation profile between these two entities. The inhibition of BET (iBET) proteins is a strategy previously demonstrated to reverse the gene upregulation pattern of ETP cells through restoration of polycomb repressive complex 2 (PRC2) activity. While CRLF2 expression was rescued by using this strategy in LOUCY (untreated vs iBET P = 0.0095, DMSO vs iBET P = 0.0286), a classical ETP-ALL cell line, PRC2 loss was not sufficient to promote CRLF2 upregulation in JURKAT, a more mature T-ALL cell line. Considering the role of IKZF1 in CRLF2 regulation and in recruitment of PCR2, we evaluated IKZF1 status according to CRLF2-expression subgroups. We identified that IKZF1 transcripts with intron retention were upregulated in the CRLF2-high subgroup. Here, we delineated the gene expression profile of CRLF2-high T-ALL samples and unravelled the crucial role of PRC2 in CRLF2 regulation in ETP-ALL.

Impacts and mechanisms of alternative mRNA splicing in cancer metabolism, immune response, and therapeutics

Alternative pre-mRNA splicing (AS) provides the potential to produce diversity at RNA and protein levels. Disruptions in the regulation of pre-mRNA splicing can lead to diseases. With the development of transcriptome and genome sequencing technology, increasing diseases have been identified to be associated with abnormal splicing of mRNAs. In tumors, abnormal alternative splicing frequently plays critical roles in cancer pathogenesis and may be considered as new biomarkers and therapeutic targets for cancer intervention. Metabolic abnormalities and immune disorders are important hallmarks of cancer. AS produces multiple different isoforms and diversifies protein expression, which is utilized by the immune and metabolic reprogramming systems to expand gene functions. The abnormal splicing events contributed to tumor progression, partially due to effects on immune response and metabolic reprogramming. Herein, we reviewed the vital role of alternative splicing in regulating cancer metabolism and immune response. We discussed how alternative splicing regulates metabolic reprogramming of cancer cells and antitumor immune response, and the possible strategies to targeting alternative splicing pathways or splicing-regulated metabolic pathway in the context of anticancer immunotherapy. Further, we highlighted the challenges and discuss the perspectives for RNA-based strategies for the treatment of cancer with abnormally alternative splicing isoforms.

Secondhand smoke affects reproductive functions by altering the mouse testis transcriptome, and leads to select intron retention in Pde1a

BACKGROUND

Human exposure to secondhand smoke (SHS) is known to result in adverse effects in multiple organ systems. However, the impact of SHS on the male reproductive system, particularly on the regulation of genes and molecular pathways that govern sperm production, maturation, and functions remains largely understudied.

OBJECTIVE

We investigated the effects of SHS on the testis transcriptome in a validated mouse model.

METHODS

Adult male mice were exposed to SHS (5 h/day, 5 days/week for 4 months) as compared to controls (clean air-exposed). RNA-seq analysis was performed on the testis of SHS-exposed mice and controls. Variant discovery and plink association analyses were also conducted to detect exposure-related transcript variants in SHS-treated mice.

RESULTS

Exposure of mice to SHS resulted in the aberrant expression of 131 testicular genes. Whilst approximately two thirds of the differentially expressed genes were protein-coding, the remaining (30.5%) comprised noncoding elements, mostly lncRNAs (19.1%). Variant discovery analysis identified a homozygous frameshift variant that is statistically significantly associated with SHS exposure (P = 7.744e-06) and is generated by retention of a short intron within Pde1a, a key regulator of spermatogenesis. Notably, this SHS-associated intron variant harbors an evolutionarily conserved, premature termination codon (PTC) that disrupts the open reading frame of Pde1a, presumably leading to its degradation via nonsense-mediated decay.

DISCUSSION

SHS alters the expression of genes involved in molecular pathways that are crucial for normal testis development and function. Preferential targeting of lncRNAs in the testis of SHS-exposed mice is especially significant considering their crucial role in the spatial and temporal modulation of spermatogenesis. Equally important is our discovery of a novel homozygous frameshift variant that is exclusively and significantly associated with SHS-exposure and is likely to represent a safeguard mechanism to regulate transcription of Pde1a and preserve normal testis function during harmful exposure to environmental agents.

Combination therapy targeting Erk1/2 and CDK4/6i in relapsed refractory multiple myeloma

Oncogenic activated RAS mutations have been detected in 50% of de novo and 70% of relapsed multiple myeloma (MM) patients. Translocation t(11;14) involving IgH/CCDN1 and overexpression of cyclin-Ds are early events in MM pathogenesis, enhancing uncontrolled MM cell growth. We hypothesized that targeting both RAS/MAPK pathway molecules including Erk1/2 along with cyclin-Ds enhances MM cytotoxicity and minimizes side effects. Recent studies have demonstrated the high potency of Erk1/2 and CDK4/6 inhibitors in metastatic relapsed cancers, and here we tested anti-MM effects of the Erk1/2 + CDK4/6 inhibitor combination. Our studies showed strong synergistic (IC < 0.5) cytotoxicity of Erk1/2i + CDK4/6i in MM-cells. Erk1/2i + CDK4/6i treatment in a dose-dependent manner arrested MM-cells in the G0/G1 phase and activated mitochondrial apoptotic signaling. Our studies showed that Erk1/2i + CDK4/6i treatment-induced inhibition of key target molecules in Erk1/2 and CDK4/6 signaling, such as c-myc, p-RSK, p-S6, p-RB, and E2F1, suggesting on-target activity of these inhibitors. We identified Erk1/2i + CDK4/6i treatment associated five-gene signature which includes SNRPB and SLC25A5; these genes are involved in RNA processing and mitochondrial metabolism, respectively. Overall, our studies provide the preclinical framework for Erk1/2i + CDK4/6i combination clinical trials to target Ras+CDK pathways to improve patient outcome in MM.

Alternative Splicing, Epigenetic Modifications and Cancer: A Dangerous Triangle, or a Hopeful One?

Simple Summary

Epigenetics studies the alteration of gene expression without changing DNA sequence and very often, epigenetic dysregulation causes cancer. Alternative splicing is a mechanism that results in the production of several mRNA isoforms from a single gene and aberrant splicing is also a frequent cause of cancer. The present review is built on the interrelations of epigenetics and alternative splicing. In an intuitive way, we say that epigenetic modifications and alternative splicing are at two vertices of a triangle, the third vertex being occupied by cancer. Interconnection between alternative splicing and epigenetic modifications occurs backward and forward and the mechanisms involved are widely reviewed. These connections also provide novel diagnostic or prognostic tools, which are listed. Finally, as epigenetic alterations are reversible and aberrant alternative splicing may be corrected, the therapeutic possibilities to break the triangle are discussed.

Abstract

The alteration of epigenetic modifications often causes cancer onset and development. In a similar way, aberrant alternative splicing may result in oncogenic products. These issues have often been individually reviewed, but there is a growing body of evidence for the interconnection of both causes of cancer. Actually, aberrant splicing may result from abnormal epigenetic signalization and epigenetic factors may be altered by alternative splicing. In this way, the interrelation between epigenetic marks and alternative splicing form the base of a triangle, while cancer may be placed at the vertex. The present review centers on the interconnections at the triangle base, i.e., between alternative splicing and epigenetic modifications, which may result in neoplastic transformations. The effects of different epigenetic factors, including DNA and histone modifications, the binding of non-coding RNAs and the alterations of chromatin organization on alternative splicing resulting in cancer are first considered. Other less-frequently considered questions, such as the epigenetic regulation of the splicing machinery, the aberrant splicing of epigenetic writers, readers and erasers, etc., are next reviewed in their connection with cancer. The knowledge of the above-mentioned relationships has allowed increasing the collection of biomarkers potentially useful as cancer diagnostic and/or prognostic tools. Finally, taking into account on one hand that epigenetic changes are reversible, and some epigenetic drugs already exist and, on the other hand, that drugs intended for reversing aberrations in alternative splicing, therapeutic possibilities for breaking the mentioned cancer-related triangle are discussed.

Characterization of Alternative Splicing Events in Porcine Skeletal Muscles with Different Intramuscular Fat Contents

Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.

TP53_PROF: a machine learning model to predict impact of missense mutations in TP53

Correctly identifying the true driver mutations in a patient’s tumor is a major challenge in precision oncology. Most efforts address frequent mutations, leaving medium- and low-frequency variants mostly unaddressed. For TP53, this identification is crucial for both somatic and germline mutations, with the latter associated with the Li-Fraumeni syndrome (LFS), a multiorgan cancer predisposition. We present TP53_PROF (prediction of functionality), a gene specific machine learning model to predict the functional consequences of every possible missense mutation in TP53, integrating human cell- and yeast-based functional assays scores along with computational scores. Variants were labeled for the training set using well-defined criteria of prevalence in four cancer genomics databases. The model’s predictions provided accuracy of 96.5%. They were validated experimentally, and were compared to population data, LFS datasets, ClinVar annotations and to TCGA survival data. Very high accuracy was shown through all methods of validation. TP53_PROF allows accurate classification of TP53 missense mutations applicable for clinical practice. Our gene specific approach integrated machine learning, highly reliable features and biological knowledge, to create an unprecedented, thoroughly validated and clinically oriented classification model. This approach currently addresses TP53 mutations and will be applied in the future to other important cancer genes.

Differential fates of introns in gene expression due to global alternative splicing

The discovery of introns over four decades ago revealed a new vision of genes and their interrupted arrangement. Throughout the years, it has appeared that introns play essential roles in the regulation of gene expression. Unique processing of excised introns through the formation of lariats suggests a widespread role for these molecules in the structure and function of cells. In addition to rapid destruction, these lariats may linger on in the nucleus or may even be exported to the cytoplasm, where they remain stable circular RNAs (circRNAs). Alternative splicing (AS) is a source of diversity in mature transcripts harboring retained introns (RI-mRNAs). Such RNAs may contain one or more entire retained intron(s) (RIs), but they may also have intron fragments resulting from sequential excision of smaller subfragments via recursive splicing (RS), which is characteristic of long introns. There are many potential fates of RI-mRNAs, including their downregulation via nuclear and cytoplasmic surveillance systems and the generation of new protein isoforms with potentially different functions. Various reports have linked the presence of such unprocessed transcripts in mammals to important roles in normal development and in disease-related conditions. In certain human neurological-neuromuscular disorders, including myotonic dystrophy type 2 (DM2), frontotemporal dementia/amyotrophic lateral sclerosis (FTD/ALS) and Duchenne muscular dystrophy (DMD), peculiar processing of long introns has been identified and is associated with their pathogenic effects. In this review, we discuss different mechanisms involved in the processing of introns during AS and the functions of these large sections of the genome in our biology.

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

Clinical and Functional Significance of TP53 Exon 4-Intron 4 Splice Junction Variants

Germline TP53 splicing variants are uncommon, and their clinical relevance is unknown. However, splice-altering variants at exon 4-intron 4 junctions are relatively enriched in pediatric adrenocortical tumors (ACT). Nevertheless, family histories of cancer compatible with classic Li-Fraumeni syndrome are rarely seen in these patients. We used conventional and in silico assays to determine protein stability, splicing, and transcriptional activity of 10 TP53 variants at exon 4-intron 4 junctions and analyzed their clinical correlates. We reviewed public databases that report the impact of TP53 variants in human cancer and examined individual reports, focusing on family history of cancer. TP53 exon 4-intron 4 junction germline variants were identified in 9 of 75 pediatric ACTs enrolled in the International Pediatric Adrenocortical Tumor Registry and Children's Oncology Group ARAR0332 study. An additional eight independent TP53 variants involving exon 4 splicing were identified in the Pediatric Cancer Genome Project (n = 5,213). These variants resulted in improper expression due to ineffective splicing, protein instability, altered subcellular localization, and loss of function. Clinical case review of carriers of TP53 exon 4-intron 4 junction variants revealed a high incidence of pediatric ACTs and atypical tumor types not consistent with classic Li-Fraumeni syndrome. Germline variants involving TP53 exon 4-intron 4 junctions are frequent in ACT and rare in other pediatric tumors. The collective impact of these germline TP53 variants on the fidelity of splicing, protein structure, and function must be considered in evaluating cancer susceptibility. IMPLICATIONS: Taken together, the data indicate that splice variants at TP53 codon 125 and surrounding bases differentially impacted p53 gene expression and function.

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Systematic analysis of the relationship between ovarian cancer prognosis and alternative splicing

Background

Ovarian cancer(OC) is the gynecological tumor with the highest mortality rate, effective biomarkers are of great significance in improving its prognosis. In recent years, there have been many studies on alternative splicing (AS) events, and the role of AS events in tumor has become a focus of attention.

Methods

Data were downloaded from the TCGA database and Univariate Cox regression analysis was performed to determine AS events associated with OC prognosis.Eight prognostic models of OC were constructed in R package, and the accuracy of the models were evaluated by the time-dependent receiver operating characteristic (ROC) curves.Eight types of survival curves were drawn to evaluate the differences between the high and low risk groups.Independent prognostic factors of OC were analyzed by single factor independent analysis and multi-factor independent prognostic analysis.Again, Univariate Cox regression analysis was used to analyze the relationship between splicing factors(SF) and AS events, and Gene Ontology(GO) and Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment analysis were performed on OS-related SFs to understand the pathways.

Results

Univariate Cox regression analysis showed that among the 15,278 genes, there were 31,286 overall survival (OS) related AS events, among which 1524 AS events were significantly correlated with OS. The area under the time-dependent receiver operating characteristic curve (AUC) of AT and ME were the largest and the RI was the smallest,which were 0.757 and 0.68 respectively. The constructed models have good value for the prognosis assessment of OC patients. Among the eight survival curves, AP was the most significant difference between the high and low risk groups, with a P value of 1.61e − 1.The results of single factor independent analysis and multi-factor independent prognostic analysis showed that risk score calculated by the model and age could be used as independent risk factors.According to univariate COX regression analysis,109 SFs were correlated with AS events and adjusted in two ways: positive and negative.

Conclusions

SFs and AS events can directly or indirectly affect the prognosis of OC patients. It is very important to find effective prognostic markers to improve the survival rate of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-021-00866-1.

Immunopeptidogenomics: Harnessing RNA-Seq to Illuminate the Dark Immunopeptidome

Human leukocyte antigen (HLA) molecules are cell-surface glycoproteins that present peptide antigens on the cell surface for surveillance by T lymphocytes, which contemporaneously seek signs of disease. Mass spectrometric analysis allows us to identify large numbers of these peptides (the immunopeptidome) following affinity purification of solubilized HLA-peptide complexes. However, in recent years, there has been a growing awareness of the "dark side" of the immunopeptidome: unconventional peptide epitopes, including neoepitopes, which elude detection by conventional search methods because their sequences are not present in reference protein databases (DBs). Here, we establish a bioinformatics workflow to aid identification of peptides generated by noncanonical translation of mRNA or by genome variants. The workflow incorporates both standard transcriptomics software and novel computer programs to produce cell line-specific protein DBs based on three-frame translation of the transcriptome. The final protein DB also includes sequences resulting from variants determined by variant calling on the same RNA-Seq data. We then searched our experimental data against both transcriptome-based and standard DBs using PEAKS Studio (Bioinformatics Solutions, Inc). Finally, further novel software helps to compare the various result sets arising for each sample, pinpoint putative genomic origins for unconventional sequences, and highlight potential neoepitopes. We applied the workflow to study the immunopeptidome of the acute myeloid leukemia cell line THP-1, using RNA-Seq and immunopeptidome data. We confidently identified over 14,000 peptides from three replicates of purified HLA peptides derived from THP-1 cells using the conventional UniProt human proteome. Using the transcriptome-based DB generated using our workflow, we recapitulated >85% of these and also identified 1029 unconventional peptides not explained by UniProt, including 16 sequences caused by nonsynonymous variants. Our workflow, which we term "immunopeptidogenomics," can provide DBs, which include pertinent unconventional sequences and allow neoepitope discovery, without becoming too large to search. Immunopeptidogenomics is a step toward unbiased search approaches that are needed to illuminate the dark side of the immunopeptidome.

The expressed mutational landscape of microsatellite stable colorectal cancers

Background

Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression.

Methods

Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids.

Results

The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models.

Conclusions

Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the “expressed mutation dose” may provide an opportunity for more fine-tuned biomarker interpretations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00955-2.

Alternative microexon splicing by RBFOX2 and PTBP1 is associated with metastasis in colorectal cancer

The splicing of microexons (very small exons) is frequently dysregulated in the brain of individuals with autism spectrum disorder. However, little is known of the patterns, regulatory mechanisms and roles of microexon splicing in cancer. We here examined the transcriptome-wide profile of microexon splicing in matched colorectal cancer (CRC) and normal tissue specimens. Out of 1492 microexons comprising 3 to 15 nucleotides, 21 (1%) manifested differential splicing between CRC and normal tissue. The 21 genes harboring the differentially spliced microexons were enriched in gene ontology terms related to cell adhesion and migration. RNA interference-mediated knockdown experiments identified two splicing factors, RBFOX2 and PTBP1, as regulators of microexon splicing in CRC cells. RBFOX2 and PTBP1 were found to directly bind to microexon-containing pre-mRNAs and to control their splicing in such cells. Differential microexon splicing was shown to be due, at least in part, to altered expression of RBFOX2 and PTBP1 in CRC tissue compared to matched normal tissue. Finally, we found that changes in the pattern of microexon splicing were associated with CRC metastasis. Our data thus suggest that altered expression of RBFOX2 and PTBP1 might influence CRC metastasis through the regulation of microexon splicing.

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

Most of the current alternative splicing (AS) analysis tools are powerless to analyse complex splicing. To address this, we developed SUVA (Splice sites Usage Variation Analysis) that decomposes complex splicing events into five types of splice junction pairs. By analysing real and simulated data, SUVA showed higher sensitivity and accuracy in detecting AS events than the compared methods. Notably, SUVA detected extensive complex AS events and screened out 69 highly conserved and dominant AS events associated with cancer. The cancer-associated complex AS events in FN1 and the co-regulated RNA-binding proteins were significantly correlated with patient survival.

HNRNPM controls circRNA biogenesis and splicing fidelity to sustain cancer cell fitness

High spliceosome activity is a dependency for cancer cells, making them more vulnerable to perturbation of the splicing machinery compared to normal cells. To identify splicing factors important for prostate cancer (PCa) fitness, we performed pooled shRNA screens in vitro and in vivo. Our screens identified heterogeneous nuclear ribonucleoprotein M (HNRNPM) as a regulator of PCa cell growth. RNA- and eCLIP-sequencing identified HNRNPM binding to transcripts of key homeostatic genes. HNRNPM binding to its targets prevents aberrant exon inclusion and backsplicing events. In both linear and circular mis-spliced transcripts, HNRNPM preferentially binds to GU-rich elements in long flanking proximal introns. Mimicry of HNRNPM-dependent linear-splicing events using splice-switching-antisense-oligonucleotides was sufficient to inhibit PCa cell growth. This suggests that PCa dependence on HNRNPM is likely a result of mis-splicing of key homeostatic coding and non-coding genes. Our results have further been confirmed in other solid tumors. Taken together, our data reveal a role for HNRNPM in supporting cancer cell fitness. Inhibition of HNRNPM activity is therefore a potential therapeutic strategy in suppressing growth of PCa and other solid tumors.

Integrative functional genomic analysis of intron retention in human and mouse brain with Alzheimer's disease

Intron retention (IR) has been implicated in the pathogenesis of complex diseases such as cancers; its association with Alzheimer's disease (AD) remains unexplored. We performed genome-wide analysis of IR through integrating genetic, transcriptomic, and proteomic data of AD subjects and mouse models from the Accelerating Medicines Partnership-Alzheimer's Disease project. We identified 4535 and 4086 IR events in 2173 human and 1736 mouse genes, respectively. Quantitation of IR enabled the identification of differentially expressed genes that conventional exon-level approaches did not reveal. There were significant correlations of intron expression within innate immune genes, like HMBOX1, with AD in humans. Peptides with a high probability of translation from intron-retained mRNAs were identified using mass spectrometry. Further, we established AD-specific intron expression Quantitative Trait Loci, and identified splicing-related genes that may regulate IR. Our analysis provides a novel resource for the search for new AD biomarkers and pathological mechanisms.

A novel class of inhibitors that target SRSF10 and promote p53-mediated cytotoxicity on human colorectal cancer cells

The elevated expression of the splicing regulator SRSF10 in metastatic colorectal cancer (CRC) stimulates the production of the pro-tumorigenic BCLAF1-L splice variant. We discovered a group of small molecules with an aminothiazole carboxamide core (GPS167, GPS192 and others) that decrease production of BCLAF1-L. While additional alternative splicing events regulated by SRSF10 are affected by GPS167/192 in HCT116 cells (e.g. in MDM4, WTAP, SLK1 and CLK1), other events are shifted in a SRSF10-independent manner (e.g. in MDM2, NAB2 and TRA2A). GPS167/192 increased the interaction of SRSF10 with the CLK1 and CLK4 kinases, leading us to show that GPS167/192 can inhibit CLK kinases preferentially impacting the activity of SRSF10. Notably, GPS167 impairs the growth of CRC cell lines and organoids, inhibits anchorage-independent colony formation, cell migration, and promotes cytoxicity in a manner that requires SRSF10 and p53. In contrast, GPS167 only minimally affects normal colonocytes and normal colorectal organoids. Thus, GPS167 reprograms the tumorigenic activity of SRSF10 in CRC cells to elicit p53-dependent apoptosis.

Biology of the mRNA Splicing Machinery and Its Dysregulation in Cancer Providing Therapeutic Opportunities

Dysregulation of messenger RNA (mRNA) processing—in particular mRNA splicing—is a hallmark of cancer. Compared to normal cells, cancer cells frequently present aberrant mRNA splicing, which promotes cancer progression and treatment resistance. This hallmark provides opportunities for developing new targeted cancer treatments. Splicing of precursor mRNA into mature mRNA is executed by a dynamic complex of proteins and small RNAs called the spliceosome. Spliceosomes are part of the supraspliceosome, a macromolecular structure where all co-transcriptional mRNA processing activities in the cell nucleus are coordinated. Here we review the biology of the mRNA splicing machinery in the context of other mRNA processing activities in the supraspliceosome and present current knowledge of its dysregulation in lung cancer. In addition, we review investigations to discover therapeutic targets in the spliceosome and give an overview of inhibitors and modulators of the mRNA splicing process identified so far. Together, this provides insight into the value of targeting the spliceosome as a possible new treatment for lung cancer.

Tumor reversion: a dream or a reality

Reversion of tumor to a normal differentiated cell once considered a dream is now at the brink of becoming a reality. Different layers of molecules/events such as microRNAs, transcription factors, alternative RNA splicing, post-transcriptional, post-translational modifications, availability of proteomics, genomics editing tools, and chemical biology approaches gave hope to manipulation of cancer cells reversion to a normal cell phenotype as evidences are subtle but definitive. Regardless of the advancement, there is a long way to go, as customized techniques are required to be fine-tuned with precision to attain more insights into tumor reversion. Tumor regression models using available genome-editing methods, followed by in vitro and in vivo proteomics profiling techniques show early evidence. This review summarizes tumor reversion developments, present issues, and unaddressed challenges that remained in the uncharted territory to modulate cellular machinery for tumor reversion towards therapeutic purposes successfully. Ongoing research reaffirms the potential promises of understanding the mechanism of tumor reversion and required refinement that is warranted in vitro and in vivo models of tumor reversion, and the potential translation of these into cancer therapy. Furthermore, therapeutic compounds were reported to induce phenotypic changes in cancer cells into normal cells, which will contribute in understanding the mechanism of tumor reversion. Altogether, the efforts collectively suggest that tumor reversion will likely reveal a new wave of therapeutic discoveries that will significantly impact clinical practice in cancer therapy.

Regulation of RNA Splicing: Aberrant Splicing Regulation and Therapeutic Targets in Cancer

RNA splicing is a critical step in the maturation of precursor mRNA (pre-mRNA) by removing introns and exons. The combination of inclusion and exclusion of introns and exons in pre-mRNA can generate vast diversity in mature mRNA from a limited number of genes. Cancer cells acquire cancer-specific mechanisms through aberrant splicing regulation to acquire resistance to treatment and to promote malignancy. Splicing regulation involves many factors, such as proteins, non-coding RNAs, and DNA sequences at many steps. Thus, the dysregulation of splicing is caused by many factors, including mutations in RNA splicing factors, aberrant expression levels of RNA splicing factors, small nuclear ribonucleoproteins biogenesis, mutations in snRNA, or genomic sequences that are involved in the regulation of splicing, such as 5' and 3' splice sites, branch point site, splicing enhancer/silencer, and changes in the chromatin status that affect the splicing profile. This review focuses on the dysregulation of RNA splicing related to cancer and the associated therapeutic methods.

A pan-cancer transcriptome analysis of exitron splicing identifies novel cancer driver genes and neoepitopes

Exitron splicing (EIS) creates a cryptic intron (called an exitron) within a protein-coding exon to increase proteome diversity. EIS is poorly characterized, but emerging evidence suggests a role for EIS in cancer. Through a systematic investigation of EIS across 33 cancers from 9,599 tumor transcriptomes, we discovered that EIS affected 63% of human coding genes and that 95% of those events were tumor specific. Notably, we observed a mutually exclusive pattern between EIS and somatic mutations in their affected genes. Functionally, we discovered that EIS altered known and novel cancer driver genes for causing gain- or loss-of-function, which promotes tumor progression. Importantly, we identified EIS-derived neoepitopes that bind to major histocompatibility complex (MHC) class I or II. Analysis of clinical data from a clear cell renal cell carcinoma cohort revealed an association between EIS-derived neoantigen load and checkpoint inhibitor response. Our findings establish the importance of considering EIS alterations when nominating cancer driver events and neoantigens.

A cell-to-patient machine learning transfer approach uncovers novel basal-like breast cancer prognostic markers amongst alternative splice variants

BACKGROUND

Breast cancer is amongst the 10 first causes of death in women worldwide. Around 20% of patients are misdiagnosed leading to early metastasis, resistance to treatment and relapse. Many clinical and gene expression profiles have been successfully used to classify breast tumours into 5 major types with different prognosis and sensitivity to specific treatments. Unfortunately, these profiles have failed to subclassify breast tumours into more subtypes to improve diagnostics and survival rate. Alternative splicing is emerging as a new source of highly specific biomarkers to classify tumours in different grades. Taking advantage of extensive public transcriptomics datasets in breast cancer cell lines (CCLE) and breast cancer tumours (TCGA), we have addressed the capacity of alternative splice variants to subclassify highly aggressive breast cancers.

RESULTS

Transcriptomics analysis of alternative splicing events between luminal, basal A and basal B breast cancer cell lines identified a unique splicing signature for a subtype of tumours, the basal B, whose classification is not in use in the clinic yet. Basal B cell lines, in contrast with luminal and basal A, are highly metastatic and express epithelial-to-mesenchymal (EMT) markers, which are hallmarks of cell invasion and resistance to drugs. By developing a semi-supervised machine learning approach, we transferred the molecular knowledge gained from these cell lines into patients to subclassify basal-like triple negative tumours into basal A- and basal B-like categories. Changes in splicing of 25 alternative exons, intimately related to EMT and cell invasion such as ENAH, CD44 and CTNND1, were sufficient to identify the basal-like patients with the worst prognosis. Moreover, patients expressing this basal B-specific splicing signature also expressed newly identified biomarkers of metastasis-initiating cells, like CD36, supporting a more invasive phenotype for this basal B-like breast cancer subtype.

CONCLUSIONS

Using a novel machine learning approach, we have identified an EMT-related splicing signature capable of subclassifying the most aggressive type of breast cancer, which are basal-like triple negative tumours. This proof-of-concept demonstrates that the biological knowledge acquired from cell lines can be transferred to patients data for further clinical investigation. More studies, particularly in 3D culture and organoids, will increase the accuracy of this transfer of knowledge, which will open new perspectives into the development of novel therapeutic strategies and the further identification of specific biomarkers for drug resistance and cancer relapse.

Alternative splicing perturbation landscape identifies RNA binding proteins as potential therapeutic targets in cancer

Alternative splicing (AS) plays an important role in gene regulation, and AS perturbations are frequently observed in cancer. RNA binding protein (RBP) is one of the molecular determinants of AS, and perturbations in RBP-gene network activity are causally associated with cancer development. Here, we performed a systematic analysis to characterize the perturbations in AS events across 18 cancer types. We showed that AS alterations were prevalent in cancer and involved in cancer-related pathways. Given that the extent of AS perturbation was associated with disease severity, we proposed a computational pipeline to identify RBP regulators. Pan-cancer analysis identified a number of conserved RBP regulators, which play important roles in regulating AS of genes involved in cancer hallmark pathways. Our application analysis revealed that the expression of 68 RBP regulators helped in cancer subtyping. Specifically, we identified four subtypes of kidney cancer with differences in cancer hallmark pathway activities and prognosis. Finally, we identified the small molecules that can potentially target the RBP genes and suggested potential candidates for cancer therapy. In summary, our comprehensive AS perturbation landscape analysis identified RBPs as potential therapeutic targets in cancer and provided novel insights into the regulatory functions of RBPs in cancer.

RNA Dysregulation: An Expanding Source of Cancer Immunotherapy Targets

Cancer transcriptomes frequently exhibit RNA dysregulation. As the resulting aberrant transcripts may be translated into cancer-specific proteins, there is growing interest in exploiting RNA dysregulation as a source of tumor antigens (TAs) and thus novel immunotherapy targets. Recent advances in high-throughput technologies and rapid accumulation of multiomic cancer profiling data in public repositories have provided opportunities to systematically characterize RNA dysregulation in cancer and identify antigen targets for immunotherapy. However, given the complexity of cancer transcriptomes and proteomes, important conceptual and technological challenges exist. Here, we highlight the expanding repertoire of TAs arising from RNA dysregulation and introduce multiomic and big data strategies for identifying optimal immunotherapy targets. We discuss extant barriers for translating these targets into effective therapies as well as the implications for future research.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Circular RNA Sequencing of Maternal Platelets: A Novel Tool for the Identification of Pregnancy-Specific Biomarkers

BACKGROUND

In the first trimester of pregnancy, the maternal platelet is directly involved in a positive feedback mechanism that facilitates invasion of the extravillous trophoblast into the maternal spiral arteries. Dysfunctional trophoblast invasion with defective deep placentation is primordial in the etiology of the "great obstetrical syndromes."

METHODS

In this proof-of-concept study, using transcriptome analysis of circular RNA (circRNA) following RNA sequencing of maternal platelets, we tested whether pregnancy-specific circRNA markers could be identified in the first trimester of normal pregnancies. Differential transcript expression analysis of circRNAs, as predicted by Accurate CircRNA Finder Suite, CircRNA Identifier (version 2), and Known and Novel Isoform Explorer, was done using thromboSeq.R with variation of multiple settings. Test performance was checked for (a) de novo circRNA identification using the novel platelet-specific Plt-circR4 as a positive control, (b) complete segregation of groups (pregnant vs nonpregnant) after heat map-dendrogram clustering, (c) identification of pregnancy-specific circRNA markers at a false discovery rate (FDR) <0.05, and (d) confirmation of differentially expressed circRNA markers with an FDR <0.05 by an independent method, reverse transcription-quantitative PCR.

RESULTS

Of the differentially expressed circRNAs with P values <0.05, 41 circRNAs were upregulated (logFC >2), and 52 circRNAs were downregulated (logFC less than -2) in first-trimester platelet RNA. Of these, nuclear receptor-interacting protein 1 circRNA covering exons 2 and 3 of the 5'-untranslated region was pregnancy specific with upregulation in first-trimester maternal platelets compared to nonpregnant controls.

CONCLUSION

CircRNA sequencing of first-trimester maternal platelets permits the identification of novel pregnancy-specific RNA biomarkers. Future use could include the assessment of maternal and fetal well-being.

Bisphenol A Exposure Changes the Transcriptomic and Proteomic Dynamics of Human Retinoblastoma Y79 Cells

Bisphenol A (BPA) is a xenoestrogen chemical commonly used to manufacture polycarbonate plastics and epoxy resin and might affect various human organs. However, the cellular effects of BPA on the eyes have not been widely investigated. This study aimed to investigate the cellular cytotoxicity by BPA exposure on human retinoblastoma cells. BPA did not show cytotoxic effects, such as apoptosis, alterations to cell viability and cell cycle regulation. Comparative analysis of the transcriptome and proteome profiles were investigated after long-term exposure of Y79 cells to low doses of BPA. Transcriptome analysis using RNA-seq revealed that mRNA expression of the post-transcriptional regulation-associated gene sets was significantly upregulated in the BPA-treated group. Cell cycle regulation-associated gene sets were significantly downregulated by exposure to BPA. Interestingly, RNA-seq analysis at the transcript level indicated that alternative splicing events, particularly retained introns, were noticeably altered by low-dose BPA treatment. Additionally, proteome profiling using MALDI-TOF-MS identified a total of nine differentially expressed proteins. These results suggest that alternative splicing events and altered gene/protein expression patterns are critical phenomena affected by long-term low-dose BPA exposure. This represents a novel marker for the detection of various diseases associated with environmental pollutants such as BPA.

Comprehensive characterisation of intronic mis-splicing mutations in human cancers

Previous studies studying mis-splicing mutations were based on exome data and thus our current knowledge is largely limited to exons and the canonical splice sites. To comprehensively characterise intronic mis-splicing mutations, we analysed 1134 pan-cancer whole genomes and transcriptomes together with 3022 normal control samples. The ratio-based splicing analysis resulted in 678 somatic intronic mutations, with 46% residing in deep introns. Among the 309 deep intronic single nucleotide variants, 245 altered core splicing codes, with 38% activating cryptic splice sites, 12% activating cryptic polypyrimidine tracts, and 36% and 12% disrupting authentic polypyrimidine tracts and branchpoints, respectively. All the intronic cryptic splice sites were created at pre-existing GT/AG dinucleotides or by GC-to-GT conversion. Notably, 85 deep intronic mutations indicated gain of splicing enhancers or loss of splicing silencers. We found that 64 tumour suppressors were affected by intronic mutations and blood cancers showed higher proportion of deep intronic mutations. In particular, a telomere maintenance gene, POT1, was recurrently mis-spliced by deep intronic mutations in blood cancers. We validated a pseudoexon activation involving a splicing silencer in POT1 by CRISPR/Cas9. Our results shed light on previously unappreciated mechanisms by which noncoding mutations acting on splicing codes in deep introns contribute to tumourigenesis.

The Duchenne muscular dystrophy gene and cancer

BACKGROUND

Mutation of the Duchenne muscular dystrophy (DMD) gene causes Duchenne and Becker muscular dystrophy, degenerative neuromuscular disorders that primarily affect voluntary muscles. However, increasing evidence implicates DMD in the development of all major cancer types. DMD is a large gene with 79 exons that codes for the essential muscle protein dystrophin. Alternative promotor usage drives the production of several additional dystrophin protein products with roles that extend beyond skeletal muscle. The importance and function(s) of these gene products outside of muscle are not well understood.

CONCLUSIONS

We highlight a clear role for DMD in the pathogenesis of several cancers, including sarcomas, leukaemia's, lymphomas, nervous system tumours, melanomas and various carcinomas. We note that the normal balance of DMD gene products is often disrupted in cancer. The short dystrophin protein Dp71 is, for example, typically maintained in cancer whilst the full-length Dp427 gene product, a likely tumour suppressor, is frequently inactivated in cancer due to a recurrent loss of 5' exons. Therefore, the ratio of short and long gene products may be important in tumorigenesis. In this review, we summarise the tumours in which DMD is implicated and provide a hypothesis for possible mechanisms of tumorigenesis, although the question of cause or effect may remain. We hope to stimulate further study into the potential role of DMD gene products in cancer and the development of novel therapeutics that target DMD.

Targeting RNA-binding proteins in acute and chronic leukemia

RNA-binding proteins (RBPs) play a crucial role in cellular physiology by regulating RNA processing, translation, and turnover. In neoplasms, RBP support of cancer-relevant expression of alternatively spliced, modified, and stabilized mRNA transcripts is essential to self-renewal, proliferation, and adaptation to stress. In this review, we assess the impact of key families of RBPs in leukemogenesis, review progress in targeting those proteins with small molecules, and discuss how multilevel composition of posttranscriptional regulation of gene expression could be used for potential therapies in acute and chronic leukemia.

SCISSOR: a framework for identifying structural changes in RNA transcripts

High-throughput sequencing protocols such as RNA-seq have made it possible to interrogate the sequence, structure and abundance of RNA transcripts at higher resolution than previous microarray and other molecular techniques. While many computational tools have been proposed for identifying mRNA variation through differential splicing/alternative exon usage, challenges in its analysis remain. Here, we propose a framework for unbiased and robust discovery of aberrant RNA transcript structures using short read sequencing data based on shape changes in an RNA-seq coverage profile. Shape changes in selecting sample outliers in RNA-seq, SCISSOR, is a series of procedures for transforming and normalizing base-level RNA sequencing coverage data in a transcript independent manner, followed by a statistical framework for its analysis ( https://github.com/hyochoi/SCISSOR ). The resulting high dimensional object is amenable to unsupervised screening of structural alterations across RNA-seq cohorts with nearly no assumption on the mutational mechanisms underlying abnormalities. This enables SCISSOR to independently recapture known variants such as splice site mutations in tumor suppressor genes as well as novel variants that are previously unrecognized or difficult to identify by any existing methods including recurrent alternate transcription start sites and recurrent complex deletions in 3' UTRs.

Genetic Etiology of Left-Sided Obstructive Heart Lesions: A Story in Development

Congenital heart disease is the most common congenital defect observed in newborns. Within the spectrum of congenital heart disease are left‐sided obstructive lesions (LSOLs), which include hypoplastic left heart syndrome, aortic stenosis, bicuspid aortic valve, coarctation of the aorta, and interrupted aortic arch. These defects can arise in isolation or as a component of a defined syndrome; however, nonsyndromic defects are often observed in multiple family members and associated with high sibling recurrence risk. This clear evidence for a heritable basis has driven a lengthy search for disease‐causing variants that has uncovered both rare and common variants in genes that, when perturbed in cardiac development, can result in LSOLs. Despite advancements in genetic sequencing platforms and broadening use of exome sequencing, the currently accepted LSOL‐associated genes explain only 10% to 20% of patients. Further, the combinatorial effects of common and rare variants as a cause of LSOLs are emerging. In this review, we highlight the genes and variants associated with the different LSOLs and discuss the strengths and weaknesses of the present genetic associations. Furthermore, we discuss the research avenues needed to bridge the gaps in our current understanding of the genetic basis of nonsyndromic congenital heart disease.

Newly Identified Members of FGFR1 Splice Variants Engage in Cross-talk with AXL/AKT Axis in Salivary Adenoid Cystic Carcinoma

Adenoid cystic carcinoma (ACC) is the second most common malignancy of the salivary gland. Although characterized as an indolent tumor, ACC often leads to incurable metastatic disease. Patients with ACC respond poorly to currently available therapeutic drugs and factors contributing to the limited response remain unknown. Determining the role of molecular alterations frequently occurring in ACC may clarify ACC tumorigenesis and advance the development of effective treatment strategies. Applying Splice Expression Variant Analysis and outlier statistics on RNA sequencing of primary ACC tumors and matched normal salivary gland tissues, we identified multiple alternative splicing events (ASE) of genes specific to ACC. In ACC cells and patient-derived xenografts, FGFR1 was a uniquely expressed ASE. Detailed PCR analysis identified three novel, truncated, intracellular domain-lacking FGFR1 variants (FGFR1v). Cloning and expression analysis suggest that the three FGFR1v are cell surface proteins, that expression of FGFR1v augmented pAKT activity, and that cells became more resistant to pharmacologic FGFR1 inhibitor. FGFR1v-induced AKT activation was associated with AXL function, and inhibition of AXL activity in FGFR1v knockdown cells led to enhanced cytotoxicity in ACC. Moreover, cell killing effect was increased by dual inhibition of AXL and FGFR1 in ACC cells. This study demonstrates that these previously undescribed FGFR1v cooperate with AXL and desensitize cells to FGFR1 inhibitor, which supports further investigation into combined FGFR1 and AXL inhibition as an effective ACC therapy.This study identifies several FGFR1 variants that function through the AXL/AKT signaling pathway independent of FGF/FGFR1, desensitizing cells to FGFR1 inhibitor suggestive of a potential resistance mechanism in ACC. SIGNIFICANCE: This study identifies several FGFR1 variants that function through the AXL/AKT signaling pathway independent of FGF/FGFR1, desensitizing cells to FGFR1 inhibitor, suggestive of a potential resistance mechanism in ACC.

Alternative splicing in aging and Alzheimer's disease: Highlighting the role of tau and estrogen receptor α isoforms in the hypothalamus

Human genes show the highest efficacy of alternative splicing (AS) in the brain as compared to other tissues. Within the brain, a remarkably rich diversity of AS events was identified in the hypothalamus. The AS frequency is increased in the aging brain. Such AS events, as intron retention and accumulation of circular RNAs, were acknowledged as some of the main hallmarks of the aging brain. In Alzheimer's disease (AD) pivotal (tau gene, in particular), risk, candidate and other genes show significant alterations in AS. Therefore AD has been suggested to be a disease of dysregulated AS. One of the reported risk factors for AD is estrogen deficiency that may interfere with the extension of neurobrillary tangles. Mounting evidence suggests that estrogens may decrease hyperphosphorylated tau deposition in the brain. Furthermore, AS of estrogen receptor α (ERα) mRNA is decreased in AD brain areas with the highest tau load. These potential interactions among tau, estrogens, and ERα AS may be important for the development of therapeutic and preventive strategies for AD. The intriguing point is that the amount of splice variants of ERα in the hypothalamus and the hippocampus is increased in aging and decreased in AD, while ERα is one of the regulators of AS and is subject to AS itself.

meCLICK-Seq, a Substrate-Hijacking and RNA Degradation Strategy for the Study of RNA Methylation

The fates of RNA species in a cell are controlled by ribonucleases, which degrade them by exploiting the universal structural 2'-OH group. This phenomenon plays a key role in numerous transformative technologies, for example, RNA interference and CRISPR/Cas13-based RNA editing systems. These approaches, however, are genetic or oligomer-based and so have inherent limitations. This has led to interest in the development of small molecules capable of degrading nucleic acids in a targeted manner. Here we describe click-degraders, small molecules that can be covalently attached to RNA species through click-chemistry and can degrade them, that are akin to ribonucleases. By using these molecules, we have developed the meCLICK-Seq (methylation CLICK-degradation Sequencing) a method to identify RNA modification substrates with high resolution at intronic and intergenic regions. The method hijacks RNA methyltransferase activity to introduce an alkyne, instead of a methyl, moiety on RNA. Subsequent copper(I)-catalyzed azide-alkyne cycloaddition reaction with the click-degrader leads to RNA cleavage and degradation exploiting a mechanism used by endogenous ribonucleases. Focusing on N6-methyladenosine (m6A), meCLICK-Seq identifies methylated transcripts, determines RNA methylase specificity, and reliably maps modification sites in intronic and intergenic regions. Importantly, we show that METTL16 deposits m6A to intronic polyadenylation (IPA) sites, which suggests a potential role for METTL16 in IPA and, in turn, splicing. Unlike other methods, the readout of meCLICK-Seq is depletion, not enrichment, of modified RNA species, which allows a comprehensive and dynamic study of RNA modifications throughout the transcriptome, including regions of low abundance. The click-degraders are highly modular and so may be exploited to study any RNA modification and design new technologies that rely on RNA degradation.

Prevalent intron retention fine-tunes gene expression and contributes to cellular senescence

Intron retention (IR) is the least well‐understood alternative splicing type in animals, and its prevalence and function in physiological and pathological processes have long been underestimated. Cellular senescence contributes to individual aging and age‐related diseases and can also serve as an important cancer prevention mechanism. Dynamic IR events have been observed in senescence models and aged tissues; however, whether and how IR impacts senescence remain unclear. Through analyzing polyA⁺ RNA‐seq data from human replicative senescence models, we found IR was prevalent and dynamically regulated during senescence and IR changes negatively correlated with expression alteration of corresponding genes. We discovered that knocking down (KD) splicing factor U2AF1, which showed higher binding density to retained introns and decreased expression during senescence, led to senescence‐associated phenotypes and global IR changes. Intriguingly, U2AF1‐KD‐induced IR changes also negatively correlated with gene expression. Furthermore, we demonstrated that U2AF1‐mediated IR of specific gene (CPNE1 as an example) contributed to cellular senescence. Decreased expression of U2AF1, higher IR of CPNE1, and reduced expression of CPNE1 were also discovered in dermal fibroblasts with age. We discovered prevalent IR could fine‐tune gene expression and contribute to senescence‐associated phenotypes, largely extending the biological significance of IR.

Identification of driver genes with aberrantly alternative splicing events in pediatric patients with retinoblastoma

Retinoblastoma (RB) is one of the most common cancer in children. However, the specific mechanism about RB tumorigenesis has not been fully understood. In this study, to comprehensively characterize the splicing alterations in the tumorigenesis of RB, we analyzed the differential alternative splicing events in RB. Specifically, the isoforms of RB1 were downregulated in the RB samples, and a large proportion of differentially expressed genes had multiple differentially expressed transcripts (64%). We identified 1453 genes with differential alternative splicing, among which, SE accounted for the majority, followed by MXE, RI, A3SS, and A5SS. Furthermore, the biological function related to the normal function of eyes, and E2F family TFs were significantly enriched by the genes with differential alternative splicing. Among the genes associated with visual sense, ABCA4 was found to have two mutually exclusive exons, resulting in two isoforms with different functionalities. Notably, DAZAP1 was identified as one of the critical splicing factors in RB, which was potentially involved in E2F and RB pathways. Functionally, differential binding sites in DAZAP1 protein were significantly observed between RB and normal samples. Based on the comprehensive analysis of the differential alternative splicing events and splicing factors, we identified some driver genes with differential alternative splicing and critical splicing factors involved in RB, which would greatly improve our understanding of the alternative splicing process in the tumorigenesis of RB.

Prognostic alternative splicing signature in cervical squamous cell carcinoma

Basing on alternative splicing events (ASEs) databases, the authors herein aim to explore potential prognostic biomarkers for cervical squamous cell carcinoma (CESC). mRNA expression profiles and relevant clinical data of 223 patients with CESC were obtained from The Cancer Genome Atlas (TCGA). Correlated genes, ASEs and percent-splice-in (PSI) were downloaded from SpliceSeq, respectively. The PSI values of survival-associated alternative splicing events (SASEs) were used to construct the basis of a prognostic index (PI). A protein-protein interaction (PPI) network of genes related to SASEs was generated by STRING and analysed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Consequently, 41,776 ASEs were discovered in 19,724 genes, 2596 of which linked with 3669 SASEs. The PPI network of SASEs related genes revealed that TP53 and UBA52 were core genes. The low-risk group had a longer survival period than high-risk counterparts, both groups being defined according to PI constructed upon the top 20 splicing events or PI on the overall splicing events. The AUC value of ROC reached up to 0.88, demonstrating the prognostic potential of PI in CESC. These findings suggested that ASEs involve in the pathogenesis of CESC and may serve as promising prognostic biomarkers for this female malignancy.