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

Predicting Nonsense-mediated mRNA Decay from Splicing Events in Sepsis using RNA-Sequencing Data

Alternative splicing (AS) and nonsense-mediated mRNA decay (NMD) are highly conserved cellular mechanisms that modulate gene expression. Here we introduce NMD pipeline that computes how splicing events introduce premature termination codons to mRNA transcripts via frameshift, then predicts the rate of PTC-dependent NMD. We utilize whole blood, deep RNA-sequencing data from critically ill patients to study gene expression in sepsis. Statistical significance was determined as adjusted p value < 0.05 and |log2foldchange| > 2 for differential gene expression and probability >= 0.9 and |DeltaPsi| > 0.1 for AS. NMD pipeline was developed based on AS data from Whippet. We demonstrate that the rate of NMD is higher in sepsis and deceased groups compared to control and survived groups, which signify purposeful downregulation of transcripts by AS-NMD or aberrant splicing due to altered physiology. Predominance of non-exon skipping events was associated with disease and mortality states. The NMD pipeline also revealed proteins with potential novel roles in sepsis. Together, these results emphasize the utility of NMD pipeline in studying AS-NMD along with differential gene expression and discovering potential protein targets in sepsis.

Decoding RAP1 's Role in Yeast mRNA Splicing

Messenger RNA (mRNA) splicing is a fundamental and tightly regulated process in eukaryotes, where the spliceosome removes non-coding sequences from pre-mRNA to produce mature mRNA for protein translation. Alternative splicing enables the generation of multiple RNA isoforms and protein products from a single gene, regulating both isoform diversity and abundance. While splicing is widespread in eukaryotes, only ∼3% of genes in Saccharomyces cerevisiae undergo splicing, with most containing a single intron. However, intron-containing genes, primarily ribosomal protein genes, are highly expressed and constitute about one-third of the total mRNA pool. These genes are transcriptionally regulated by Repressor Activator Protein 1 ( RAP1 ), prompting us to investigate whether RAP1 influences mRNA splicing. Using RNA sequencing, we identified a novel role for RAP1 in alternative splicing, particularly in intron retention (IR) while minor effects were observed on alternative 3' and 5' splice site usage. Many IR-containing transcripts introduced premature termination codons, likely leading to degradation via nonsense-mediated decay (NMD). Consistent with previous literature, genes with predicted NMD in our study also had reduced overall expression levels suggesting that RAP1 plays an important role in this understudied mechanism of gene expression regulation.

Germline-specific RNA helicase DDX4 forms cytoplasmic granules in cancer cells and promotes tumor growth

Cancer cells undergo major epigenetic alterations and transcriptomic changes, including ectopic expression of tissue- and cell-type-specific genes. Here, we show that the germline-specific RNA helicase DDX4 forms germ-granule-like cytoplasmic ribonucleoprotein granules in various human tumors, but not in cultured cancer cells. These cancerous DDX4 complexes contain RNA-binding proteins and splicing regulators, including many known germ granule components. The deletion of DDX4 in cancer cells induces transcriptomic changes and affects the alternative splicing landscape of a number of genes involved in cancer growth and invasiveness, leading to compromised capability of DDX4-null cancer cells to form xenograft tumors in immunocompromised mice. Importantly, the occurrence of DDX4 granules is associated with poor survival in patients with head and neck squamous cell carcinoma and higher histological grade of prostate cancer. Taken together, these results show that the germ-granule-resembling cancerous DDX4 granules control gene expression and promote malignant and invasive properties of cancer cells.

Histone H3K9 Lactylation Confers Temozolomide Resistance in Glioblastoma via LUC7L2-Mediated MLH1 Intron Retention

Temozolomide (TMZ) resistance remains the major obstacle in the treatment of glioblastoma (GBM). Lactylation is a novel post-translational modification that is involved in various tumors. However, whether lactylation plays a role in GBM TMZ resistance remains unclear. Here it is found that histone H3K9 lactylation (H3K9la) confers TMZ resistance in GBM via LUC7L2-mediated intron 7 retention of MLH1. Mechanistically, lactylation is upregulated in recurrent GBM tissues and TMZ-resistant cells, and is mainly concentrated in histone H3K9. Combined multi-omics analysis, including CUT&Tag, SLAM-seq, and RNA-seq, reveals that H3K9 lactylation is significantly enriched in the LUC7L2 promoter and activates LUC7L2 transcription to promote its expression. LUC7L2 mediates intron 7 retention of MLH1 to reduce MLH1 expression, and thereby inhibit mismatch repair (MMR), ultimately leading to GBM TMZ resistance. Of note, it is identified that a clinical anti-epileptic drug, stiripentol, which can cross the blood-brain barrier and inhibit lactate dehydrogenase A/B (LDHA/B) activity, acts as a lactylation inhibitor and renders GBM cells more sensitive to TMZ in vitro and in vivo. These findings not only shed light on the mechanism of lactylation in GBM TMZ resistance but also provide a potential combined therapeutic strategy for clinical GBM treatment.

LINC00460-FUS-MYC feedback loop drives breast cancer metastasis and doxorubicin resistance

Therapeutic resistance and metastasis largely contribute to mortality from breast cancer and therefore understanding the underlying mechanisms of such remains an urgent challenge. By cross-analysis of TCGA and GEO databases, LINC00460 was identified as an oncogenic long non-coding RNA, highly expressed in Doxorubicin resistant breast cancer. LINC00460 was further demonstrated to promote stem cell-like and epithelial-mesenchymal transition (EMT) characteristics in breast cancer cells. LINC00460 interacts with FUS protein with consequent enhanced stabilization, which further promotes MYC mRNA maturation. LINC00460 expression was transcriptionally enhanced by c-MYC protein, forming a positive feedback loop to promote metastasis and Doxorubicin resistance. LINC00460 depletion in Doxorubicin-resistant breast cancer cells restored sensitivity to Doxorubicin and increased the efficacy of c-MYC inhibitor therapy. Collectively, these findings implicate LINC00460 as a promising prognostic biomarker and potential therapeutic target to overcome Doxorubicin resistance in breast cancer.

Advancements and hurdles in the development of a vaccine for triple-negative breast cancer: A comprehensive review of multi-omics and immunomics strategies

Triple-Negative Breast Cancer (TNBC) presents a significant challenge in oncology due to its aggressive behavior and limited therapeutic options. This review explores the potential of immunotherapy, particularly vaccine-based approaches, in addressing TNBC. It delves into the role of immunoinformatics in creating effective vaccines against TNBC. The review first underscores the distinct attributes of TNBC and the importance of tumor antigens in vaccine development. It then elaborates on antigen detection techniques such as exome sequencing, HLA typing, and RNA sequencing, which are instrumental in identifying TNBC-specific antigens and selecting vaccine candidates. The discussion then shifts to the in-silico vaccine development process, encompassing antigen selection, epitope prediction, and rational vaccine design. This process merges computational simulations with immunological insights. The role of Artificial Intelligence (AI) in expediting the prediction of antigens and epitopes is also emphasized. The review concludes by encapsulating how Immunoinformatics can augment the design of TNBC vaccines, integrating tumor antigens, advanced detection methods, in-silico strategies, and AI-driven insights to advance TNBC immunotherapy. This could potentially pave the way for more targeted and efficacious treatments.

Comprehensive analysis of alternative splicing across multiple transcriptomic cohorts reveals prognostic signatures in prostate cancer

BACKGROUND

Alternative splicing (AS) plays a crucial role in transcriptomic diversity and is a hallmark of cancer that profoundly influences the development and progression of prostate cancer (PCa), a prevalent and potentially life-limiting cancer among men. Accumulating evidence has highlighted the association between AS dysregulation and the onset and progression of PCa. However, a comprehensive and integrative analysis of AS profiles at the event level, utilising data from multiple high-throughput cohorts and evaluating the prognosis of PCa progression, remains lacking and calls for thorough exploration.

RESULTS

We identified a differentially expressed retained intron event in ZWINT across three distinct cohorts, encompassing an original array-based dataset profiled by us previously and two RNA sequencing (RNA-seq) datasets. Subsequent in-depth analyses of these RNA-seq datasets revealed  141 altered events, of which 21 demonstrated a significant association with patients' biochemical recurrence-free survival (BCRFS). We formulated an AS event-based prognostic signature, capturing six pivotal events in genes CYP4F12, NFATC4, PIGO, CYP3A5, ALS2CL, and FXYD3. This signature effectively differentiated  high-risk patients diagnosed with PCa, who experienced shorter BCRFS, from their low-risk counterparts. Notably, the signature's predictive power surpassed traditional clinicopathological markers in forecasting 5-year BCRFS, demonstrating robust performance in both internal and external validation sets. Lastly, we constructed a novel nomogram that integrates patients' Gleason scores with pathological tumour stages, demonstrating improved prognostication of BCRFS.

CONCLUSIONS

Prediction of clinical progression remains elusive in PCa. This research uncovers novel splicing events associated with BCRFS, augmenting existing prognostic tools, thus potentially refining clinical decision-making.

Identification of Intron Retention in the Slc16a3 Gene Transcript Encoding the Transporter MCT4 in the Brain of Aged and Alzheimer-Disease Model (APPswePS1dE9) Mice

The monocarboxylate transporter 4 (MCT4; Slc16a3) is expressed in the central nervous system, notably by astrocytes. It is implicated in lactate release and the regulation of glycolytic flux. Whether its expression varies during normal and/or pathological aging is unclear. As the presence of its mature transcript in the brain of young and old mice was determined, an unexpectedly longer RT-PCR fragment was detected in the mouse frontal cortex and hippocampus at 12 vs. 3 months of age. Cultured astrocytes expressed the expected 516 base pair (bp) fragment but treatment with IL-1β to mimic inflammation as can occur during aging led to the additional expression of a 928 bp fragment like that seen in aged mice. In contrast, cultured pericytes (a component of the blood-brain barrier) only exhibited the 516 bp fragment. Intriguingly, cultured endothelial cells constitutively expressed both fragments. When RT-PCR was performed on brain subregions of an Alzheimer mouse model (APPswePS1dE9), no fragment was detected at 3 months, while only the 928 bp fragment was present at 12 months. Sequencing of MCT4 RT-PCR products revealed the presence of a remaining intron between exon 2 and 3, giving rise to the longer fragment detected by RT-PCR. These results unravel the existence of intron retention for the MCT4 gene in the central nervous system. Such alternative splicing appears to increase with age in the brain and might be prominent in neurodegenerative diseases such as Alzheimer's disease. Hence, further studies in vitro and in vivo of intron 2 retention in the Slc16a3 gene transcript are required for adequate characterization concerning the biological roles of Slc16a3 isoforms in the context of aging and Alzheimer's disease pathology.

The application of long-read sequencing in clinical settings

Long-read DNA sequencing technologies have been rapidly evolving in recent years, and their ability to assess large and complex regions of the genome makes them ideal for clinical applications in molecular diagnosis and therapy selection, thereby providing a valuable tool for precision medicine. In the third-generation sequencing duopoly, Oxford Nanopore Technologies and Pacific Biosciences work towards increasing the accuracy, throughput, and portability of long-read sequencing methods while trying to keep costs low. These trades have made long-read sequencing an attractive tool for use in research and clinical settings. This article provides an overview of current clinical applications and limitations of long-read sequencing and explores its potential for point-of-care testing and health care in remote settings.

Splicing Characterization and Isoform Switch Events in Human Keratinocytes Carrying Oncogenes from High-Risk HPV-16 and Low-Risk HPV-84

Infection of epithelial cells with high-risk HPV (HR-HPV) types, followed by expression of virus oncogenic proteins (E5, E6, and E7), leads to genomic imbalance, suppression of tumor inhibitors, and induction of oncogenes. Low-risk HPV (LR-HPV) may slow the rate at which cervical cancer spreads to an invasive stage since co-infection with LR-HPV is linked to a decreased risk of future invasive cancer than infection with HR-HPV alone. We then propose that cancer-progressing changes may be distinguished through identifying the functional differences between LR-HPV and HR-HPV. Lentiviral strategies were followed to establish HaCaT cells with constitutive expression of HPV oncogenes. RNAseq experiments were designed to analyze the transcriptome modulations caused by each of the E5, E6, and E7 oncogenes of HPV-16 and HPV-84 in HaCaT cells. We identified enhanced RNA degradation, spliceosome, and RNA polymerase pathways related to mRNA processing. ATTS (alternative transcription termination site) was discovered to be more prevalent in cells with HPV-16E5 than HPV-84E5. In HPV-16E6-infected cells, ATTS gain was significantly higher than ATTS loss. Cells with HPV-16E7 had more isoforms with intron retention (IR) than those with HPV-84E7. We identified switches in ADAM10, CLSPN, and RNPS1 that led to greater expression of the coding isoforms in HR-HPV. The results of this work highlight differences between LR-HPV and HR-HPV in mRNA processing. Moreover, crucial cervical cancer-related switch events were detected.

Increased chromatin accessibility facilitates intron retention in specific cell differentiation states

Dynamic intron retention (IR) in vertebrate cells is of widespread biological importance. Aberrant IR is associated with numerous human diseases including several cancers. Despite consistent reports demonstrating that intrinsic sequence features can help introns evade splicing, conflicting findings about cell type- or condition-specific IR regulation by trans-regulatory and epigenetic mechanisms demand an unbiased and systematic analysis of IR in a controlled experimental setting. We integrated matched mRNA sequencing (mRNA-Seq), whole-genome bisulfite sequencing (WGBS), nucleosome occupancy methylome sequencing (NOMe-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) data from primary human myeloid and lymphoid cells. Using these multi-omics data and machine learning, we trained two complementary models to determine the role of epigenetic factors in the regulation of IR in cells of the innate immune system. We show that increased chromatin accessibility, as revealed by nucleosome-free regions, contributes substantially to the retention of introns in a cell-specific manner. We also confirm that intrinsic characteristics of introns are key for them to evade splicing. This study suggests an important role for chromatin architecture in IR regulation. With an increasing appreciation that pathogenic alterations are linked to RNA processing, our findings may provide useful insights for the development of novel therapeutic approaches that target aberrant splicing.

Long-Read Sequencing Annotation of the Transcriptome in DNA-PK Inactivated Cells

The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) with a Ku70/Ku80 heterodimer constitutes the intact DNA-PK kinase, which is an upstream component of the DNA repair machinery that signals the DNA damage, orchestrates the DNA repair, and serves to maintain genome integrity. Beyond its role in DNA damage repair, the DNA-PK kinase is also implicated in transcriptional regulation and RNA metabolism, with an illuminated impact on tumor progression and therapeutic responses. However, the efforts to identify DNA-PK regulated transcriptomes are limited by short-read sequencing to resolve the full complexity of the transcriptome. Therefore, we leveraged the PacBio Single Molecule, Real-Time (SMRT) Sequencing platform to study the transcriptome after DNA-PK inactivation to further underscore the importance of its role in diseases. Our analysis revealed additional novel transcriptome and complex gene structures in the DNA-PK inactivated cells, identifying 8,355 high-confidence new isoforms from 3,197 annotated genes and 523 novel genes. Among them, 380 lncRNAs were identified. We validated these findings using computational approaches and confirmatory transcript quantification with short-read sequencing. Several novel isoforms representing distinct splicing events have been validated through PCR experiments. Our analyses provide novel insights into DNA-PK function in transcriptome regulation and RNA metabolism.

Meta-Analysis Suggests That Intron Retention Can Affect Quantification of Transposable Elements from RNA-Seq Data

Simple Summary

Transposable elements (TEs) are repetitive sequences comprising more than one third of the human genome with the original ability to change their location within the genome. Owing to their repetitive nature, the quantification of TEs results often challenging. RNA-seq is a useful tool for genome-wide TEs quantification, nevertheless it also presents technical issues, including low reads mappability and erroneous quantification derived from the transcription of TEs fragments embedded in canonical transcripts. Fragments derived from TEs are found within the introns of most genes, which led to the hypothesis that intron retention (IR) can affect the unbiased quantification of TEs expression. Performing meta-analysis of public RNA-seq datasets, here we observe that IR can indeed impact the quantification of TEs by increasing the number of reads mapped on intronic TE copies. Our work highlights a correlation between IR and TEs expression measurement by RNA-seq that should be taken into account to achieve reliable TEs quantification, especially in samples characterized by extensive IR, because differential IR might be confused with differential TEs expression.

Abstract

Transposable elements (TEs), also known as “jumping genes”, are repetitive sequences with the capability of changing their location within the genome. They are key players in many different biological processes in health and disease. Therefore, a reliable quantification of their expression as transcriptional units is crucial to distinguish between their independent expression and the transcription of their sequences as part of canonical transcripts. TEs quantification faces difficulties of different types, the most important one being low reads mappability due to their repetitive nature preventing an unambiguous mapping of reads originating from their sequences. A large fraction of TEs fragments localizes within introns, which led to the hypothesis that intron retention (IR) can be an additional source of bias, potentially affecting accurate TEs quantification. IR occurs when introns, normally removed from the mature transcript by the splicing machinery, are maintained in mature transcripts. IR is a widespread mechanism affecting many different genes with cell type-specific patterns. We hypothesized that, in an RNA-seq experiment, reads derived from retained introns can introduce a bias in the detection of overlapping, independent TEs RNA expression. In this study we performed meta-analysis using public RNA-seq data from lymphoblastoid cell lines and show that IR can impact TEs quantification using established tools with default parameters. Reads mapped on intronic TEs were indeed associated to the expression of TEs and influence their correct quantification as independent transcriptional units. We confirmed these results using additional independent datasets, demonstrating that this bias does not appear in samples where IR is not present and that differential TEs expression does not impact on IR quantification. We concluded that IR causes the over-quantification of intronic TEs and differential IR might be confused with differential TEs expression. Our results should be taken into account for a correct quantification of TEs expression from RNA-seq data, especially in samples in which IR is abundant.

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.

Prospects and challenges of cancer systems medicine: from genes to disease networks

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

Neo-Antigen mRNA Vaccines

The interest in therapeutic cancer vaccines has caught enormous attention in recent years due to several breakthroughs in cancer research, among which the finding that successful checkpoint blockade treatments reinvigorate neo-antigen-specific T cells and that successful adoptive cell therapies are directed towards neo-antigens. Neo-antigens are cancer-specific antigens, which develop from somatic mutations in the cancer cell genome that can be highly immunogenic and are not subjected to central tolerance. As the majority of neo-antigens are unique to each patient's cancer, a vaccine technology that is flexible and potent is required to develop personalized neo-antigen vaccines. In vitro transcribed mRNA is such a technology platform and has been evaluated for delivery of neo-antigens to professional antigen-presenting cells both ex vivo and in vivo. In addition, strategies that support the activity of T cells in the tumor microenvironment have been developed. These represent a unique opportunity to ensure durable T cell activity upon vaccination. Here, we comprehensively review recent progress in mRNA-based neo-antigen vaccines, summarizing critical milestones that made it possible to bring the promise of therapeutic cancer vaccines within reach.

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.