Widespread intron retention diversifies most cancer transcriptomes

Distinct Minor Splicing Patterns across Cancers

In human cells, the U12 spliceosome, also known as the minor spliceosome, is responsible for the splicing of 0.5% of introns, while the major U2 spliceosome is responsible for the other 99.5%. While many studies have been done to characterize and understand splicing dysregulation in cancer, almost all of them have focused on U2 splicing and ignored U12 splicing, despite evidence suggesting minor splicing is involved in cell cycle regulation. In this study, we analyzed RNA-seq data from The Cancer Genome Atlas for 14 different cohorts to determine differential splicing of minor introns in tumor and adjacent normal tissue. We found that in some cohorts, such as breast cancer, there was a strong skew towards minor introns showing increased splicing in the tumor; in others, such as the renal chromophobe cell carcinoma cohort, the opposite pattern was found, with minor introns being much more likely to have decreased splicing in the tumor. Further analysis of gene expression did not reveal any candidate regulatory mechanisms that could cause these different minor splicing phenotypes between cohorts. Our data suggest context-dependent roles of the minor spliceosome in tumorigenesis and provides a foundation for further investigation of minor splicing in cancer, which could then serve as a basis for novel therapeutic strategies.

Intron-Retention Neoantigen Load Predicts Favorable Prognosis in Pancreatic Cancer

PURPOSE

High tumor mutation burden (TMB) in many cancer types is associated with the production of tumor-specific neoantigens, a favorable outcome and response to immune checkpoint blockade (ICB) therapy. Besides mutation-derived neoantigens, aberrant intron retention also produces tumor neopeptides that could trigger an immune response. The relationship between intron-retention-derived tumor neoantigens (IR-neoAg) and clinical outcomes in pancreatic cancer remains uncertain. Here, we quantify IR-neoAg in pancreatic cancer and evaluate whether IR-neoAg load might serve as a biomarker for selecting patients who may benefit from ICB therapy.

METHODS

We developed a computational approach to estimate patient-specific IR-neoAg load from transcriptome data available in The Cancer Genome Atlas pancreatic cancer cohort. Associations between IR-neoAg load and patient overall survival were evaluated using Kaplan-Meier estimates and Cox regression. Differential expression of immune checkpoint and HLA-I genes was evaluated in tumors with high IR-neoAg load.

RESULTS

High IR-neoAg load predicted better overall survival in pancreatic cancer, although no association was found for TMB. IR-neoAg load remained a significant prognostic factor after adjusting for patient age, sex, tumor stage and grade, and TMB. Moreover, pancreatic tumors with both high IR-neoAg load and high HLA-I gene expression had similar gene expression profiles as other tumor types that showed response to anti-programmed cell death protein 1 therapy.

CONCLUSION

IR-neoAg load is associated with favorable survival in pancreatic cancer. These findings provide strong evidence for considering IR-neoAgs when selecting patients who might benefit from ICB therapy.

The role of splicing factor PRPF8 in breast cancer

BACKGROUND:

Alternative splicing is a mechanism to produce different proteins with diverse functions from one gene. Many splicing factors play an important role in cancer progression. PRPF8 is a core protein component of the spliceosome complex, U4/U6-U5 tri-snRNP.

OBJECTIVE:

However, PRPF8 involved in mRNA alternative splicing are rarely included in the prognosis.

METHODS:

We found that PRPF8 was expressed in all examined cancer types. Further analyses found that PRPF8 expression was significantly different between the breast cancer and paracancerous tissues.

RESULTS:

Survival analyses showed that PRPF8-high patients had a poor prognosis, and the expression of PRPF8 is associated with distant metastasis-free survival (DMFS) and post progression survival (PPS). Gene Set Enrichment Analysis (GSEA) has revealed that PRPF8 expression is correlated with TGF-β, JAK-STAT, and cell cycle control pathways. Consistent with these results, upon PRPF8 silencing, the growth of MCF-7 cells was reduced, the ability of cell clone formation was weakened, and p⁢21 expression was increased.

CONCLUSIONS:

These results have revealed that PRPF8 is a significant factor for splicing in breast cancer progression.

Type I and II PRMTs inversely regulate post-transcriptional intron detention through Sm and CHTOP methylation

Protein arginine methyltransferases (PRMTs) are required for the regulation of RNA processing factors. Type I PRMT enzymes catalyze mono- and asymmetric dimethylation; Type II enzymes catalyze mono- and symmetric dimethylation. To understand the specific mechanisms of PRMT activity in splicing regulation, we inhibited Type I and II PRMTs and probed their transcriptomic consequences. Using the newly developed Splicing Kinetics and Transcript Elongation Rates by Sequencing (SKaTER-seq) method, analysis of co-transcriptional splicing demonstrated that PRMT inhibition resulted in altered splicing rates. Surprisingly, co-transcriptional splicing kinetics did not correlate with final changes in splicing of polyadenylated RNA. This was particularly true for retained introns (RI). By using actinomycin D to inhibit ongoing transcription, we determined that PRMTs post-transcriptionally regulate RI. Subsequent proteomic analysis of both PRMT-inhibited chromatin and chromatin-associated polyadenylated RNA identified altered binding of many proteins, including the Type I substrate, CHTOP, and the Type II substrate, SmB. Targeted mutagenesis of all methylarginine sites in SmD3, SmB, and SmD1 recapitulated splicing changes seen with Type II PRMT inhibition, without disrupting snRNP assembly. Similarly, mutagenesis of all methylarginine sites in CHTOP recapitulated the splicing changes seen with Type I PRMT inhibition. Examination of subcellular fractions further revealed that RI were enriched in the nucleoplasm and chromatin. Taken together, these data demonstrate that, through Sm and CHTOP arginine methylation, PRMTs regulate the post-transcriptional processing of nuclear, detained introns.

Towards resolution of the intron retention paradox in breast cancer

BACKGROUND

After many years of neglect in the field of alternative splicing, the importance of intron retention (IR) in cancer has come into focus following landmark discoveries of aberrant IR patterns in cancer. Many solid and liquid tumours are associated with drastic increases in IR, and such patterns have been pursued as both biomarkers and therapeutic targets. Paradoxically, breast cancer (BrCa) is the only tumour type in which IR is reduced compared to adjacent normal breast tissue.

METHODS

In this study, we have conducted a pan-cancer analysis of IR with emphasis on BrCa and its subtypes. We explored mechanisms that could cause aberrant and pathological IR and clarified why normal breast tissue has unusually high IR.

RESULTS

Strikingly, we found that aberrantly decreasing IR in BrCa can be largely attributed to normal breast tissue having the highest occurrence of IR events compared to other healthy tissues. Our analyses suggest that low numbers of IR events in breast tumours are associated with poor prognosis, particularly in the luminal B subtype. Interestingly, we found that IR frequencies negatively correlate with cell proliferation in BrCa cells, i.e. rapidly dividing tumour cells have the lowest number of IR events. Aberrant RNA-binding protein expression and changes in tissue composition are among the causes of aberrantly decreasing IR in BrCa.

CONCLUSIONS

Our results suggest that IR should be considered for therapeutic manipulation in BrCa patients with aberrantly low IR levels and that further work is needed to understand the cause and impact of high IR in other tumour types.

The RNA-binding protein IGF2BP3 is critical for MLL-AF4-mediated leukemogenesis

Despite recent advances in therapeutic approaches, patients with MLL-rearranged leukemia still have poor outcomes. Here, we find that the RNA-binding protein IGF2BP3, which is overexpressed in MLL-translocated leukemia, strongly amplifies MLL-Af4-mediated leukemogenesis. Deletion of Igf2bp3 significantly increases the survival of mice with MLL-Af4-driven leukemia and greatly attenuates disease, with a minimal impact on baseline hematopoiesis. At the cellular level, MLL-Af4 leukemia-initiating cells require Igf2bp3 for their function in leukemogenesis. At the molecular level, IGF2BP3 regulates a complex posttranscriptional operon governing leukemia cell survival and proliferation. IGF2BP3-targeted mRNA transcripts include important MLL-Af4-induced genes, such as those in the Hoxa locus, and the Ras signaling pathway. Targeting of transcripts by IGF2BP3 regulates both steady-state mRNA levels and, unexpectedly, pre-mRNA splicing. Together, our findings show that IGF2BP3 represents an attractive therapeutic target in this disease, providing important insights into mechanisms of posttranscriptional regulation in leukemia.

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.

Integrated multi-omics analysis of the clinical relevance and potential regulatory mechanisms of splicing factors in hepatocellular carcinoma

Splicing factors (SFs) have been increasingly documented to perturb the genome of cancers. However, little is known about the alterations of SFs in hepatocellular carcinoma (HCC). This study comprehensively delineated the genomic and epigenomic characteristics of 404 SFs in HCC based on the multi-omics data from the Cancer Genome Atlas database. The analysis revealed several clinically relevant SFs that could be effective biomarkers for monitoring the onset and prognosis of HCC (such as, HSPB1, DDX39A, and NELFE, which were the three most significant clinically relevant SFs). Functional enrichment analysis of these indicators showed the enrichment of pathways related to splicing and mRNA processes. Furthermore, the study found that SF copy number variation is common in HCC and could be a typical characteristic of hepato-carcinogenesis; the complex expression regulation of SFs was significantly affected by copy number variant and methylation. Several SFs with significant mutation patterns were identified (such as, RNF213, SF3B1, SPEN, NOVA1, and EEF1A1), and the potential regulatory network of SFs was constructed to identify their potential mechanisms for regulating clinically relevant alternative splicing events. Therefore, this study established a foundation to uncover the broad molecular spectrum of SFs for future functional and therapeutic studies of HCC.

Oncogenic dysregulation of pre-mRNA processing by protein kinases: challenges and therapeutic opportunities

Alternative splicing and polyadenylation represent two major steps in pre-mRNA-processing, which ensure proper gene expression and diversification of human transcriptomes. Deregulation of these processes contributes to oncogenic programmes involved in the onset, progression and evolution of human cancers, which often result in the acquisition of resistance to existing therapies. On the other hand, cancer cells frequently increase their transcriptional rate and develop a transcriptional addiction, which imposes a high stress on the pre-mRNA-processing machinery and establishes a therapeutically exploitable vulnerability. A prominent role in fine-tuning pre-mRNA-processing mechanisms is played by three main families of protein kinases: serine arginine protein kinase (SRPK), CDC-like kinase (CLK) and cyclin-dependent kinase (CDK). These kinases phosphorylate the RNA polymerase, splicing factors and regulatory proteins involved in cleavage and polyadenylation of the nascent transcripts. The activity of SRPKs, CLKs and CDKs can be altered in cancer cells, and their inhibition was shown to exert anticancer effects. In this review, we describe key findings that have been reported on these topics and discuss challenges and opportunities of developing therapeutic approaches targeting splicing factor kinases.

Comprehensive Analysis of Splicing Factor and Alternative Splicing Event to Construct Subtype-Specific Prognosis-Predicting Models for Breast Cancer

Recent evidence suggests that splicing factors (SFs) and alternative splicing (AS) play important roles in cancer progression. We constructed four SF-risk-models using 12 survival-related SFs. In Luminal-A, Luminal-B, Her-2, and Basal-Like BRCA, SF-risk-models for three genes (PAXBP1, NKAP, and NCBP2), four genes (RBM15B, PNN, ACIN1, and SRSF8), three genes (LSM3, SNRNP200, and SNU13), and three genes (SRPK3, PUF60, and PNN) were constructed. These models have a promising prognosis-predicting power. The co-expression and protein-protein interaction analysis suggest that the 12 SFs are highly functional-connected. Pathway analysis and gene set enrichment analysis suggests that the functional role of the selected 12 SFs is highly context-dependent among different BRCA subtypes. We further constructed four AS-risk-models with good prognosis predicting ability in four BRCA subtypes by integrating the four SF-risk-models and 21 survival-related AS-events. This study proposed that SFs and ASs were potential multidimensional biomarkers for the diagnosis, prognosis, and treatment of BRCA.

ISOTOPE: ISOform-guided prediction of epiTOPEs in cancer

Immunotherapies provide effective treatments for previously untreatable tumors and identifying tumor-specific epitopes can help elucidate the molecular determinants of therapy response. Here, we describe a pipeline, ISOTOPE (ISOform-guided prediction of epiTOPEs In Cancer), for the comprehensive identification of tumor-specific splicing-derived epitopes. Using RNA sequencing and mass spectrometry for MHC-I associated proteins, ISOTOPE identified neoepitopes from tumor-specific splicing events that are potentially presented by MHC-I complexes. Analysis of multiple samples indicates that splicing alterations may affect the production of self-epitopes and generate more candidate neoepitopes than somatic mutations. Although there was no difference in the number of splicing-derived neoepitopes between responders and non-responders to immune therapy, higher MHC-I binding affinity was associated with a positive response. Our analyses highlight the diversity of the immunogenic impacts of tumor-specific splicing alterations and the importance of studying splicing alterations to fully characterize tumors in the context of immunotherapies. ISOTOPE is available at https://github.com/comprna/ISOTOPE.

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.

Splicing factor gene mutations in acute myeloid leukemia offer additive value if incorporated in current risk classification

Splicing factor (SF) mutations are important contributors to the pathogenesis of hematological malignancies; however, their relevance in risk classification of acute myeloid leukemia (AML) warrants further investigation. To gain more insight into the characteristics of patients with AML carrying SF mutations, we studied their association with clinical features, cytogenetic and molecular abnormalities, and clinical outcome in a large cohort of 1447 patients with AML and high-risk myelodysplastic syndrome. SF mutations were identified in 22% of patients and were associated with multiple unfavorable clinical features, such as older age, antecedent myeloid disorders, and adverse risk factors (mutations in RUNX1 and ASXL1). Furthermore, they had significantly shorter event-free and overall survival. Notably, in European LeukemiaNet (ELN) 2017 favorable- and intermediate-risk groups, SF3B1 mutations were indicative of relatively poor prognosis. In addition, patients carrying concomitant SF mutations and RUNX1 mutations had a particularly adverse prognosis. In patients without any of the 4 most common SF mutations, RUNX1 mutations were associated with relatively good outcome, which was comparable to that of intermediate-risk patients. In this study, we propose that SF mutations be considered for incorporation into prognostic classification systems. First, SF3B1 mutations could be considered an intermediate prognostic factor when co-occurring with favorable risk features and as an adverse prognostic factor for patients currently categorized as having intermediate risk, according to the ELN 2017 classification. Second, the prognostic value of the current adverse factor RUNX1 mutations seems to be limited to its co-occurrence with SF mutations.

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.

Maternal cecal microbiota transfer rescues early-life antibiotic-induced enhancement of type 1 diabetes in mice

Early-life antibiotic exposure perturbs the intestinal microbiota and accelerates type 1 diabetes (T1D) development in the NOD mouse model. Here, we found that maternal cecal microbiota transfer (CMT) to NOD mice after early-life antibiotic perturbation largely rescued the induced T1D enhancement. Restoration of the intestinal microbiome was significant and persistent, remediating the antibiotic-depleted diversity, relative abundance of particular taxa, and metabolic pathways. CMT also protected against perturbed metabolites and normalized innate and adaptive immune effectors. CMT restored major patterns of ileal microRNA and histone regulation of gene expression. Further experiments suggest a gut-microbiota-regulated T1D protection mechanism centered on Reg3γ, in an innate intestinal immune network involving CD44, TLR2, and Reg3γ. This regulation affects downstream immunological tone, which may lead to protection against tissue-specific T1D injury.

NSrp70 is a lymphocyte-essential splicing factor that controls thymocyte development

Alternative pre-mRNA splicing is a critical step to generate multiple transcripts, thereby dramatically enlarging the proteomic diversity. Thus, a common feature of most alternative splicing factor knockout models is lethality. However, little is known about lineage-specific alternative splicing regulators in a physiological setting. Here, we report that NSrp70 is selectively expressed in developing thymocytes, highest at the double-positive (DP) stage. Global splicing and transcriptional profiling revealed that NSrp70 regulates the cell cycle and survival of thymocytes by controlling the alternative processing of various RNA splicing factors, including the oncogenic splicing factor SRSF1. A conditional-knockout of Nsrp1 (NSrp70-cKO) using CD4Cre developed severe defects in T cell maturation to single-positive thymocytes, due to insufficient T cell receptor (TCR) signaling and uncontrolled cell growth and death. Mice displayed severe peripheral lymphopenia and could not optimally control tumor growth. This study establishes a model to address the function of lymphoid-lineage-specific alternative splicing factor NSrp70 in a thymic T cell developmental pathway.

The RNA-binding protein LUC7L2 mediates MITA/STING intron retention to negatively regulate innate antiviral response

MITA (also known as STING) is an ER-located adaptor protein, which mediates DNA-triggered innate immune response and is critically involved in autoimmune diseases and tumorigenesis. MITA is regulated by post-translational modifications, but how post-transcriptional mechanisms are involved in the regulation of MITA is still largely unknown. Here, we identified the RNA-binding protein LUC7L2 as a negative regulator of DNA virus-triggered innate immune response. LUC7L2-deficient mice exhibited resistance to lethal herpes simplex virus 1 (HSV-1) infection and reduced HSV-1 loads in the brain. Mechanistically, LUC7L2 directly bound to intron 3 of MITA precursor messenger RNA, inhibited its splicing and promoted its nonsense-mediated decay, leading to its downregulation at protein level. LUC7L2-deficient cells had markedly increased MITA level, leading to heightened innate antiviral response. Finally, LUC7L2 was induced following HSV-1 infection. Our findings reveal a feedback negative post-transcriptional regulatory mechanism for regulation of MITA-mediated innate immune response to viral and aberrant cellular DNA.

HER2-PI9 and HER2-I12: two novel and functionally active splice variants of the oncogene HER2 in breast cancer

In this study, two novel alternative splice variants of HER2, named HER2-PI9 and HER2-I12, were identified in breast cancer cell lines and breast tumour tissues. Whilst HER2-P19 arises from the inclusion of an 117 bp cassette-exon of intron 9 of HER2, HER2-I12 results from intron 12 inclusion. In silico analyses were performed to predict the amino acid sequences of these two HER2 novel variants. To confirm their protein expression, plasmid vectors were generated and transfected into the HER2 negative breast cancer cell line, MCF-7. Additionally, their functional properties in oncogenic signalling were confirmed. Expression of HER2-PI9 and HER2-I12 was successful and matched the in silico predictions. Importantly, these splice variants can modulate the phosphorylation levels of extracellular signal-related kinase 1/2 (ERK1/2) and Akt/protein kinase B (Akt) signalling in MCF-7 breast cancer cells. Enhanced cellular proliferation, migration and invasion were observed in the case of the HER2-I12 expressing model. In human tissues and breast carcinoma tumours both variants were present. This study reveals two novel splice variants of HER2. Additionally, the potential biological activity for HER2-PI9 and HER2-I12 in breast cancer cells is also reported..

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.

Intron retention as a new pre-symptomatic marker of aging and its recovery to the normal state by a traditional Japanese multi-herbal medicine

Intron retention (IR) is an important regulatory mechanism that affects gene expression and protein functions. Using klotho mice at the pre-symptomatic state, we discovered that retained-introns accumulated in several organs including the liver and that among these retained introns in the liver a subset was recovered to the normal state by a Japanese traditional herbal medicine. This is the first report of IR recovery by a medicine. IR-recovered genes fell into two categories: those involved in liver-specific metabolism and in splicing. Metabolome analysis of the liver showed that the klotho mice were under starvation stress. In addition, our differentially expressed gene analysis showed that liver metabolism was actually recovered by the herbal medicine at the transcriptional level. By analogy with the widespread accumulation of intron-retained pre-mRNAs induced by heat shock stress, we propose a model in which retained-introns in klotho mice were induced by an aging stress and in which this medicine-related IR recovery is indicative of the actual recovery of liver-specific metabolic function to the healthy state. Accumulation of retained-introns was also observed at the pre-symptomatic state of aging in wild-type mice and may be an excellent marker for this state in general.

Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity

RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.

Targeting Epigenetics to Cure HIV-1: Lessons From (and for) Cancer Treatment

The Human Immunodeficiency Virus type 1 (HIV-1) integrates in the host genome as a provirus resulting in a long-lived reservoir of infected CD4 cells. As a provirus, HIV-1 has several aspects in common with an oncogene. Both the HIV-1 provirus and oncogenes only cause disease when expressed. A successful cure of both cancer and HIV-1 includes elimination of all cells with potential to regenerate the disease. For over two decades, epigenetic drugs developed against cancer have been used in the HIV-1 field to modulate the state of the proviral chromatin. Cells with an intact HIV-1 provirus exist in three states of infection: productive, inducible latent, and non-inducible latent. Here focus is on HIV-1, transcription control and chromatin structure; how the inducible proviruses are maintained in a chromatin structure that allows reactivation of transcription; and how transcription switches between different stages to allow for an abundance of different transcripts from a single promoter. Recently it was shown that a functional cure of HIV can be achieved by encapsulating all intact HIV-1 proviruses in heterochromatin, giving hope that epigenetic interventions may be used to end the HIV-1 epidemic.

Targeting the transcription cycle and RNA processing in cancer treatment

The transcriptional program and RNA splicing machinery are highly and frequently dysregulated in human cancers due to genomic and epigenomic alterations during tumorigenesis. This leads to cancer-specific dependencies on components of the transcriptional program and RNA splicing machinery, providing alternative and targetable 'Achilles' heels' for cancer treatment in the clinic. To target these vulnerabilities in cancer cells, potent and specific transcriptional CDK inhibitors and chemical compounds that impair splicing have been developed and evaluated in preclinical cancer models. Several novel combination approaches with immune or targeted therapies have also been proposed for cancer treatment. More recently, inhibitors targeting transcriptional CDKs, splicing, or PRMT5 have shown promising therapeutic potential in preclinical studies, and many of them have rapidly advanced into early clinical trials for treatment of human cancer.

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.

Tracking pre-mRNA maturation across subcellular compartments identifies developmental gene regulation through intron retention and nuclear anchoring

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Splice it up: Atypical transcripts to boost leukemia immunotherapy

Neoantigens are prime targets for cancer immunotherapy, but their identification in low mutational burden malignancies remains challenging. In this issue of Immunity, Ehx et al. show that atypical transcripts, and particularly retained introns, expand the spectrum of leukemia immunotherapy targets.

Atypical acute myeloid leukemia-specific transcripts generate shared and immunogenic MHC class-I-associated epitopes

Acute myeloid leukemia (AML) has not benefited from innovative immunotherapies, mainly because of the lack of actionable immune targets. Using an original proteogenomic approach, we analyzed the major histocompatibility complex class I (MHC class I)-associated immunopeptidome of 19 primary AML samples and identified 58 tumor-specific antigens (TSAs). These TSAs bore no mutations and derived mainly (86%) from supposedly non-coding genomic regions. Two AML-specific aberrations were instrumental in the biogenesis of TSAs, intron retention, and epigenetic changes. Indeed, 48% of TSAs resulted from intron retention and translation, and their RNA expression correlated with mutations of epigenetic modifiers (e.g., DNMT3A). AML TSA-coding transcripts were highly shared among patients and were expressed in both blasts and leukemic stem cells. In AML patients, the predicted number of TSAs correlated with spontaneous expansion of cognate T cell receptor clonotypes, accumulation of activated cytotoxic T cells, immunoediting, and improved survival. These TSAs represent attractive targets for AML immunotherapy.

easyCLIP analysis of RNA-protein interactions incorporating absolute quantification

Quantitative criteria to identify proteins as RNA-binding proteins (RBPs) are presently lacking, as are criteria to define RBP target RNAs. Here, we develop an ultraviolet (UV) cross-linking immunoprecipitation (CLIP)-sequencing method, easyCLIP. easyCLIP provides absolute cross-link rates, as well as increased simplicity, efficiency, and capacity to visualize RNA libraries during sequencing library preparation. Measurement of >200 independent cross-link experiments across >35 proteins identifies an RNA cross-link rate threshold that distinguishes RBPs from non-RBPs and defines target RNAs as those with a complex frequency unlikely for a random protein. We apply easyCLIP to the 33 most recurrent cancer mutations across 28 RBPs, finding increased RNA binding per RBP molecule for KHDRBS2 R168C, A1CF E34K and PCBP1 L100P/Q cancer mutations. Quantitating RBP-RNA interactions can thus nominate proteins as RBPs and define the impact of specific disease-associated RBP mutations on RNA association.

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.

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.

Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

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

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

Alternative Splicing and Cleavage of GLUT8

The GLUT (SLC2) family of membrane-associated transporters are described as glucose transporters. However, this family is divided into three classes and, though the regulated transporter activity of class I proteins is becoming better understood, class III protein functions continue to be obscure. We have cataloged the relative expression and splicing of SLC2 mRNA isomers in tumors and normal tissues, with a focus on breast tumors and cell lines. mRNA for the class III protein GLUT8 is the predominant SLC2 species expressed alongside GLUT1 in many tissues, but GLUT8 mRNA exists mostly as an untranslated splice form in tumors. We confirm that GLUT8 is not presented at the cell surface and does not transport glucose directly. However, we reveal a lysosome-dependent reaction that cleaves the GLUT8 protein and releases the carboxy-terminal peptide to a separate vesicle population. Given the localization of GLUT8 at a major metabolic hub (the late endosomal/lysosomal interface) and its regulated cleavage reaction, we evaluated TXNIP-mediated hexosamine homeostasis and speculate that GLUT8 may function as a sensory component of this reaction.

Widespread Aberrant Alternative Splicing despite Molecular Remission in Chronic Myeloid Leukaemia Patients

Simple Summary

This study provides new insights into the changing transcriptomic and epigenomic landscapes in chronic myeloid leukaemia (CML) patients who are receiving tyrosine kinase inhibitor (TKI) therapy (often life-long). Alternative splicing, vital for cellular homeostasis, is dysregulated in human cancers. Remarkably, we found abnormal splicing patterns despite molecular remission in peripheral blood cells of chronic-phase CML patients. This phenomenon is independent of the TKI drug used and in striking contrast to the normalisation of gene expression and DNA methylation patterns.

Abstract

Vast transcriptomics and epigenomics changes are characteristic of human cancers, including leukaemia. At remission, we assume that these changes normalise so that omics-profiles resemble those of healthy individuals. However, an in-depth transcriptomic and epigenomic analysis of cancer remission has not been undertaken. A striking exemplar of targeted remission induction occurs in chronic myeloid leukaemia (CML) following tyrosine kinase inhibitor (TKI) therapy. Using RNA sequencing and whole-genome bisulfite sequencing, we profiled samples from chronic-phase CML patients at diagnosis and remission and compared these to healthy donors. Remarkably, our analyses revealed that abnormal splicing distinguishes remission samples from normal controls. This phenomenon is independent of the TKI drug used and in striking contrast to the normalisation of gene expression and DNA methylation patterns. Most remarkable are the high intron retention (IR) levels that even exceed those observed in the diagnosis samples. Increased IR affects cell cycle regulators at diagnosis and splicing regulators at remission. We show that aberrant splicing in CML is associated with reduced expression of specific splicing factors, histone modifications and reduced DNA methylation. Our results provide novel insights into the changing transcriptomic and epigenomic landscapes of CML patients during remission. The conceptually unanticipated observation of widespread aberrant alternative splicing after remission induction warrants further exploration. These results have broad implications for studying CML relapse and treating minimal residual disease.

Intron retention is a robust marker of intertumoral heterogeneity in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer with a 5-year survival rate of <8%. Unsupervised clustering of 76 PDAC patients based on intron retention (IR) events resulted in two clusters of tumors (IR-1 and IR-2). While gene expression-based clusters are not predictive of patient outcome in this cohort, the clusters we developed based on intron retention were associated with differences in progression-free interval. IR levels are lower and clinical outcome is worse in IR-1 compared with IR-2. Oncogenes were significantly enriched in the set of 262 differentially retained introns between the two IR clusters. Higher IR levels in IR-2 correlate with higher gene expression, consistent with detention of intron-containing transcripts in the nucleus in IR-2. Out of 258 genes encoding RNA-binding proteins (RBP) that were differentially expressed between IR-1 and IR-2, the motifs for seven RBPs were significantly enriched in the 262-intron set, and the expression of 25 RBPs were highly correlated with retention levels of 139 introns. Network analysis suggested that retention of introns in IR-2 could result from disruption of an RBP protein-protein interaction network previously linked to efficient intron removal. Finally, IR-based clusters developed for the majority of the 20 cancer types surveyed had two clusters with asymmetrical distributions of IR events like PDAC, with one cluster containing mostly intron loss events. Taken together, our findings suggest IR may be an important biomarker for subclassifying tumors.

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.

Alternative Splicing: Expanding the Landscape of Cancer Biomarkers and Therapeutics

Alternative splicing (AS) is a critical post-transcriptional regulatory mechanism used by more than 95% of transcribed human genes and responsible for structural transcript variation and proteome diversity. In the past decade, genome-wide transcriptome sequencing has revealed that AS is tightly regulated in a tissue- and developmental stage-specific manner, and also frequently dysregulated in multiple human cancer types. It is currently recognized that splicing defects, including genetic alterations in the spliced gene, altered expression of both core components or regulators of the precursor messenger RNA (pre-mRNA) splicing machinery, or both, are major drivers of tumorigenesis. Hence, in this review we provide an overview of our current understanding of splicing alterations in cancer, and emphasize the need to further explore the cancer-specific splicing programs in order to obtain new insights in oncology. Furthermore, we also discuss the recent advances in the identification of dysregulated splicing signatures on a genome-wide scale and their potential use as biomarkers. Finally, we highlight the therapeutic opportunities arising from dysregulated splicing and summarize the current approaches to therapeutically target AS in cancer.

HER2 splice variants in breast cancer: investigating their impact on diagnosis and treatment outcomes

Overexpression of the HER2 receptor occurs in approximately 20% of breast cancer patients. HER2 positivity is associated with poor prognosis and aggressive tumour phenotypes, which led to rapid progress in HER2 targeted therapeutics and diagnostic testing. Whilst these advances have greatly increased patients' chances of survival, resistance to HER2 targeted therapies, be that intrinsic or acquired, remains a problem. Different forms of the HER2 protein exist within tumours in tandem and can display altered biological activities. Interest in HER2 variants in breast cancer increased when links between resistance to anti-HER2 therapies and a particular variant, Δ16-HER2, were identified. Moreover, the P100 variant potentially reduces the efficacy of the anti-HER2 therapy trastuzumab. Another variant, Herstatin, exhibits 'auto-inhibitory' behaviour. More recently, new HER2 variants have been identified and are currently being assessed for their pro- and anti-cancer properties. It is important when directing the care of patients to consider HER2 variants collectively. This review considers HER2 variants in the context of the tumour environment where multiple variants are co-expressed at altered ratios. This study also provides an up to date account of the landscape of HER2 variants and links this to patterns of resistance against HER2 therapies and treatment plans.

Intron retention and its impact on gene expression and protein diversity: A review and a practical guide

Intron retention (IR) occurs when a complete and unspliced intron remains in mature mRNA. An increasing body of literature has demonstrated a major role for IR in numerous biological functions, including several that impact human health and disease. Although experimental technologies used to study other forms of mRNA splicing can also be used to investigate IR, a specialized downstream computational analysis is optimal for IR discovery and analysis. Here we provide a review of IR and its biological implications, as well as a practical guide for how to detect and analyze it. Several methods, including long read third generation direct RNA sequencing, are described. We have developed an R package, FakIR, to facilitate the execution of the bioinformatic tasks recommended in this review and a tutorial on how to fit them to users aims. Additionally, we provide guidelines and experimental protocols to validate IR discovery and to evaluate the potential impact of IR on gene expression and protein output. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Splicing Regulation/Alternative Splicing RNA Methods > RNA Analyses in vitro and In Silico.

Widespread intron retention impairs protein homeostasis in C9orf72 ALS brains

The GGGGCC hexanucleotide expansion in C9orf72 (C9) is the most frequent known cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD), yet a clear understanding of how C9 fits into the broader context of ALS/FTD pathology has remained lacking. The repetitive RNA derived from the C9 repeat is known to sequester hnRNPH, a splicing regulator, into insoluble aggregates, resulting in aberrant alternative splicing. Furthermore, hnRNPH insolubility and altered splicing of a robust set of targets have been observed to correlate in C9 and sporadic ALS/FTD patients alike, suggesting that changes along this axis are a core feature of disease pathogenesis. Here, we characterize previously uncategorized RNA splicing defects involving widespread intron retention affecting almost 2000 transcripts in C9ALS/FTD brains exhibiting a high amount of sequestered, insoluble hnRNPH. These intron retention events appear not to alter overall expression levels of the affected transcripts but rather the protein-coding regions. These retained introns affect transcripts in multiple cellular pathways predicted to be involved in C9 as well as sporadic ALS/FTD etiology, including the proteasomal and autophagy systems. The retained intron pre-mRNAs display a number of characteristics, including enrichment of hnRNPH-bound splicing enhancer motifs and a propensity for G-quadruplex (G-Q) formation, linking the defective splicing directly to high amounts of sequestered hnRNPH. Together, our results reveal previously undetected splicing defects in high insoluble hnRNPH-associated C9ALS brains, suggesting a feedback between effective RNA-binding protein dosage and protein quality control in C9, and perhaps all, ALS/FTD.

Holding on to Junk Bonds: Intron Retention in Cancer and Therapy

Intron retention (IR) in cancer was for a long time overlooked by the scientific community, as it was previously considered to be an artifact of a dysfunctional spliceosome. Technological advancements made in the last decade offer unique opportunities to explore the role of IR as a widespread phenomenon that contributes to the transcriptional diversity of many cancers. Numerous studies in cancer have shed light on dysregulation of cellular mechanisms that lead to aberrant and pathologic IR. IR is not merely a mechanism of gene regulation, but rather it can mediate cancer pathogenesis and therapeutic resistance in various human diseases. The burden of IR in cancer is governed by perturbations to mechanisms known to regulate this phenomenon and include epigenetic variation, mutations within the gene body, and splicing factor dysregulation. This review summarizes possible causes for aberrant IR and discusses the role of IR in therapy or as a consequence of disease treatment. As neoepitopes originating from retained introns can be presented on the cancer cell surface, the development of personalized cancer vaccines based on IR-derived neoepitopes should be considered. Ultimately, a deeper comprehension about the origins and consequences of aberrant IR may aid in the development of such personalized cancer vaccines.

The role of alternative splicing in cancer: From oncogenesis to drug resistance

Alternative splicing is a tightly regulated process whereby non-coding sequences of pre-mRNA are removed and protein-coding segments are assembled in diverse combinations, ultimately giving rise to proteins with distinct or even opposing functions. In the past decade, whole genome/transcriptome sequencing studies revealed the high complexity of splicing regulation, which occurs co-transcriptionally and is influenced by chromatin status and mRNA modifications. Consequently, splicing profiles of both healthy and malignant cells display high diversity and alternative splicing was shown to be widely deregulated in multiple cancer types. In particular, mutations in pre-mRNA regulatory sequences, splicing regulators and chromatin modifiers, as well as differential expression of splicing factors are important contributors to cancer pathogenesis. It has become clear that these aberrations contribute to many facets of cancer, including oncogenic transformation, cancer progression, response to anticancer drug treatment as well as resistance to therapy. In this respect, alternative splicing was shown to perturb the expression a broad spectrum of relevant genes involved in drug uptake/metabolism (i.e. SLC29A1, dCK, FPGS, and TP), activation of nuclear receptor pathways (i.e. GR, AR), regulation of apoptosis (i.e. MCL1, BCL-X, and FAS) and modulation of response to immunotherapy (CD19). Furthermore, aberrant splicing constitutes an important source of novel cancer biomarkers and the spliceosome machinery represents an attractive target for a novel and rapidly expanding class of therapeutic agents. Small molecule inhibitors targeting SF3B1 or splice factor kinases were highly cytotoxic against a wide range of cancer models, including drug-resistant cells. Importantly, these effects are enhanced in specific cancer subsets, such as splicing factor-mutated and c-MYC-driven tumors. Furthermore, pre-clinical studies report synergistic effects of spliceosome modulators in combination with conventional antitumor agents. These strategies based on the use of low dose splicing modulators could shift the therapeutic window towards decreased toxicity in healthy tissues. Here we provide an extensive overview of the latest findings in the field of regulation of splicing in cancer, including molecular mechanisms by which cancer cells harness alternative splicing to drive oncogenesis and evade anticancer drug treatment as well as splicing-based vulnerabilities that can provide novel treatment opportunities. Furthermore, we discuss current challenges arising from genome-wide detection and prediction methods of aberrant splicing, as well as unravelling functional relevance of the plethora of cancer-related splicing alterations.

A new mechanism for a familiar mutation - bovine DGAT1 K232A modulates gene expression through multi-junction exon splice enhancement

BACKGROUND

The DGAT1 gene encodes an enzyme responsible for catalysing the terminal reaction in mammary triglyceride synthesis, and underpins a well-known pleiotropic quantitative trait locus (QTL) with a large influence on milk composition phenotypes. Since first described over 15 years ago, a protein-coding variant K232A has been assumed as the causative variant underlying these effects, following in-vitro studies that demonstrated differing levels of triglyceride synthesis between the two protein isoforms.

RESULTS

We used a large RNAseq dataset to re-examine the underlying mechanisms of this large milk production QTL, and hereby report novel expression-based functions of the chr14 g.1802265AA > GC variant that encodes the DGAT1 K232A substitution. Using expression QTL (eQTL) mapping, we demonstrate a highly-significant mammary eQTL for DGAT1, where the K232A mutation appears as one of the top associated variants for this effect. By conducting in vitro expression and splicing experiments in bovine mammary cell culture, we further show modulation of splicing efficiency by this mutation, likely through disruption of an exon splice enhancer as a consequence of the allele encoding the 232A variant.

CONCLUSIONS

The relative contributions of the enzymatic and transcription-based mechanisms now attributed to K232A remain unclear; however, these results suggest that transcriptional impacts contribute to the diversity of lactation effects observed at the DGAT1 locus.

Androgen receptor signaling regulates the transcriptome of prostate cancer cells by modulating global alternative splicing

Androgen receptor (AR), is a transcription factor and a member of a hormone receptor superfamily. AR plays a vital role in the progression of prostate cancer and is a crucial target for therapeutic interventions. While the majority of advanced-stage prostate cancer patients will initially respond to the androgen deprivation, the disease often progresses to castrate-resistant prostate cancer (CRPC). Interestingly, CRPC tumors continue to depend on hyperactive AR signaling and will respond to potent second-line antiandrogen therapies, including bicalutamide (CASODEX^®) and enzalutamide (XTANDI^®). However, the progression-free survival rate for the CRPC patients on antiandrogen therapies is only 8–19 months. Hence, there is a need to understand the mechanisms underlying CRPC progression and eventual treatment resistance. Here, we have leveraged next-generation sequencing and newly developed analytical methodologies to evaluate the role of AR signaling in regulating the transcriptome of prostate cancer cells. The genomic and pharmacologic stimulation and inhibition of AR activity demonstrates that AR regulates alternative splicing within cancer-relevant genes. Furthermore, by integrating transcriptomic data from in vitro experiments and in prostate cancer patients, we found that a significant number of AR-regulated splicing events are associated with tumor progression. For example, we found evidence for an inadvertent AR-antagonist-mediated switch in IDH1 and PL2G2A isoform expression, which is associated with a decrease in overall survival of patients. Mechanistically, we discovered that the epithelial-specific splicing regulators (ESRP1 and ESRP2), flank many AR-regulated alternatively spliced exons. And, using 2D invasion assays, we show that the inhibition of ESRPs can suppress AR-antagonist-driven tumor invasion. Our work provides evidence for a new mechanism by which AR alters the transcriptome of prostate cancer cells by modulating alternative splicing. As such, our work has important implications for CRPC progression and development of resistance to treatment with bicalutamide and enzalutamide.

Ctcf haploinsufficiency mediates intron retention in a tissue-specific manner

CTCF is a master regulator of gene transcription and chromatin organisation with occupancy at thousands of DNA target sites genome-wide. While CTCF is essential for cell survival, CTCF haploinsufficiency is associated with tumour development and hypermethylation. Increasing evidence demonstrates CTCF as a key player in several mechanisms regulating alternative splicing (AS), however, the genome-wide impact of Ctcf dosage on AS has not been investigated.

We examined the effect of Ctcf haploinsufficiency on gene expression and AS in five tissues from Ctcf hemizygous (Ctcf^+/-) mice. Reduced Ctcf levels caused distinct tissue-specific differences in gene expression and AS in all tissues. An increase in intron retention (IR) was observed in Ctcf^+/- liver and kidney. In liver, this specifically impacted genes associated with cytoskeletal organisation, splicing and metabolism. Strikingly, most differentially retained introns were short, with a high GC content and enriched in Ctcf binding sites in their proximal upstream genomic region. This study provides new insights into the effects of CTCF haploinsufficiency on organ transcriptomes and the role of CTCF in AS regulation.

Spliceosomal factor mutations and mis-splicing in MDS

Somatic, heterozygous missense and nonsense mutations in at least seven proteins that function in the spliceosome are found at high frequency in MDS patients. These proteins act at various steps in the process of splicing by the spliceosome and lead to characteristic alterations in the alternative splicing of a subset of genes. Several studies have investigated the effects of these mutations and have attempted to identify a commonly affected gene or pathway. Here, we summarize what is known about the normal function of these proteins and how the mutations alter the splicing landscape of the genome. We also summarize the commonly mis-spliced gene targets and discuss the state of mechanistic unification that has been achieved. Finally, we discuss alternative mechanisms by which these mutations may lead to disease.