Widespread intron retention diversifies most cancer transcriptomes

eIF4E orchestrates mRNA processing, RNA export and translation to modify specific protein production

The eukaryotic translation initiation factor eIF4E acts as a multifunctional factor that simultaneously influences mRNA processing, export, and translation in many organisms. Its multifactorial effects are derived from its capacity to bind to the methyl-7-guanosine cap on the 5'end of mRNAs and thus can act as a cap chaperone for transcripts in the nucleus and cytoplasm. In this review, we describe the multifactorial roles of eIF4E in major mRNA-processing events including capping, splicing, cleavage and polyadenylation, nuclear export and translation. We discuss the evidence that eIF4E acts at two levels to generate widescale changes to processing, export and ultimately the protein produced. First, eIF4E alters the production of components of the mRNA processing machinery, supporting a widescale reprogramming of multiple mRNA processing events. In this way, eIF4E can modulate mRNA processing without physically interacting with target transcripts. Second, eIF4E also physically interacts with both capped mRNAs and components of the RNA processing or translation machineries. Further, specific mRNAs are sensitive to eIF4E only in particular mRNA processing events. This selectivity is governed by the presence of cis-acting elements within mRNAs known as USER codes that recruit relevant co-factors engaging the appropriate machinery. In all, we describe the molecular bases for eIF4E's multifactorial function and relevant regulatory pathways, discuss the basis for selectivity, present a compendium of ~80 eIF4E-interacting factors which play roles in these activities and provide an overview of the relevance of its functions to its oncogenic potential. Finally, we summarize early-stage clinical studies targeting eIF4E in cancer.

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.

PRPF40A induces inclusion of exons in GC-rich regions important for human myeloid cell differentiation

We characterized the regulatory mechanisms and role in human myeloid cell survival and differentiation of PRPF40A, a splicing factor lacking a canonical RNA Binding Domain. Upon PRPF40A knockdown, HL-60 cells displayed increased cell death, decreased proliferation and slight differentiation phenotype with upregulation of immune activation genes. Suggestive of both redundant and specific functions, cell death but not proliferation was rescued by overexpression of its paralog PRPF40B. Transcriptomic analysis revealed the predominant role of PRPF40A as an activator of cassette exon inclusion of functionally relevant splicing events. Mechanistically, the exons exclusively upregulated by PRPF40A are flanked by short and GC-rich introns which tend to localize to nuclear speckles in the nucleus center. These PRPF40A regulatory features are shared with other splicing regulators such as SRRM2, SON, PCBP1/2, and to a lesser extent TRA2B and SRSF2, as a part of a functional network that regulates splicing partly via co-localization in the nucleus.

Aberrant spliceosome activity via elevated intron retention and upregulation and phosphorylation of SF3B1 in chronic lymphocytic leukemia

Splicing factor 3b subunit 1 (SF3B1) is the largest subunit and core component of the spliceosome. Inhibition of SF3B1 was associated with an increase in broad intron retention (IR) on most transcripts, suggesting that IR can be used as a marker of spliceosome inhibition in chronic lymphocytic leukemia (CLL) cells. Furthermore, we separately analyzed exonic and intronic mapped reads on annotated RNA-sequencing transcripts obtained from B cells (n = 98 CLL patients) and healthy volunteers (n = 9). We measured intron/exon ratio to use that as a surrogate for alternative RNA splicing (ARS) and found that 66% of CLL-B cell transcripts had significant IR elevation compared with normal B cells (NBCs) and that correlated with mRNA downregulation and low expression levels. Transcripts with the highest IR levels belonged to biological pathways associated with gene expression and RNA splicing. A >2-fold increase of active pSF3B1 was observed in CLL-B cells compared with NBCs. Additionally, when the CLL-B cells were treated with macrolides (pladienolide-B), a significant decrease in pSF3B1, but not total SF3B1 protein, was observed. These findings suggest that IR/ARS is increased in CLL, which is associated with SF3B1 phosphorylation and susceptibility to SF3B1 inhibitors. These data provide additional support to the relevance of ARS in carcinogenesis and evidence of pSF3B1 participation in this process.

The Many Roads from Alternative Splicing to Cancer: Molecular Mechanisms Involving Driver Genes

Cancer driver genes are either oncogenes or tumour suppressor genes that are classically activated or inactivated, respectively, by driver mutations. Alternative splicing-which produces various mature mRNAs and, eventually, protein variants from a single gene-may also result in driving neoplastic transformation because of the different and often opposed functions of the variants of driver genes. The present review analyses the different alternative splicing events that result in driving neoplastic transformation, with an emphasis on their molecular mechanisms. To do this, we collected a list of 568 gene drivers of cancer and revised the literature to select those involved in the alternative splicing of other genes as well as those in which its pre-mRNA is subject to alternative splicing, with the result, in both cases, of producing an oncogenic isoform. Thirty-one genes fall into the first category, which includes splicing factors and components of the spliceosome and splicing regulators. In the second category, namely that comprising driver genes in which alternative splicing produces the oncogenic isoform, 168 genes were found. Then, we grouped them according to the molecular mechanisms responsible for alternative splicing yielding oncogenic isoforms, namely, mutations in cis splicing-determining elements, other causes involving non-mutated cis elements, changes in splicing factors, and epigenetic and chromatin-related changes. The data given in the present review substantiate the idea that aberrant splicing may regulate the activation of proto-oncogenes or inactivation of tumour suppressor genes and details on the mechanisms involved are given for more than 40 driver genes.

Biological relevance of alternative splicing in hematologic malignancies

Alternative splicing (AS) is a strictly regulated process that generates multiple mRNA variants from a single gene, thus contributing to proteome diversity. Transcriptome-wide sequencing studies revealed networks of functionally coordinated splicing events, which produce isoforms with distinct or even opposing functions. To date, several mechanisms of AS are deregulated in leukemic cells, mainly due to mutations in splicing and/or epigenetic regulators and altered expression of splicing factors (SFs). In this review, we discuss aberrant splicing events induced by mutations affecting SFs (SF3B1, U2AF1, SRSR2, and ZRSR2), spliceosome components (PRPF8, LUC7L2, DDX41, and HNRNPH1), and epigenetic modulators (IDH1 and IDH2). Finally, we provide an extensive overview of the biological relevance of aberrant isoforms of genes involved in the regulation of apoptosis (e. g. BCL-X, MCL-1, FAS, and c-FLIP), activation of key cellular signaling pathways (CASP8, MAP3K7, and NOTCH2), and cell metabolism (PKM).

Identification and Functional Characterization of Alternative Transcripts of LncRNA HNF1A-AS1 and Their Impacts on Cell Growth, Differentiation, Liver Diseases, and in Response to Drug Induction

The long non-coding RNA (lncRNA) hepatocyte nuclear factor-1 alpha (HNF1A) antisense RNA 1 (HNF1A-AS1) is an important lncRNA for liver growth, development, cell differentiation, and drug metabolism. Like many lncRNAs, HNF1A-AS1 has multiple annotated alternative transcripts in the human genome. Several fundamental biological questions are still not solved: (1) How many transcripts really exist in biological samples, such as liver samples and liver cell lines? (2) What are the expression patterns of different alternative HNF1A-AS1 transcripts at different conditions, including during cell growth and development, after exposure to xenobiotics (such as drugs), and in disease conditions, such as metabolic dysfunction-associated steatotic liver disease (MASLD), alcohol-associated liver disease (ALD) cirrhosis, and obesity? (3) Does the siRNA used in previous studies knock down one or multiple transcripts? (4) Do different transcripts have the same or different functions for gene regulation? The presented data confirm the existence of several annotated HNF1A-AS1 transcripts in liver samples and cell lines, but also identify some new transcripts, which are not annotated in the Ensembl genome database. Expression patterns of the identified HNF1A-AS1 transcripts are highly correlated with the cell differentiation of matured hepatocyte-like cells from human embryonic stem cells (hESC), growth and differentiation of HepaRG cells, in response to rifampicin induction, and in various liver disease conditions. The expression levels of the HNF1A-AS1 transcripts are also highly correlated to the expression of cytochrome P450 enzymes, such as CYP3A4, during HepaRG growth, differentiation, and in response to rifampicin induction.

ASAPA: a bioinformatic pipeline based on Iso-Seq that identifies the links among alternative splicing, alternative transcription initiation and alternative polyadenylation

BACKGROUND

Although the events associated with alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI) can be identified by many approaches based on isoform sequencing (Iso-Seq), these analyses are generally independent of each other and the links between these events are still rarely mentioned. However, an interdependency analysis can be achieved because the transcriptional start site, splice sites and polyA site could be simultaneously included in a long, full-length read from Iso-Seq.

RESULTS

We create ASAPA pipeline that enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data. We tested this pipeline using Arabidopsis data and found some interesting results: some adjacent introns tend to be simultaneously spliced or retained; coupling between AS and ATI or APA is limited to the initial or terminal intron; and ATI and APA are potentially linked in some special cases.

CONCLUSION

Our pipeline enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data, which is conducive to a better understanding of transcription landscape generation.

betAS: intuitive analysis and visualization of differential alternative splicing using beta distributions

Next-generation RNA sequencing allows alternative splicing (AS) quantification with unprecedented resolution, with the relative inclusion of an alternative sequence in transcripts being commonly quantified by the proportion of reads supporting it as percent spliced-in (PSI). However, PSI values do not incorporate information about precision, proportional to the respective AS events' read coverage. Beta distributions are suitable to quantify inclusion levels of alternative sequences, using reads supporting their inclusion and exclusion as surrogates for the two distribution shape parameters. Each such beta distribution has the PSI as its mean value and is narrower when the read coverage is higher, facilitating the interpretability of its precision when plotted. We herein introduce a computational pipeline, based on beta distributions accurately modeling PSI values and their precision, to quantitatively and visually compare AS between groups of samples. Our methodology includes a differential splicing significance metric that compromises the magnitude of intergroup differences, the estimation uncertainty in individual samples, and the intragroup variability, being therefore suitable for multiple-group comparisons. To make our approach accessible and clear to both noncomputational and computational biologists, we developed betAS, an interactive web app and user-friendly R package for visual and intuitive differential splicing analysis from read count data.

Replicative aging in yeast involves dynamic intron retention patterns associated with mRNA processing/export and protein ubiquitination

Saccharomyces cerevisiae (baker's yeast) has yielded relevant insights into some of the basic mechanisms of organismal aging. Among these are genomic instability, oxidative stress, caloric restriction and mitochondrial dysfunction. Several genes are known to have an impact on the aging process, with corresponding mutants exhibiting short- or long-lived phenotypes. Research dedicated to unraveling the underlying cellular mechanisms can support the identification of conserved mechanisms of aging in other species. One of the hitherto less studied fields in yeast aging is how the organism regulates its gene expression at the transcriptional level. To our knowledge, we present the first investigation into alternative splicing, particularly intron retention, during replicative aging of S. cerevisiae. This was achieved by utilizing the IRFinder algorithm on a previously published RNA-seq data set by Janssens et al. (2015). In the present work, 44 differentially retained introns in 43 genes were identified during replicative aging. We found that genes with altered intron retention do not display significant changes in overall transcript levels. It was possible to functionally assign distinct groups of these genes to the cellular processes of mRNA processing and export (e.g., YRA1) in early and middle-aged yeast, and protein ubiquitination (e.g., UBC5) in older cells. In summary, our work uncovers a previously unexplored layer of the transcriptional program of yeast aging and, more generally, expands the knowledge on the occurrence of alternative splicing in baker's yeast.

SETDB1 promotes progression through upregulation of SF3B4 expression and regulates the immunity in ovarian cancer

BACKGROUND

Ovarian cancer (OC) is the most lethal gynecologic malignant tumour. The mechanism promoting OC initiation and progression remains unclear. SET domain bifurcated histone lysine methyltransferase 1(SETDB1) acts as an oncogene in a variety of tumours. This study aims to explore the role of SETDB1 in OC.

METHODS

GEO, TCGA, CSIOVDB and CPTAC databases jointly analysed SETDB1 mRNA and protein expression. Effect of SETDB1 expression on the clinical prognosis of OC patients was analysed through online Kaplan‒Meier plotter and CSIOVDB database. Then, the effect of SETDB1 in OC cells progression and mobility was examined using MTT, EdU, colony formation and transwell assay. Additionally, Cistrome DB database was used to visualize the binding of SETDB1 protein and splicing factor 3b subunit 4 (SF3B4) promoter, and dual-luciferase reporter gene assay was performed to confirm the interaction. Finally, bioinformatics analysis was employed to reveal the relationship between SETDB1 and the microenvironment of OC.

RESULTS

In the present study, we found that SETDB1 was obviously upregulated in OC and its overexpression predicted poor prognosis of OC patients. Then, we verified that SETDB1 promoted the progression and motility of OC cells in vitro. Knockdown of SETDB1 had the opposite effect. Further research showed that SETDB1 acted as a transcription factor to activate SF3B4 expression. SF3B4 knockdown impaired the effect of SETDB1 to promote the proliferative capacity and motility of OC cells. Finally, the results of bioinformatics analysis confirmed that SETDB1 regulated the immune microenvironment of ovarian cancer.

CONCLUSION

SETDB1 promoted ovarian cancer progression by upregulating the expression of SF3B4 and inhibiting the tumour immunity. SETDB1 may be a promising prognostic and therapeutic marker for OC.

TRIM25 promotes glioblastoma cell growth and invasion via regulation of the PRMT1/c-MYC pathway by targeting the splicing factor NONO

BACKGROUND

Ubiquitination plays an important role in proliferating and invasive characteristic of glioblastoma (GBM), similar to many other cancers. Tripartite motif 25 (TRIM25) is a member of the TRIM family of proteins, which are involved in tumorigenesis through substrate ubiquitination.

METHODS

Difference in TRIM25 expression levels between nonneoplastic brain tissue samples and primary glioma samples was demonstrated using publicly available glioblastoma database, immunohistochemistry, and western blotting. TRIM25 knockdown GBM cell lines (LN229 and U251) and patient derived GBM stem-like cells (GSCs) GBM#021 were used to investigate the function of TRIM25 in vivo and in vitro. Co-immunoprecipitation (Co-IP) and mass spectrometry analysis were performed to identify NONO as a protein that interacts with TRIM25. The molecular mechanisms underlying the promotion of GBM development by TRIM25 through NONO were investigated by RNA-seq and validated by qRT-PCR and western blotting.

RESULTS

We observed upregulation of TRIM25 in GBM, correlating with enhanced glioblastoma cell growth and invasion, both in vitro and in vivo. Subsequently, we screened a panel of proteins interacting with TRIM25; mass spectrometry and co-immunoprecipitation revealed that NONO was a potential substrate of TRIM25. TRIM25 knockdown reduced the K63-linked ubiquitination of NONO, thereby suppressing the splicing function of NONO. Dysfunctional NONO resulted in the retention of the second intron in the pre-mRNA of PRMT1, inhibiting the activation of the PRMT1/c-MYC pathway.

CONCLUSIONS

Our study demonstrates that TRIM25 promotes glioblastoma cell growth and invasion by regulating the PRMT1/c-MYC pathway through mediation of the splicing factor NONO. Targeting the E3 ligase activity of TRIM25 or the complex interactions between TRIM25 and NONO may prove beneficial in the treatment of GBM.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity

BACKGROUND

Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms.

RESULTS

In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms.

CONCLUSION

This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

Transcriptome reprogramming through alternative splicing triggered by apigenin drives cell death in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is characterized by its aggressiveness and resistance to cancer-specific transcriptome alterations. Alternative splicing (AS) is a major contributor to the diversification of cancer-specific transcriptomes. The TNBC transcriptome landscape is characterized by aberrantly spliced isoforms that promote tumor growth and resistance, underscoring the need to identify approaches that reprogram AS circuitry towards transcriptomes, favoring a delay in tumorigenesis or responsiveness to therapy. We have previously shown that flavonoid apigenin is associated with splicing factors, including heterogeneous nuclear ribonucleoprotein A2 (hnRNPA2). Here, we showed that apigenin reprograms TNBC-associated AS transcriptome-wide. The AS events affected by apigenin were statistically enriched in hnRNPA2 substrates. Comparative transcriptomic analyses of human TNBC tumors and non-tumor tissues showed that apigenin can switch cancer-associated alternative spliced isoforms (ASI) to those found in non-tumor tissues. Apigenin preferentially affects the splicing of anti-apoptotic and proliferation factors, which are uniquely observed in cancer cells, but not in non-tumor cells. Apigenin switches cancer-associated aberrant ASI in vivo in TNBC xenograft mice by diminishing proliferation and increasing pro-apoptotic ASI. In accordance with these findings, apigenin increased apoptosis and reduced tumor proliferation, thereby halting TNBC growth in vivo. Our results revealed that apigenin reprograms transcriptome-wide TNBC-specific AS, thereby inducing apoptosis and hindering tumor growth. These findings underscore the impactful effects of nutraceuticals in altering cancer transcriptomes, offering new options to influence outcomes in TNBC treatments.

Unleashing the potential of catalytic RNAs to combat mis-spliced transcripts

Human transcriptome can undergo RNA mis-splicing due to spliceopathies contributing to the increasing number of genetic diseases including muscular dystrophy (MD), Alzheimer disease (AD), Huntington disease (HD), myelodysplastic syndromes (MDS). Intron retention (IR) is a major inducer of spliceopathies where two or more introns remain in the final mature mRNA and account for many intronic expansion diseases. Potential removal of such introns for therapeutic purposes can be feasible when utilizing bioinformatics, catalytic RNAs, and nano-drug delivery systems. Overcoming delivery challenges of catalytic RNAs was discussed in this review as a future perspective highlighting the significance of utilizing synthetic biology in addition to high throughput deep sequencing and computational approaches for the treatment of mis-spliced transcripts.

Tumor-wide RNA splicing aberrations generate immunogenic public neoantigens

T-cell-mediated immunotherapies are limited by the extent to which cancer-specific antigens are homogenously expressed throughout a tumor. We reasoned that recurrent splicing aberrations in cancer represent a potential source of tumor-wide and public neoantigens, and to test this possibility, we developed a novel pipeline for identifying neojunctions expressed uniformly within a tumor across diverse cancer types. Our analyses revealed multiple neojunctions that recur across patients and either exhibited intratumor heterogeneity or, in some cases, were tumor-wide. We identified CD8+ T-cell clones specific for neoantigens derived from tumor-wide and conserved neojunctions in GNAS and RPL22 , respectively. TCR-engineered CD8 + T-cells targeting these mutations conferred neoantigen-specific tumor cell eradication. Furthermore, we revealed that cancer-specific dysregulation in splicing factor expression leads to recurrent neojunction expression. Together, these data reveal that a subset of neojunctions are both intratumorally conserved and public, providing the molecular basis for novel T-cell-based immunotherapies that address intratumoral heterogeneity.

Functions of SRPK, CLK and DYRK kinases in stem cells, development, and human developmental disorders

Human developmental disorders encompass a wide range of debilitating physical conditions and intellectual disabilities. Perturbation of protein kinase signalling underlies the development of some of these disorders. For example, disrupted SRPK signalling is associated with intellectual disabilities, and the gene dosage of DYRKs can dictate the pathology of disorders including Down's syndrome. Here, we review the emerging roles of the CMGC kinase families SRPK, CLK, DYRK, and sub-family HIPK during embryonic development and in developmental disorders. In particular, SRPK, CLK, and DYRK kinase families have key roles in developmental signalling and stem cell regulation, and can co-ordinate neuronal development and function. Genetic studies in model organisms reveal critical phenotypes including embryonic lethality, sterility, musculoskeletal errors, and most notably, altered neurological behaviours arising from defects of the neuroectoderm and altered neuronal signalling. Further unpicking the mechanisms of specific kinases using human stem cell models of neuronal differentiation and function will improve our understanding of human developmental disorders and may provide avenues for therapeutic strategies.

Targeted splicing therapy: new strategies for colorectal cancer

RNA splicing is the process of forming mature mRNA, which is an essential phase necessary for gene expression and controls many aspects of cell proliferation, survival, and differentiation. Abnormal gene-splicing events are closely related to the development of tumors, and the generation of oncogenic isoform in splicing can promote tumor progression. As a main process of tumor-specific splicing variants, alternative splicing (AS) can promote tumor progression by increasing the production of oncogenic splicing isoforms and/or reducing the production of normal splicing isoforms. This is the focus of current research on the regulation of aberrant tumor splicing. So far, AS has been found to be associated with various aspects of tumor biology, including cell proliferation and invasion, resistance to apoptosis, and sensitivity to different chemotherapeutic drugs. This article will review the abnormal splicing events in colorectal cancer (CRC), especially the tumor-associated splicing variants arising from AS, aiming to offer an insight into CRC-targeted splicing therapy.

DHX38 restricts chemoresistance by regulating the alternative pre-mRNA splicing of RELL2 in pancreatic ductal adenocarcinoma

Intron retention plays an important role in cancer progression and chemotherapy resistance and seems to be essential for the maintenance of genome stability in cancer. Here, our goal was to analyze the role of receptor expressed in lymphoid tissue (Relt)-like 2 (RELL2) intron 4 retention in promoting pancreatic ductal adenocarcinoma (PDAC) progression. Our results showed that intron retention (IR) occurs at the fourth intron of RELL2 transcript in gemcitabine resistant PDAC cells, however, the regulatory mechanism and the clinical implications of IR of RELL2 are unclear. Firstly, we found that RELL2 plays an anti-oncogenic role in PDAC by performing in vitro functional assays including cell proliferation, GEM cytotoxicity assay and apoptosis. Subsequently, we identified the upstream gene of RELL2, DEAH-Box Helicase 38 (DHX38), and demonstrated the direct interaction between DHX38 and RELL2 by RIP-qPCR. We also found that altered expression of DHX38 resulted in corresponding changes in intron 4 retention of RELL2. Importantly, we unveiled that overexpression of DHX38 on the basis of knocking down of the fourth intron of RELL2 resulted in an impaired intron 4 intention. Overall, our study identified a new IR site in PDAC, which could be a possible target for PDAC therapy.

Cytoplasmic aggregation of mutant FUS causes multistep RNA splicing perturbations in the course of motor neuron pathology

Dysfunction of the RNA-binding protein (RBP) FUS implicated in RNA metabolism can cause amyotrophic lateral sclerosis (ALS) and other neurodegenerative diseases. Mutations affecting FUS nuclear localization can drive RNA splicing defects and stimulate the formation of non-amyloid inclusions in affected neurons. However, the mechanism by which FUS mutations contribute to the development of ALS remains uncertain. Here we describe a pattern of RNA splicing changes in the dynamics of the continuous proteinopathy induced by mislocalized FUS. We show that the decrease in intron retention of FUS-associated transcripts represents the hallmark of the pathogenesis of ALS and is the earliest molecular event in the course of progression of the disease. As FUS aggregation increases, the pattern of RNA splicing changes, becoming more complex, including a decrease in the inclusion of neuron-specific microexons and induction of cryptic exon splicing due to the sequestration of additional RBPs into FUS aggregates. Crucially, the identified features of the pathological splicing pattern are also observed in ALS patients in both sporadic and familial cases. Our data provide evidence that both a loss of nuclear FUS function due to mislocalization and the subsequent cytoplasmic aggregation of mutant protein lead to the disruption of RNA splicing in a multistep fashion during FUS aggregation.

Nuclear OsFKBP20-1b maintains SR34 stability and promotes the splicing of retained introns upon ABA exposure in rice

Alternative splicing (AS) is a critical means by which plants respond to changes in the environment, but few splicing factors contributing to AS have been reported and functionally characterized in rice (Oryza sativa L.). Here, we explored the function and molecular mechanism of the spliceosome-associated protein OsFKBP20-1b during AS. We determined the AS landscape of wild-type and osfkbp20-1b knockout plants upon abscisic acid (ABA) treatment by transcriptome deep sequencing. To capture the dynamics of translating intron-containing mRNAs, we blocked transcription with cordycepin and performed polysome profiling. We also analyzed whether OsFKBP20-1b and the splicing factors OsSR34 and OsSR45 function together in AS using protoplast transfection assays. We show that OsFKBP20-1b interacts with OsSR34 and regulates its stability, suggesting a role as a chaperone-like protein in the spliceosome. OsFKBP20-1b facilitates the splicing of mRNAs with retained introns after ABA treatment; some of these mRNAs are translatable and encode functional transcriptional regulators of stress-responsive genes. In addition, interacting proteins, OsSR34 and OsSR45, regulate the splicing of the same retained introns as OsFKBP20-1b after ABA treatment. Our findings reveal that spliceosome-associated immunophilin functions in alternative RNA splicing in rice by positively regulating the splicing of retained introns to limit ABA response.

Aberrant DNA Methylation, Expression, and Occurrence of Transcript Variants of the ABC Transporter ABCA7 in Breast Cancer

The ABC transporter ABCA7 has been found to be aberrantly expressed in a variety of cancer types, including breast cancer. We searched for specific epigenetic and genetic alterations and alternative splicing variants of ABCA7 in breast cancer and investigated whether these alterations are associated with ABCA7 expression. By analyzing tumor tissues from breast cancer patients, we found CpGs at the exon 5-intron 5 boundary aberrantly methylated in a molecular subtype-specific manner. The detection of altered DNA methylation in tumor-adjacent tissues suggests epigenetic field cancerization. In breast cancer cell lines, DNA methylation levels of CpGs in promoter-exon 1, intron 1, and at the exon 5-intron 5 boundary were not correlated with ABCA7 mRNA levels. By qPCR involving intron-specific and intron-flanking primers, we identified intron-containing ABCA7 mRNA transcripts. The occurrence of intron-containing transcripts was neither molecular subtype-specific nor directly correlated with DNA methylation at the respective exon-intron boundaries. Treatment of breast cancer cell lines MCF-7, BT-474, SK-BR3, and MDA-MB-231 with doxorubicin or paclitaxel for 72 h resulted in altered ABCA7 intron levels. Shotgun proteomics revealed that an increase in intron-containing transcripts was associated with significant dysregulation of splicing factors linked to alternative splicing.

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.

Noncoding translation mitigation

Translation is pervasive outside of canonical coding regions, occurring in long noncoding RNAs, canonical untranslated regions and introns1-4, especially in ageing4-6, neurodegeneration5,7 and cancer8-10. Notably, the majority of tumour-specific antigens are results of noncoding translation11-13. Although the resulting polypeptides are often nonfunctional, translation of noncoding regions is nonetheless necessary for the birth of new coding sequences14,15. The mechanisms underlying the surveillance of translation in diverse noncoding regions and how escaped polypeptides evolve new functions remain unclear10,16-19. Functional polypeptides derived from annotated noncoding sequences often localize to membranes20,21. Here we integrate massively parallel analyses of more than 10,000 human genomic sequences and millions of random sequences with genome-wide CRISPR screens, accompanied by in-depth genetic and biochemical characterizations. Our results show that the intrinsic nucleotide bias in the noncoding genome and in the genetic code frequently results in polypeptides with a hydrophobic C-terminal tail, which is captured by the ribosome-associated BAG6 membrane protein triage complex for either proteasomal degradation or membrane targeting. By contrast, canonical proteins have evolved to deplete C-terminal hydrophobic residues. Our results reveal a fail-safe mechanism for the surveillance of unwanted translation from diverse noncoding regions and suggest a possible biochemical route for the preferential membrane localization of newly evolved proteins.

Alternative splicing: a new breakthrough for understanding tumorigenesis and potential clinical applications

BACKGROUND

Alternative splicing (AS) is a post-transcriptional process that produces transcript variants, thus leading to transcriptome complexity. Recently, the scope of AS studies has been greatly expanded toward clinical applications owing to the abundance of RNA sequencing data.

OBJECTIVE

This review consists of two parts. We first summarize bioinformatic resources that are useful for large-scale cancer-related AS studies. We then highlight the research efforts to utilize AS events for predicting clinical outcomes and planning therapeutic strategies.

RESULTS

Computational approaches to interrogate AS events have been reviewed under three categories: (1) databases to provide functional and clinical annotation of AS events, (2) analytical tools to identify cancer-associated AS event, and (3) methods to identify splicing-related DNA variants and splicing-derived neoantigens. We also present the recent progress in exploring the clinical utility of AS under four categories: (1) identification of AS events for cancer prognosis, (2) utilization of AS events in molecular classification of various cancers, (3) regulatory mechanisms of AS underlying drug resistance, and (4) potential use of AS in cancer therapy.

CONCLUSION

This review will be helpful for understanding the biological implications of AS in cancer and facilitate the development of AS markers for cancer prognosis and treatment. We anticipate that future studies will lead to the application of genome-wide AS profiles in cancer precision medicine.

Abnormalities in intron retention characterize patients with systemic lupus erythematosus

Regulation of intron retention (IR), a form of alternative splicing, is a newly recognized checkpoint in gene expression. Since there are numerous abnormalities in gene expression in the prototypic autoimmune disease systemic lupus erythematosus (SLE), we sought to determine whether IR was intact in patients with this disease. We, therefore, studied global gene expression and IR patterns of lymphocytes in SLE patients. We analyzed RNA-seq data from peripheral blood T cell samples from 14 patients suffering from systemic lupus erythematosus (SLE) and 4 healthy controls and a second, independent data set of RNA-seq data from B cells from16 SLE patients and 4 healthy controls. We identified intron retention levels from 26,372 well annotated genes as well as differential gene expression and tested for differences between cases and controls using unbiased hierarchical clustering and principal component analysis. We followed with gene-disease enrichment analysis and gene-ontology enrichment analysis. Finally, we then tested for significant differences in intron retention between cases and controls both globally and with respect to specific genes. Overall decreased IR was found in T cells from one cohort and B cells from another cohort of patients with SLE and was associated with increased expression of numerous genes, including those encoding spliceosome components. Different introns within the same gene displayed both up- and down-regulated retention profiles indicating a complex regulatory mechanism. These results indicate that decreased IR in immune cells is characteristic of patients with active SLE and may contribute to the abnormal expression of specific genes in this autoimmune disease.

m6A-driven SF3B1 translation control steers splicing to direct genome integrity and leukemogenesis

SF3B1 is the most mutated splicing factor (SF) in myelodysplastic syndromes (MDSs), which are clonal hematopoietic disorders with variable risk of leukemic transformation. Although tumorigenic SF3B1 mutations have been extensively characterized, the role of "non-mutated" wild-type SF3B1 in cancer remains largely unresolved. Here, we identify a conserved epitranscriptomic program that steers SF3B1 levels to counteract leukemogenesis. Our analysis of human and murine pre-leukemic MDS cells reveals dynamic regulation of SF3B1 protein abundance, which affects MDS-to-leukemia progression in vivo. Mechanistically, ALKBH5-driven 5' UTR m⁶A demethylation fine-tunes SF3B1 translation directing splicing of central DNA repair and epigenetic regulators during transformation. This impacts genome stability and leukemia progression in vivo, supporting an integrative analysis in humans that SF3B1 molecular signatures may predict mutational variability and poor prognosis. These findings highlight a post-transcriptional gene expression nexus that unveils unanticipated SF3B1-dependent cancer vulnerabilities.

RNA splicing dysregulation and the hallmarks of cancer

Dysregulated RNA splicing is a molecular feature that characterizes almost all tumour types. Cancer-associated splicing alterations arise from both recurrent mutations and altered expression of trans-acting factors governing splicing catalysis and regulation. Cancer-associated splicing dysregulation can promote tumorigenesis via diverse mechanisms, contributing to increased cell proliferation, decreased apoptosis, enhanced migration and metastatic potential, resistance to chemotherapy and evasion of immune surveillance. Recent studies have identified specific cancer-associated isoforms that play critical roles in cancer cell transformation and growth and demonstrated the therapeutic benefits of correcting or otherwise antagonizing such cancer-associated mRNA isoforms. Clinical-grade small molecules that modulate or inhibit RNA splicing have similarly been developed as promising anticancer therapeutics. Here, we review splicing alterations characteristic of cancer cell transcriptomes, dysregulated splicing's contributions to tumour initiation and progression, and existing and emerging approaches for targeting splicing for cancer therapy. Finally, we discuss the outstanding questions and challenges that must be addressed to translate these findings into the clinic.

The eukaryotic translation initiation factor eIF4E reprograms alternative splicing

Aberrant splicing is typically attributed to splice-factor (SF) mutation and contributes to malignancies including acute myeloid leukemia (AML). Here, we discovered a mutation-independent means to extensively reprogram alternative splicing (AS). We showed that the dysregulated expression of eukaryotic translation initiation factor eIF4E elevated selective splice-factor production, thereby impacting multiple spliceosome complexes, including factors mutated in AML such as SF3B1 and U2AF1. These changes generated a splicing landscape that predominantly supported altered splice-site selection for ~800 transcripts in cell lines and ~4,600 transcripts in specimens from high-eIF4E AML patients otherwise harboring no known SF mutations. Nuclear RNA immunoprecipitations, export assays, polysome analyses, and mutational studies together revealed that eIF4E primarily increased SF production via its nuclear RNA export activity. By contrast, eIF4E dysregulation did not induce known SF mutations or alter spliceosome number. eIF4E interacted with the spliceosome and some pre-mRNAs, suggesting its direct involvement in specific splicing events. eIF4E induced simultaneous effects on numerous SF proteins, resulting in a much larger range of splicing alterations than in the case of mutation or dysregulation of individual SFs and providing a novel paradigm for splicing control and dysregulation.

Noncanonical splicing junctions between exons and transposable elements represent a source of immunogenic recurrent neo-antigens in patients with lung cancer

Although most characterized tumor antigens are encoded by canonical transcripts (such as differentiation or tumor-testis antigens) or mutations (both driver and passenger mutations), recent results have shown that noncanonical transcripts including long noncoding RNAs and transposable elements (TEs) can also encode tumor-specific neo-antigens. Here, we investigate the presentation and immunogenicity of tumor antigens derived from noncanonical mRNA splicing events between coding exons and TEs. Comparing human non-small cell lung cancer (NSCLC) and diverse healthy tissues, we identified a subset of splicing junctions that is both tumor specific and shared across patients. We used HLA-I peptidomics to identify peptides encoded by tumor-specific junctions in primary NSCLC samples and lung tumor cell lines. Recurrent junction-encoded peptides were immunogenic in vitro, and CD8+ T cells specific for junction-encoded epitopes were present in tumors and tumor-draining lymph nodes from patients with NSCLC. We conclude that noncanonical splicing junctions between exons and TEs represent a source of recurrent, immunogenic tumor-specific antigens in patients with NSCLC.

Non-canonical functions of spliceosome components in cancer progression

Dysregulation of pre-mRNA splicing is a common hallmark of cancer cells and it is associated with altered expression, localization, and mutations of the components of the splicing machinery. In the last few years, it has been elucidated that spliceosome components can also influence cellular processes in a splicing-independent manner. Here, we analyze open source data to understand the effect of the knockdown of splicing factors in human cells on the expression and splicing of genes relevant to cell proliferation, migration, cell cycle regulation, DNA repair, and cell death. We supplement this information with a comprehensive literature review of non-canonical functions of splicing factors linked to cancer progression. We also specifically discuss the involvement of splicing factors in intercellular communication and known autoregulatory mechanisms in restoring their levels in cells. Finally, we discuss strategies to target components of the spliceosome machinery that are promising for anticancer therapy. Altogether, this review greatly expands understanding of the role of spliceosome proteins in cancer progression.

MS0621, a novel small-molecule modulator of Ewing sarcoma chromatin accessibility, interacts with an RNA-associated macromolecular complex and influences RNA splicing

Ewing sarcoma is a cancer of children and young adults characterized by the critical translocation-associated fusion oncoprotein EWSR1::FLI1. EWSR1::FLI1 targets characteristic genetic loci where it mediates aberrant chromatin and the establishment of de novo enhancers. Ewing sarcoma thus provides a model to interrogate mechanisms underlying chromatin dysregulation in tumorigenesis. Previously, we developed a high-throughput chromatin-based screening platform based on the de novo enhancers and demonstrated its utility in identifying small molecules capable of altering chromatin accessibility. Here, we report the identification of MS0621, a molecule with previously uncharacterized mechanism of action, as a small molecule modulator of chromatin state at sites of aberrant chromatin accessibility at EWSR1::FLI1-bound loci. MS0621 suppresses cellular proliferation of Ewing sarcoma cell lines by cell cycle arrest. Proteomic studies demonstrate that MS0621 associates with EWSR1::FLI1, RNA binding and splicing proteins, as well as chromatin regulatory proteins. Surprisingly, interactions with chromatin and many RNA-binding proteins, including EWSR1::FLI1 and its known interactors, were RNA-independent. Our findings suggest that MS0621 affects EWSR1::FLI1-mediated chromatin activity by interacting with and altering the activity of RNA splicing machinery and chromatin modulating factors. Genetic modulation of these proteins similarly inhibits proliferation and alters chromatin in Ewing sarcoma cells. The use of an oncogene-associated chromatin signature as a target allows for a direct approach to screen for unrecognized modulators of epigenetic machinery and provides a framework for using chromatin-based assays for future therapeutic discovery efforts.

Transcripts of the Prostate Cancer-Associated Gene ANO7 Are Retained in the Nuclei of Prostatic Epithelial Cells

Prostate cancer affects millions of men globally. The prostate cancer-associated gene ANO7 is downregulated in advanced prostate cancer, whereas benign tissue and low-grade cancer display varying expression levels. In this study, we assess the spatial correlation between ANO7 mRNA and protein using fluorescent in situ hybridization and immunohistochemistry for the detection of mRNA and protein in parallel sections of tissue microarrays prepared from radical prostatectomy samples. We show that ANO7 mRNA and protein expression correlate in prostate tissue. Furthermore, we show that ANO7 mRNA is enriched in the nuclei of the luminal cells at 89% in benign ducts and low-grade cancer, and at 78% in high-grade cancer. The nuclear enrichment of ANO7 mRNA was validated in prostate cancer cell lines 22Rv1 and MDA PCa 2b using droplet digital polymerase chain reaction (ddPCR) on RNA isolated from nuclear and cytoplasmic fractions of the cells. The nuclear enrichment of ANO7 mRNA was compared to the nuclearly-enriched lncRNA MALAT1, confirming the surprisingly high nuclear retention of ANO7 mRNA. ANO7 has been suggested to be used as a diagnostic marker and a target for immunotherapy, but a full comprehension of its role in prostate cancer progression is currently lacking. Our results contribute to a better understanding of the dynamics of ANO7 expression in prostatic tissue.

Endogenous retroelements as alarms for disruptions to cellular homeostasis

Endogenous retroelements are DNA sequences which can duplicate and move to new locations in the genome. Actively moving endogenous retroelements can be disruptive to the host, and their expression is therefore often repressed. Interestingly, drugs that disrupt the repression of endogenous retroelements show promise for treating cancer. Expressed endogenous retroelements can activate innate immune receptors that activate the antiviral response, potentially leading to the death of cancer cells. We discuss disruptions to cellular processes which can lead to activation of the antiviral state from endogenous retroelements, and present the 'fire alarm hypothesis', where we argue that endogenous retroelements act as alarms for disruptions to these cellular processes. Furthermore, we discuss the properties of endogenous retroelements which make them suitable as alarms.

Exploration of prognostic genes and risk signature in breast cancer patients based on RNA binding proteins associated with ferroptosis

Background: Breast cancer (BRCA) is a life-threatening malignancy in women with an unsatisfactory prognosis. The purpose of this study was to explore the prognostic biomarkers and a risk signature based on ferroptosis-related RNA-binding proteins (FR-RBPs). Methods: FR-RBPs were identified using Spearman correlation analysis. Differentially expressed genes (DEGs) were identified by the "limma" R package. The univariate Cox and multivariate Cox analyses were executed to determine the prognostic genes. The risk signature was constructed and verified with the training set, testing set, and validation set. Mutation analysis, immune checkpoint expression analysis in high- and low-risk groups, and correlation between risk signature and chemotherapeutic agents were conducted using the "maftools" package, "ggplot2" package, and the CellMiner database respectively. The Human Protein Atlas (HPA) database was employed to confirm protein expression trends of prognostic genes in BRCA and normal tissues. The expression of prognostic genes in cell lines was verified by Real-time quantitative polymerase chain reaction (RT-qPCR). Kaplan-meier (KM) plotter database analysis was applied to predict the correlation between the expression levels of signature genes and survival statuses. Results: Five prognostic genes (GSPT2, RNASE1, TIPARP, TSEN54, and SAMD4A) to construct an FR-RBPs-related risk signature were identified and the risk signature was validated by the International Cancer Genome Consortium (ICGC) cohort. Univariate and multivariate Cox regression analysis demonstrated the risk score was a robust independent prognostic factor in overall survival prediction. The Tumor Mutational Burden (TMB) analysis implied that the high- and low-risk groups responded differently to immunotherapy. Drug sensitivity analysis suggested that the risk signature may serve as a chemosensitivity predictor. The results of GSEA suggested that five prognostic genes might be related to DNA replication and the immune-related pathways. RT-qPCR results demonstrated that the expression trends of prognostic genes in cell lines were consistent with the results from public databases. KM plotter database analysis suggested that high expression levels of GSPT2, RNASE1, and SAMD4A contributed to poor prognoses. Conclusion: In conclusion, this study identified the FR-RBPs-related prognostic genes and developed an FR-RBPs-related risk signature for the prognosis of BRCA, which will be of great significance in developing new therapeutic targets and prognostic molecular biomarkers for BRCA.

Alternative splicing and genetic variation of mhc-e: implications for rhesus cytomegalovirus-based vaccines

Rhesus cytomegalovirus (RhCMV)-based vaccination against Simian Immunodeficiency virus (SIV) elicits MHC-E-restricted CD8+ T cells that stringently control SIV infection in ~55% of vaccinated rhesus macaques (RM). However, it is unclear how accurately the RM model reflects HLA-E immunobiology in humans. Using long-read sequencing, we identified 16 Mamu-E isoforms and all Mamu-E splicing junctions were detected among HLA-E isoforms in humans. We also obtained the complete Mamu-E genomic sequences covering the full coding regions of 59 RM from a RhCMV/SIV vaccine study. The Mamu-E gene was duplicated in 32 (54%) of 59 RM. Among four groups of Mamu-E alleles: three ~5% divergent full-length allele groups (G1, G2, G2_LTR) and a fourth monomorphic group (G3) with a deletion encompassing the canonical Mamu-E exon 6, the presence of G2_LTR alleles was significantly (p = 0.02) associated with the lack of RhCMV/SIV vaccine protection. These genomic resources will facilitate additional MHC-E targeted translational research.

RNA and RNA Derivatives: Light and Dark Sides in Cancer Immunotherapy

Significance: Immunotherapy, which utilizes the patient's immune system to fight tumor cells, has been approved for the treatment of some types of advanced cancer. Recent Advances: The complexity and diversity of tumor immunity are responsible for the varying response rates toward current immunotherapy strategies and highlight the importance of exploring regulators in tumor immunotherapy. Several genetic factors have proved to be critical regulators of tumor immunotherapy. RNAs, including messenger RNAs and non-coding RNAs, play vital and diverse roles in tumorigenesis, metastasis, drug resistance, and immunotherapy response. RNA modifications, including N6-methyladenosine methylation, are involved in tumor immunity. Critical Issues: A critical issue is the lack of summary of the regulatory RNA molecules and their derivatives in mediating immune activities in human cancers that could provide potential applications for tumor immunotherapeutic strategy. Future Directions: This review summarizes the dual roles (the light and dark sides) of RNA and its derivatives in tumor immunotherapy and discusses the development of RNA-based therapies as novel immunotherapeutic strategies for cancer treatment. Antioxid. Redox Signal. 37, 1266-1290.

Retained introns in long RNA-seq reads are not reliably detected in sample-matched short reads

BACKGROUND

There is growing interest in retained introns in a variety of disease contexts including cancer and aging. Many software tools have been developed to detect retained introns from short RNA-seq reads, but reliable detection is complicated by overlapping genes and transcripts as well as the presence of unprocessed or partially processed RNAs.

RESULTS

We compared introns detected by 8 tools using short RNA-seq reads with introns observed in long RNA-seq reads from the same biological specimens. We found significant disagreement among tools (Fleiss' [Formula: see text]) such that 47.7% of all detected intron retentions were not called by more than one tool. We also observed poor performance of all tools, with none achieving an F1-score greater than 0.26, and qualitatively different behaviors between general-purpose alternative splicing detection tools and tools confined to retained intron detection.

CONCLUSIONS

Short-read tools detect intron retention with poor recall and precision, calling into question the completeness and validity of a large percentage of putatively retained introns called by commonly used methods.

The role of NPM1 alternative splicing in patients with chronic lymphocytic leukemia

OBJECTIVES

Chronic lymphocytic leukemia (CLL) is a lymphoproliferative disease with heterogeneous clinical course. Recent studies revealed a link between NOTCH1 mutation and the overexpression of MYC and MYC-related genes involved in ribosome biogenesis and protein biosynthesis, such as nucleophosmin-1 (NPM1), in CLL cells. In the present study, we aim to evaluate the impact of the NOTCH1 mutation on the MYC and MYC induced NPM1 expression in CLL cells via quantification of their transcripts.

METHODS

Using qRT-PCR, we analyzed the levels of MYC and three main NPM1 splice variants in 214 samples collected from CLL patients. We assessed the impact of each splice variant on CLL prognostic markers, including the IGHV, TP53, NOTCH1, SF3B1, and MYD88 mutational status, cytogenetic aberrations, and laboratory features.

RESULTS

Significantly higher levels of NPM1.R1 transcripts in patients with unmutated compared to mutated IGHV status were found. The median time to first treatment (TTFT) in patients with a high level of NPM1.R1 was significantly shorter compared to the group with low NPM1.R1 levels (1.5 vs 33 months, p = 0.0002). Moreover, in Multivariate Cox Proportional Hazard Regression Model NPM1.R1 splice variant provided an independent prognostic value for TTFT.

CONCLUSION

In conclusion, our study indicates the prognostic significance of the level of NPM1.R1 expression and suggests the importance of splicing alterations in the pathogenesis of CLL.

Physiological intron retaining transcripts in the cytoplasm abound during human motor neurogenesis

Intron retention (IR) is now recognized as a dominant splicing event during motor neuron (MN) development; however, the role and regulation of intron-retaining transcripts (IRTs) localized to the cytoplasm remain particularly understudied. Here we show that IR is a physiological process that is spatiotemporally regulated during MN lineage restriction and that IRTs in the cytoplasm are detected in as many as 13% (n = 2297) of the genes expressed during this process. We identify a major class of cytoplasmic IRTs that are not associated with reduced expression of their own genes but instead show a high capacity for RNA-binding protein and miRNA occupancy. Finally, we show that ALS-causing VCP mutations lead to a selective increase in cytoplasmic abundance of this particular class of IRTs, which in turn temporally coincides with an increase in the nuclear expression level of predicted miRNA target genes. Altogether, our study identifies a previously unrecognized class of cytoplasmic intronic sequences with potential regulatory function beyond gene expression.

RNA velocity unraveled

We perform a thorough analysis of RNA velocity methods, with a view towards understanding the suitability of the various assumptions underlying popular implementations. In addition to providing a self-contained exposition of the underlying mathematics, we undertake simulations and perform controlled experiments on biological datasets to assess workflow sensitivity to parameter choices and underlying biology. Finally, we argue for a more rigorous approach to RNA velocity, and present a framework for Markovian analysis that points to directions for improvement and mitigation of current problems.

Single-cell sequencing data are snapshots of biological processes, making it challenging to infer dynamic relationships between cell types. RNA velocity attempts to bypass this challenge by treating the unspliced RNA content as a proxy for spliced RNA content in the near future, and using this “extrapolation” to build directional relationships. However, the method, as implemented in several software packages, is not yet reliable enough to be actionable, in part due to the large number of arbitrary, user-set hyperparameters, as well as fundamental incompatibilities between the biophysics of transcription in the living cell and the models used throughout the velocity workflows. In this study, we review these issues, and use existing results from the fields of stochastic modeling and fluorescence transcriptomics to develop an alternative theoretical framework. We show that our framework can facilitate the development and inference of physically consistent models for sequencing data, as well as the unification of single-cell analyses to self-consistently treat variation due to cell type dynamics and identities, the stochasticity inherent to single-molecule processes, and the uncertainty introduced by sequencing experiments.

Splicing dysregulation in human hematologic malignancies: beyond splicing mutations

Splicing is a fundamental process in pre-mRNA maturation. Whereas alternative splicing (AS) enriches the diversity of the proteome, its aberrant regulation can drive oncogenesis. So far, most attention has been given to spliceosome mutations (SMs) in the context of splicing dysregulation in hematologic diseases. However, in recent years, post-translational modifications (PTMs) and transcriptional alterations of splicing factors (SFs), just as epigenetic signatures, have all been shown to contribute to global splicing dysregulation as well. In addition, the contribution of aberrant splicing to the neoantigen repertoire of cancers has been recognized. With the pressing need for novel therapeutics to combat blood cancers, this article provides an overview of emerging mechanisms that contribute to aberrant splicing, as well as their clinical potential.

3’UTR-Seq analysis of chicken abdominal adipose tissue reveals widespread intron retention in 3’UTR and provides insight into molecular basis of feed efficiency

Feed efficiency (FE) is an important trait in the broiler industry due to its direct correlation to efficient muscle growth instead of fat deposition. The present study characterized and compared gene expression profiles in abdominal fat from broiler chickens of different FE levels to enhance the understanding of FE biology. Specifically, traditional whole-transcript RNA-sequencing (RNA-seq) and 3’ UTR-sequencing (3’ UTR-seq) were applied to 22 and 61 samples, respectively. Overall, these two sequencing techniques shared a high correlation (0.76) between normalized counts, although 3’ UTR-seq showed a higher variance in sequencing and mapping performance statistics across samples and a lower rate of uniquely mapped reads. A higher percentage of 3’ UTR-seq reads mapped to introns suggested the frequent presence of cleavage sites in introns, thus warranting future research to study its regulatory function. Differential expression analysis identified 1198 differentially expressed genes (DEGs) between high FE (HFE) and intermediate FE (IFE) chickens with False Discovery Rate < 0.05 and fold change > 1.2. The processes that were significantly enriched by the DEGs included extracellular matrix remodeling and mechanisms impacting gene expression at the transcriptional and translational levels. Gene ontology enrichment analysis suggested that the divergence in fat deposition and FE in broiler chickens could be associated with peroxisome and lipid metabolism possibly regulated by G0/G1 switch gene 2 (G0S2).

Selective ablation of 3' RNA ends and processive RTs facilitate direct cDNA sequencing of full-length host cell and viral transcripts

Alternative splicing (AS) is necessary for viral proliferation in host cells and a critical regulatory component of viral gene expression. Conventional RNA-seq approaches provide incomplete coverage of AS due to their short read lengths and are susceptible to biases and artifacts introduced in prevailing library preparation methodologies. Moreover, viral splicing studies are often conducted separately from host cell transcriptome analysis, precluding an assessment of the viral manipulation of host splicing machinery. To address current limitations, we developed a quantitative full-length direct cDNA sequencing strategy to simultaneously profile viral and host cell transcripts. This nanopore-based approach couples processive reverse transcriptases with a novel one-step chemical ablation of 3' RNA ends (termed CASPR), which decreases ribosomal RNA reads and enriches polyadenylated coding sequences. We extensively validate our approach using synthetic reference transcripts and show that CASPR doubles the breadth of coverage per transcript and increases detection of long transcripts (>4 kb), while being functionally equivalent to PolyA+ selection for transcript quantification. We used our approach to interrogate host cell and HIV-1 transcript dynamics during viral reactivation and identified novel putative HIV-1 host factors containing exon skipping or novel intron retentions and delineated the HIV-1 transcriptional state associated with these differentially regulated host factors.

Alternative Splicing of the Aryl Hydrocarbon Receptor Nuclear Translocator (ARNT) Is Regulated by RBFOX2 in Lymphoid Malignancies

Aberrant alternative splicing (AS) of pre-mRNAs promotes the development and proliferation of cancerous cells. Accordingly, we had previously observed higher levels of the aryl hydrocarbon receptor nuclear translocator (ARNT) spliced variant isoform 1 in human lymphoid malignancies compared to that in normal lymphoid cells, which is a consequence of increased inclusion of alternative exon 5. ARNT is a transcription factor that has been implicated in the survival of various cancers. Notably, we found that ARNT isoform 1 promoted the growth and survival of lymphoid malignancies, but the regulatory mechanism controlling ARNT AS is unclear. Here, we report cis- and trans-regulatory elements which are important for the inclusion of ARNT exon 5. Specifically, we identified recognition motifs for the RNA-binding protein RBFOX2, which are required for RBFOX2-mediated exon 5 inclusion. RBFOX2 upregulation was observed in lymphoid malignancies, correlating with the observed increase in ARNT exon 5 inclusion. Moreover, suppression of RBFOX2 significantly reduced ARNT isoform 1 levels and cell growth. These observations reveal RBFOX2 as a critical regulator of ARNT AS in lymphoid malignancies and suggest that blocking the ARNT-specific RBFOX2 motifs to decrease ARNT isoform 1 levels is a viable option for targeting the growth of lymphoid malignancies.

Can aggressive cancers be identified by the "aggressiveness" of their chromatin?

Phenotypic plasticity is a crucial feature of aggressive cancer, providing the means for cancer progression. Stochastic changes in tumor cell transcriptional programs increase the chances of survival under any condition. I hypothesize that unstable chromatin permits stochastic transitions between transcriptional programs in aggressive cancers and supports non-genetic heterogeneity of tumor cells as a basis for their adaptability. I present a mechanistic model for unstable chromatin which includes destabilized nucleosomes, mobile chromatin fibers and random enhancer-promoter contacts, resulting in stochastic transcription. I suggest potential markers for "unsettled" chromatin in tumors associated with poor prognosis. Although many of the characteristics of unstable chromatin have been described, they were mostly used to explain changes in the transcription of individual genes. I discuss approaches to evaluate the role of unstable chromatin in non-genetic tumor cell heterogeneity and suggest using the degree of chromatin instability and transcriptional noise in tumor cells to predict cancer aggressiveness.

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

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

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

BACKGROUND

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

OBJECTIVE

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

METHODS

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

RESULTS

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

DISCUSSION

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