CDK9 inhibition activates innate immune response through viral mimicry

Cancer cells frequently exhibit hyperactivation of transcription, which can lead to increased sensitivity to compounds targeting the transcriptional kinases, in particular CDK9. However, mechanistic details of CDK9 inhibition-induced cancer cell-selective anti-proliferative effects remain largely unknown. Here, we discover that CDK9 inhibition activates the innate immune response through viral mimicry in cancer cells. In MYC over-expressing prostate cancer cells, CDK9 inhibition leads to the gross accumulation of mis-spliced RNA. Double-stranded RNA (dsRNA)-activated kinase can recognize these mis-spliced RNAs, and we show that the activity of this kinase is required for the CDK9 inhibitor-induced anti-proliferative effects. Using time-resolved transcriptional profiling (SLAM-seq), targeted proteomics, and ChIP-seq, we show that, similar to viral infection, CDK9 inhibition significantly suppresses transcription of most genes but allows selective transcription and translation of cytokines related to the innate immune response. In particular, CDK9 inhibition activates NFκB-driven cytokine signaling at the transcriptional and secretome levels. The transcriptional signature induced by CDK9 inhibition identifies prostate cancers with a high level of genome instability. We propose that it is possible to induce similar effects in patients using CDK9 inhibition, which, we show, causes DNA damage in vitro. In the future, it is important to establish whether CDK9 inhibitors can potentiate the effects of immunotherapy against late-stage prostate cancer, a currently lethal disease.

The Arabidopsis SR45 splicing factor bridges the splicing machinery and the exon-exon junction complex

The Arabidopsis splicing factor serine/arginine-rich 45 (SR45) contributes to several biological processes. The sr45-1 loss-of-function mutant exhibits delayed root development, late flowering, unusual numbers of floral organs, shorter siliques with decreased seed sets, narrower leaves and petals, and altered metal distribution. SR45 bears a unique RNA recognition motif (RRM) flanked by one serine/arginine-rich (RS) domain on both sides. Here, we studied the function of each SR45 domains by examining their involvement in: (i) the spatial distribution of SR45; (ii) the establishment of a protein-protein interaction network including spliceosomal and exon-exon junction complex (EJC) components; and (iii) the RNA binding specificity. We report that the endogenous SR45 promoter is active during vegetative and reproductive growth, and that the SR45 protein localizes in the nucleus. We demonstrate that the C-terminal arginine/serine-rich domain is a determinant of nuclear localization. We show that the SR45 RRM domain specifically binds purine-rich RNA motifs via three residues (H101, H141, and Y143), and is also involved in protein-protein interactions. We further show that SR45 bridges both mRNA splicing and surveillance machineries as a partner of EJC core components and peripheral factors, which requires phosphoresidues probably phosphorylated by kinases from both the CLK and SRPK families. Our findings provide insights into the contribution of each SR45 domain to both spliceosome and EJC assemblies.

Splice-Switching Antisense Oligonucleotides Correct Phenylalanine Hydroxylase Exon 11 Skipping Defects and Rescue Enzyme Activity in Phenylketonuria

The PAH gene encodes the hepatic enzyme phenylalanine hydroxylase (PAH), and its deficiency, known as phenylketonuria (PKU), leads to neurotoxic high levels of phenylalanine. PAH exon 11 is weakly defined, and several missense and intronic variants identified in patients affect the splicing process. Recently, we identified a novel intron 11 splicing regulatory element where U1snRNP binds, participating in exon 11 definition. In this work, we describe the implementation of an antisense strategy targeting intron 11 sequences to correct the effect of PAH mis-splicing variants. We used an in vitro assay with minigenes and identified splice-switching antisense oligonucleotides (SSOs) that correct the exon skipping defect of PAH variants c.1199+17G>A, c.1199+20G>C, c.1144T>C, and c.1066-3C>T. To examine the functional rescue induced by the SSOs, we generated a hepatoma cell model with variant c.1199+17G>A using CRISPR/Cas9. The edited cell line reproduces the exon 11 skipping pattern observed from minigenes, leading to reduced PAH protein levels and activity. SSO transfection results in an increase in exon 11 inclusion and corrects PAH deficiency. Our results provide proof of concept of the potential therapeutic use of a single SSO for different exonic and intronic splicing variants causing PAH exon 11 skipping in PKU.

TERT promoter mutations increase tumor aggressiveness by altering TERT mRNA splicing in papillary thyroid carcinoma

CONTEXT

Telomerase reverse transcriptase promoter (TERT-p) mutations, which upregulate TERT expression, are strongly associated with tumor aggressiveness and worse prognosis in papillary thyroid carcinomas (PTCs). TERT expression is also observed in a proportion of PTCs without TERT-p mutations, but such tumors show less aggressiveness and better prognosis compared with TERT-p mutation-positive tumors.

OBJECTIVE

TERT has multiple splicing variants whose relationships with the TERT-p status and clinicopathological characteristics remain poorly understood. We examined the relationship between the TERT-p mutational status, the TERT splicing pattern, and clinicopathological features.

METHODS

We investigated the expression of two major variants, α deletion (dA) and β deletion (dB), in a series of 207 PTCs operated between November 2001 and March 2020 in Nagasaki University Hospital and Kuma Hospital.

RESULTS

The TERT-p mutations were found in 33 cases, and among 174 mutation-negative cases, 24 showed TERT expression. All cases were classified into three groups: the TERT-p mutation-negative/expression-negative group (mut-/exp-), the TERT-p mutation-negative/expression-positive group (mut-/exp+), and the TERT-p mutation-positive group (mut+/exp+). The +A + B/dB ratio in mut+/exp + was significantly higher than that in mut-/exp + PTCs. Analysis with clinicopathological data revealed that +A + B expression was associated with higher PTC aggressiveness, whereas dB expression counteracted this effect. Functional in vitro study demonstrated that dB strongly inhibited cell growth, migration, and clonogenicity, suggesting its tumor suppressive role.

CONCLUSION

These results provide evidence that the TERT-p mutations alter the expression of different TERT splice variants, which, in turn, associates with different tumor aggressiveness.

Post-transcriptional splicing can occur in a slow-moving zone around the gene

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

Deciphering the impact of coding and non-coding SCN1A gene variants on RNA splicing

Variants that disrupt normal pre-mRNA splicing are increasingly being recognized as a major cause of monogenic disorders. The SCN1A gene, a key epilepsy gene that is linked to various epilepsy phenotypes, is no exception. Approximately 10% of all reported variants in the SCN1A gene are designated as splicing variants, with many located outside of the canonical donor and acceptor splice sites, and most have not been functionally investigated. However, given its restricted expression pattern, functional analysis of splicing variants in the SCN1A gene could not be routinely performed. In this study, we conducted a comprehensive analysis of all reported SCN1A variants and their potential to impact SCN1A splicing and conclude that splicing variants are substantially misannotated and under-represented. We created a splicing reporter system consisting of 18 splicing vectors covering all 26 protein-coding exons with different genomic contexts and several promoters of varying strengths in order to reproduce the wild-type splicing pattern of the SCN1A gene, revealing cis-regulatory elements essential for proper recognition of SCN1A exons. Functional analysis of 95 SCN1A variants was carried out, including all 68 intronic variants reported in the literature, located outside of the splice sites canonical dinucleotides; 21 exonic variants of different classes (synonymous, missense, nonsense and in-frame deletion) and six variants observed in patients with epilepsy. Interestingly, almost 20% of tested intronic variants had no influence on SCN1A splicing, despite being reported as causative in the literature. Moreover, we confirmed that the majority of predicted exonic variants affect splicing unravelling their true molecular mechanism. We used functional data to perform genotype-phenotype correlation, revealing distinct distribution patterns for missense and splice-affecting 'missense' variants and observed no difference in the phenotype severity of variants leading to in-frame and out-of-frame isoforms, indicating that the Nav1.1 protein is highly intolerant to structural variations. Our work demonstrates the importance of functional analysis in proper variant annotation and provides a tool for high-throughput delineation of splice-affecting variants in SCN1A in a whole-gene manner.

Molecular and Cellular Characterization of Primary Endothelial Cells from a Familial Cavernomatosis Patient

Cerebral cavernous malformation (CCM) or familial cavernomatosis is a rare, autosomal dominant, inherited disease characterized by the presence of vascular malformations consisting of blood vessels with an abnormal structure in the form of clusters. Based on the altered gene (CCM1/Krit1, CCM2, CCM3) and its origin (spontaneous or familial), different types of this disease can be found. In this work we have isolated and cultivated primary endothelial cells (ECs) from peripheral blood of a type 1 CCM patient. Differential functional and gene expression profiles of these cells were analyzed and compared to primary ECs from a healthy donor. The mutation of the familial index case consisted of a heterozygous point mutation in the position +1 splicing consensus between exons 15 and 16, causing failure in RNA processing and in the final protein. Furthermore, gene expression analysis by quantitative PCR revealed a decreased expression of genes involved in intercellular junction formation, angiogenesis, and vascular homeostasis. Cell biology analysis showed that CCM1 ECs were impaired in angiogenesis and cell migration. Taken together, the results obtained suggest that the alterations found in CCM1 ECs are already present in the heterozygous condition, suffering from vascular impairment and somewhat predisposed to vascular damage.

Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

CIRI-Deep Enables Single-Cell and Spatial Transcriptomic Analysis of Circular RNAs with Deep Learning

Circular RNAs (circRNAs) are a crucial yet relatively unexplored class of transcripts known for their tissue- and cell-type-specific expression patterns. Despite the advances in single-cell and spatial transcriptomics, these technologies face difficulties in effectively profiling circRNAs due to inherent limitations in circRNA sequencing efficiency. To address this gap, a deep learning model, CIRI-deep, is presented for comprehensive prediction of circRNA regulation on diverse types of RNA-seq data. CIRI-deep is trained on an extensive dataset of 25 million high-confidence circRNA regulation events and achieved high performances on both test and leave-out data, ensuring its accuracy in inferring differential events from RNA-seq data. It is demonstrated that CIRI-deep and its adapted version enable various circRNA analyses, including cluster- or region-specific circRNA detection, BSJ ratio map visualization, and trans and cis feature importance evaluation. Collectively, CIRI-deep's adaptability extends to all major types of RNA-seq datasets including single-cell and spatial transcriptomic data, which will undoubtedly broaden the horizons of circRNA research.

Long-read direct RNA sequencing of the mitochondrial transcriptome of Saccharomyces cerevisiae reveals condition-dependent intron abundance

Mitochondria fulfil many essential roles and have their own genome, which is expressed as polycistronic transcripts that undergo co- or posttranscriptional processing and splicing. Due to the inherent complexity and limited technical accessibility of the mitochondrial transcriptome, fundamental questions regarding mitochondrial gene expression and splicing remain unresolved, even in the model eukaryote Saccharomyces cerevisiae. Long-read sequencing could address these fundamental questions. Therefore, a method for the enrichment of mitochondrial RNA and sequencing using Nanopore technology was developed, enabling the resolution of splicing of polycistronic genes and the quantification of spliced RNA. This method successfully captured the full mitochondrial transcriptome and resolved RNA splicing patterns with single-base resolution and was applied to explore the transcriptome of S. cerevisiae grown with glucose or ethanol as the sole carbon source, revealing the impact of growth conditions on mitochondrial RNA expression and splicing. This study uncovered a remarkable difference in the turnover of Group II introns between yeast grown in either mostly fermentative or fully respiratory conditions. Whether this accumulation of introns in glucose medium has an impact on mitochondrial functions remains to be explored. Combined with the high tractability of the model yeast S. cerevisiae, the developed method enables to monitor mitochondrial transcriptome responses in a broad range of relevant contexts, including oxidative stress, apoptosis and mitochondrial diseases.

Exitron splicing of odor receptor genes in Drosophila

Proper expression of odor receptor genes is critical for the function of olfactory systems. In this study, we identified exitrons (exonic introns) in four of the 39 Odorant receptor (Or) genes expressed in the Drosophila antenna. Exitrons are sequences that can be spliced out from within a protein-coding exon, thereby altering the encoded protein. We focused on Or88a, which encodes a pheromone receptor, and found that exitron splicing of Or88a is conserved across five Drosophila species over 20 My of evolution. The exitron was spliced out in 15% of Or88a transcripts. Removal of this exitron creates a non-coding RNA rather than an RNA that encodes a stable protein. Our results suggest the hypothesis that in the case of Or88a, exitron splicing could act in neuronal modulation by decreasing the level of functional Or transcripts. Activation of Or88a-expressing olfactory receptor neurons via either optogenetics or pheromone stimulation increased the level of exitron-spliced transcripts, with optogenetic activation leading to a 14-fold increase. A fifth Or can also undergo an alternative splicing event that eliminates most of the canonical open reading frame. Besides these cases of alternative splicing, we found alternative polyadenylation of four Ors, and exposure of Or67c to its ligand ethyl lactate in the antenna downregulated all of its 3' isoforms. Our study reveals mechanisms by which neuronal activity could be modulated via regulation of the levels of Or isoforms.

Kin17 regulates proper cortical localization of Miranda in Drosophila neuroblasts by regulating Flfl expression

During asymmetric division of Drosophila larval neuroblasts, the fate determinant Prospero (Pros) and its adaptor Miranda (Mira) are segregated to the basal cortex through atypical protein kinase C (aPKC) phosphorylation of Mira and displacement from the apical cortex, but Mira localization after aPKC phosphorylation is not well understood. We identify Kin17, a DNA replication and repair protein, as a regulator of Mira localization during asymmetric cell division. Loss of Kin17 leads to aberrant localization of Mira and Pros to the centrosome, cytoplasm, and nucleus. We provide evidence to show that the mislocalization of Mira and Pros is likely due to reduced expression of Falafel (Flfl), a component of protein phosphatase 4 (PP4), and defects in dephosphorylation of serine-96 of Mira. Our work reveals that Mira is likely dephosphorylated by PP4 at the centrosome to ensure proper basal localization of Mira after aPKC phosphorylation and that Kin17 regulates PP4 activity by regulating Flfl expression.

RNA analysis and computer-aided facial phenotyping help to classify a novel TRIO splice site variant

Pathogenic variants in TRIO, encoding the guanine nucleotide exchange factor, are associated with two distinct neurodevelopmental delay phenotypes: gain-of-function missense mutations within the spectrin repeats are causative for a severe developmental delay with macrocephaly (MIM: 618825), whereas loss-of-function missense variants in the GEF1 domain and truncating variants throughout the gene lead to a milder developmental delay and microcephaly (MIM: 617061). In three affected family members with mild intellectual disability/NDD and microcephaly, we detected a novel heterozygous TRIO variant at the last coding base of exon 31 (NM_007118.4:c.4716G>A). RNA analysis from patient-derived lymphoblastoid cells confirmed aberrant splicing resulting in the skipping of exon 31 (r.4615_4716del), leading to an in-frame deletion in the first Pleckstrin homology subdomain of the GEF1 domain: p.(Thr1539_Lys1572del). To test for a distinct gestalt, facial characteristics of the family members and 41 previously published TRIO cases were systematically evaluated via GestaltMatcher. Computational analysis of the facial gestalt suggests a distinguishable facial TRIO-phenotype not outlined in the existing literature.

Spliceosomal helicases DDX41/SACY-1 and PRP22/MOG-5 both contribute to proofreading against proximal 3' splice site usage

RNA helicases drive necessary rearrangements and ensure fidelity during the pre-mRNA splicing cycle. DEAD-box helicase DDX41 has been linked to human disease and has recently been shown to interact with DEAH-box helicase PRP22 in the spliceosomal C* complex, yet its function in splicing remains unknown. Depletion of DDX41 homolog SACY-1 from somatic cells has been previously shown to lead to changes in alternative 3' splice site (3'ss) usage. Here, we show by transcriptomic analysis of published and novel data sets that SACY-1 perturbation causes a previously unreported pattern in alternative 3' splicing in introns with pairs of 3' splice sites ≤18 nt away from each other. We find that both SACY-1 depletion and the allele sacy-1(G533R) lead to a striking unidirectional increase in the usage of the proximal (upstream) 3'ss. We previously discovered a similar alternative splicing pattern between germline tissue and somatic tissue, in which there is a unidirectional increase in proximal 3'ss usage in the germline for ∼200 events; many of the somatic SACY-1 alternative 3' splicing events overlap with these developmentally regulated events. We generated targeted mutant alleles of the Caenorhabditis elegans homolog of PRP22, mog-5, in the region of MOG-5 that is predicted to interact with SACY-1 based on the human C* structure. These viable alleles, and a mimic of the myelodysplastic syndrome-associated allele DDX41(R525H), all promote the usage of proximal alternative adjacent 3' splice sites. We show that PRP22/MOG-5 and DDX41/SACY-1 have overlapping roles in proofreading the 3'ss.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing associated with therapeutic response to splicing inhibitor

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that jumonji domain containing 6, arginine demethylase, and lysine hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven human neuroblastoma. JMJD6 cooperates with MYC in cellular transformation of murine neural crest cells by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a 'molecular glue' that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

Altered CELF4 splicing factor enhances pancreatic neuroendocrine tumors aggressiveness influencing mTOR and everolimus response

Pancreatic neuroendocrine tumors (PanNETs) comprise a heterogeneous group of tumors with growing incidence. Recent molecular analyses provided a precise picture of their genomic and epigenomic landscape. Splicing dysregulation is increasingly regarded as a novel cancer hallmark influencing key tumor features. We have previously demonstrated that splicing machinery is markedly dysregulated in PanNETs. Here, we aimed to elucidate the molecular and functional implications of CUGBP ELAV-like family member 4 (CELF4), one of the most altered splicing factors in PanNETs. CELF4 expression was determined in 20 PanNETs, comparing tumor and non-tumoral adjacent tissue. An RNA sequencing (RNA-seq) dataset was analyzed to explore CELF4-linked interrelations among clinical features, gene expression, and splicing events. Two PanNET cell lines were employed to assess CELF4 function in vitro and in vivo. PanNETs display markedly upregulated CELF4 expression, which is closely associated with malignancy features, altered expression of key tumor players, and distinct splicing event profiles. Modulation of CELF4 influenced proliferation in vitro and reduced in vivo xenograft tumor growth. Interestingly, functional assays and RNA-seq analysis revealed that CELF4 silencing altered mTOR signaling pathway, enhancing the effect of everolimus. We demonstrate that CELF4 is dysregulated in PanNETs, where it influences tumor development and aggressiveness, likely by modulating the mTOR pathway, suggesting its potential as therapeutic target.

Splicing defects and CRISPR-Cas9 correction in isogenic homozygous photoreceptor precursors harboring clustered deep-intronic ABCA4 variants

Splicing defects from deep-intronic variants significantly contribute to the mutational spectrum in ABCA4-associated inherited retinal diseases, necessitating functional validation for their pathological classification. Typically, minigene assays in HEK293(T) can qualitatively assess splicing defects, yet they often fail to quantitatively reproduce the resulting mis-splicing patterns, leaving uncertainty on severity and pathogenicity. As a potential cellular model derived from patient cells, photoreceptor precursor cells (PPCs) play a pivotal role in assessing the severity of specific splicing mutations. Nevertheless, the accessibility of biosamples is commonly constrained, and their establishment is costly and laborious. In this study, we combined and investigated the use of a minigene assay and isogenic PPCs, as superior qualitative and more accessible cellular models for the assessment of splicing defects. Specifically, we focused on the clustered c.5196+1013A>G, c.5196+1056A>G, and c.5196+1216C>A deep-intronic variants in intron 36 of ABCA4, comparing their resulting (mis)splicing patterns in minigene-transfected cells and isogenic CRISPR-Cas9-knocked-in PPCs harboring these pathogenic variants in homozygous state. Moreover, we demonstrate the successful correction of these three splicing defects in homozygous mutant PPCs using a single pair of guide RNAs to target Cas9 cleavage, thereby identifying an efficient gene editing strategy for therapeutic applications.

Autophagy-dependent alternative splicing of ribosomal protein S24 produces a more stable isoform that aids in hypoxic cell survival

Cells remodel splicing and translation machineries to mount specialized gene expression responses to stress. Here, we show that hypoxic human cells in 2D and 3D culture models increase the relative abundance of a longer mRNA variant of ribosomal protein S24 (RPS24L) compared to a shorter mRNA variant (RPS24S) by favoring the inclusion of a 22 bp cassette exon. Mechanistically, RPS24L and RPS24S are induced and repressed, respectively, by distinct pathways in hypoxia: RPS24L is induced in an autophagy-dependent manner, while RPS24S is reduced by mTORC1 repression in a hypoxia-inducible factor-dependent manner. RPS24L produces a more stable protein isoform that aids in hypoxic cell survival and growth, which could be exploited by cancer cells in the tumor microenvironment.

RNA-binding proteins in pain

RNAs are meticulously controlled by proteins. Through direct and indirect associations, every facet in the brief life of an mRNA is subject to regulation. RNA-binding proteins (RBPs) permeate biology. Here, we focus on their roles in pain. Chronic pain is among the largest challenges facing medicine and requires new strategies. Mounting pharmacologic and genetic evidence obtained in pre-clinical models suggests fundamental roles for a broad array of RBPs. We describe their diverse roles that span RNA modification, splicing, stability, translation, and decay. Finally, we highlight opportunities to expand our understanding of regulatory interactions that contribute to pain signaling. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Regulation RNA in Disease and Development > RNA in Disease.

A novel KNL1 intronic splicing variant likely destabilizes the KMN complex, causing primary microcephaly

Primary microcephaly (MCPH) is an autosomal recessive disorder characterized by head circumference of at least two standard deviations below the mean. Biallelic variants in the kinetochore gene KNL1 is a known cause of MCPH4. KNL1 is the central component of the KNL1-MIS12-NSL1 (KMN) network, which acts as the signaling hub of the kinetochore and is required for correct chromosomal segregation during mitosis. We identified biallelic KNL1 variants in two siblings from a non-consanguineous family with microcephaly and intellectual disability. The two siblings carry a frameshift variant predicted to prematurely truncate the transcript and undergo nonsense mediated decay, and an intronic single nucleotide variant (SNV) predicted to disrupt splicing. An in vitro splicing assay and qPCR from blood-derived RNA confirmed that the intronic variant skips exon 23, significantly reducing levels of the canonical transcript. Protein modeling confirmed that absence of exon 23, an inframe exon, would disrupt a key interaction within the KMN network and likely destabilize the kinetochore signaling hub, disrupting mitosis. Therefore, this splicing variant is pathogenic and, in trans with a frameshift variant, causes the MCPH phenotype associated with KLN1. This finding furthers the association of splicing variants as a common pathogenic variant class for KNL1.

DDX41: exploring the roles of a versatile helicase

DDX41 is a DEAD-box helicase and is conserved across species. Mutations in DDX41 have been associated with myeloid neoplasms, including myelodysplastic syndrome and acute myeloid leukemia. Though its pathogenesis is not completely known, DDX41 has been shown to have many cellular roles, including in pre-mRNA splicing, innate immune sensing, ribosome biogenesis, translational regulation, and R-loop metabolism. In this review, we will summarize the latest understandings regarding the various roles of DDX41, as well as highlight challenges associated with drug development to target DDX41. Overall, understanding the molecular and cellular mechanisms of DDX41 could help develop novel therapeutic options for DDX41 mutation-related hematologic malignancies.

Identification of a novel splice-site WWOX variant with paternal uniparental isodisomy in a patient with infantile epileptic encephalopathy

WOREE syndrome is an early infantile epileptic encephalopathy characterized by drug-resistant seizures and severe psychomotor developmental delays. We report a case of a WWOX splice-site mutation with uniparental isodisomy. A 1-year and 7-month-old girl presented with nystagmus and epileptic seizures from early infancy, with no fixation or pursuit of vision. Physical examination revealed small deformities, such as swelling of both cheeks, folded fingers, rocking feet, and scoliosis. Brain imaging revealed slight hypoplasia of the cerebrum. Electroencephalogram showed focal paroxysmal discharges during the interictal phase of seizures. Vitamin B6 and zonisamide were administered for early infantile epileptic encephalopathy; however, the seizures were not relieved. Despite altering the type and dosage of antiepileptic drugs and ACTH therapy, the seizures were intractable. Whole-exome analysis revealed the homozygosity of WWOX(NM_016373.4):c.516+1G>A. The WWOX mRNA sequencing using peripheral blood RNA confirmed that exon 5 was homozygously deleted. Based on these results, the patient was diagnosed with WOREE syndrome at 5 months. The WWOX variant found in this study is novel and has never been reported before. WOREE syndrome being extremely rare, further case series and analyses of its pathophysiology are warranted.

Molecular Targets of the 5-Amido-Carboxamide Bumped Kinase Inhibitor BKI-1748 in Cryptosporidium parvum and HCT-8 Host Cells

Cryptosporidium parvum is an apicomplexan parasite causing persistent diarrhea in humans and animals. Issuing from target-based drug development, calcium-dependent protein kinase 1 inhibitors, collectively named bumped kinase inhibitors (BKIs), with excellent efficacies in vitro and in vivo have been generated. Some BKIs including BKI-1748 share a core structure with similarities to the first-generation antiprotozoal drug quinine, which is known to exert notorious side effects. Unlike quinine, BKI-1748 rapidly interfered with C. parvum proliferation in the human colon tumor (HCT) cell line HCT-8 cells and caused dramatic effects on the parasite ultrastructure. To identify putative BKI targets in C. parvum and in host cells, we performed differential affinity chromatography with cell-free extracts from non-infected and infected HCT-8 cells using BKI-1748 and quinine epoxy-activated sepharose columns followed by mass spectrometry. C. parvum proteins of interest were identified in eluates from columns coupled to BKI-1748, or in eluates from both BKI-1748 and quinine columns. However, no C. parvum proteins could be identified binding exclusively to BKI-1748. In contrast, 25 BKI-1748-specific binding proteins originating from HCT-8 cells were detected. Moreover, 29 C. parvum and 224 host cell proteins were identified in both BKI-1748 as well as in quinine eluates. In both C. parvum and host cells, the largest subset of binding proteins was involved in RNA binding and modification, with a focus on ribosomal proteins and proteins involved in RNA splicing. These findings extend previous results, showing that BKI-1748 interacts with putative targets involved in common, essential pathways such as translation and RNA processing.

Biallelic variants in SNUPN cause a limb girdle muscular dystrophy with myofibrillar-like features

Alterations in RNA-splicing are a molecular hallmark of several neurological diseases, including muscular dystrophies where mutations in genes involved in RNA metabolism or characterised by alterations in RNA splicing have been described. Here, we present five patients from two unrelated families with a limb-girdle muscular dystrophy (LGMD) phenotype carrying a biallelic variant in SNUPN gene. Snurportin-1, the protein encoded by SNUPN, plays an important role in the nuclear transport of small nuclear ribonucleoproteins (snRNPs), essential components of the spliceosome. We combine deep phenotyping, including clinical features, histopathology and muscle magnetic resonance image (MRI), with functional studies in patient-derived cells and muscle biopsies to demonstrate that variants in SNUPN are the cause of a new type of LGMD according to current definition. Moreover, an in vivo model in Drosophila melanogaster further supports the relevance of Snurportin-1 in muscle. SNUPN patients show a similar phenotype characterised by proximal weakness starting in childhood, restrictive respiratory dysfunction and prominent contractures, although interindividual variability in terms of severity even in individuals from the same family was found. Muscle biopsy showed myofibrillar-like features consisting of myotilin deposits and Z-disc disorganisation. MRI showed predominant impairment of paravertebral, vasti, sartorius, gracilis, peroneal and medial gastrocnemius muscles. Conservation and structural analyses of Snurportin-1 p.Ile309Ser variant suggest an effect in nuclear-cytosol snRNP trafficking. In patient-derived fibroblasts and muscle, cytoplasmic accumulation of snRNP components is observed, while total expression of Snurportin-1 and snRNPs remains unchanged, which demonstrates a functional impact of SNUPN variant in snRNP metabolism. Furthermore, RNA-splicing analysis in patients' muscle showed widespread splicing deregulation, in particular in genes relevant for muscle development and splicing factors that participate in the early steps of spliceosome assembly. In conclusion, we report that SNUPN variants are a new cause of limb girdle muscular dystrophy with specific clinical, histopathological and imaging features, supporting SNUPN as a new gene to be included in genetic testing of myopathies. These results further support the relevance of splicing-related proteins in muscle disorders.

Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

Sf3b4 mutation in Xenopus tropicalis causes RNA splicing defects followed by massive gene dysregulation that disrupt cranial neural crest development

Nager syndrome is a rare craniofacial and limb disorder characterized by midface retrusion, micrognathia, absent thumbs, and radial hypoplasia. This disorder results from haploinsufficiency of SF3B4 (splicing factor 3b, subunit 4) a component of the pre-mRNA spliceosomal machinery. The spliceosome is a complex of RNA and proteins that function together to remove introns and join exons from transcribed pre-mRNA. While the spliceosome is present and functions in all cells of the body, most spliceosomopathies - including Nager syndrome - are cell/tissue-specific in their pathology. In Nager syndrome patients, it is the neural crest (NC)-derived craniofacial skeletal structures that are primarily affected. To understand the pathomechanism underlying this condition, we generated a Xenopus tropicalis sf3b4 mutant line using the CRISPR/Cas9 gene editing technology. Here we describe the sf3b4 mutant phenotype at neurula, tail bud, and tadpole stages, and performed temporal RNA-sequencing analysis to characterize the splicing events and transcriptional changes underlying this phenotype. Our data show that while loss of one copy of sf3b4 is largely inconsequential in Xenopus tropicalis, homozygous deletion of sf3b4 causes major splicing defects and massive gene dysregulation, which disrupt cranial NC cell migration and survival, thereby pointing at an essential role of Sf3b4 in craniofacial development.

Recurrent "outsider" intronic variation in the SLC5A6 gene causes severe mixed axonal and demyelinating neuropathy, cyclic vomiting and optic atrophy in 3 families from Maghreb

Sodium dependent multivitamin transporter (SMVT) deficiency is a very rare autosomal recessive disorder characterized by multisystemic clinical manifestations due to combined biotin, panthotenic acid and lipoic acid deficiency. About 10 families have been described so far. Accurate diagnosis is crucial because of the possibility of a supplementation treatment with proven efficacy. Here we describe 4 new patients (3 additional families) originating from the same world region (Algeria, Maghreb). All patients, born form consanguineous parents, were homozygous carriers of the same intronic variation, outside of canonical sites, in the SLC5A6 gene encoding SMVT. RNA study in one family allowed confirming the pathogenic effect of the variation and re-classifying this variant of uncertain significance as pathogenic, opening the possibility of genetic counseling and treatment. The identification of the same variation in three distinct and apparently unrelated families is suggestive of a founder effect. The phenotype of all patients was very similar, with systematic optic atrophy (initially considered as a very rare sign), severe cyclic vomiting, and rapidly progressive mixed axonal and demyelinating sensory motor neuropathy.

Collisions of RNA polymerases behind the replication fork promote alternative RNA splicing in newly replicated chromatin

DNA replication produces a global disorganization of chromatin structure that takes hours to be restored. However, how these chromatin rearrangements affect the regulation of gene expression and the maintenance of cell identity is not clear. Here, we use ChOR-seq and ChrRNA-seq experiments to analyze RNA polymerase II (RNAPII) activity and nascent RNA synthesis during the first hours after chromatin replication in human cells. We observe that transcription elongation is rapidly reactivated in nascent chromatin but that RNAPII abundance and distribution are altered, producing heterogeneous changes in RNA synthesis. Moreover, this first wave of transcription results in RNAPII blockages behind the replication fork, leading to changes in alternative splicing. Altogether, our results deepen our understanding of how transcriptional programs are regulated during cell division and uncover molecular mechanisms that explain why chromatin replication is an important source of gene expression variability.

mRNA nuclear export: how mRNA identity features distinguish functional RNAs from junk transcripts

The division of the cellular space into nucleoplasm and cytoplasm promotes quality control mechanisms that prevent misprocessed mRNAs and junk RNAs from gaining access to the translational machinery. Here, we explore how properly processed mRNAs are distinguished from both misprocessed mRNAs and junk RNAs by the presence or absence of various 'identity features'.

Why U matters: detection and functions of pseudouridine modifications in mRNAs

The uridine modifications pseudouridine (Ψ), dihydrouridine, and 5-methyluridine are present in eukaryotic mRNAs. Many uridine-modifying enzymes are associated with human disease, underscoring the importance of uncovering the functions of uridine modifications in mRNAs. These modified uridines have chemical properties distinct from those of canonical uridines, which impact RNA structure and RNA-protein interactions. Ψ, the most abundant of these uridine modifications, is present across (pre-)mRNAs. Recent work has shown that many Ψs are present at intermediate to high stoichiometries that are likely conducive to function and at locations that are poised to influence pre-/mRNA processing. Technological innovations and mechanistic investigations are unveiling the functions of uridine modifications in pre-mRNA splicing, translation, and mRNA stability, which are discussed in this review.

RNA processing and degradation mechanisms shaping the mitochondrial transcriptome of budding yeasts

Yeast mitochondrial genes are expressed as polycistronic transcription units that contain RNAs from different classes and show great evolutionary variability. The promoters are simple, and transcriptional control is rudimentary. Posttranscriptional mechanisms involving RNA maturation, stability, and degradation are thus the main force shaping the transcriptome and determining the expression levels of individual genes. Primary transcripts are fragmented by tRNA excision by RNase P and tRNase Z, additional processing events occur at the dodecamer site at the 3' end of protein-coding sequences. groups I and II introns are excised in a self-splicing reaction that is supported by protein splicing factors encoded by the nuclear genes, or by the introns themselves. The 3'-to-5' exoribonucleolytic complex called mtEXO is the main RNA degradation activity involved in RNA turnover and processing, supported by an auxiliary 5'-to-3' exoribonuclease Pet127p. tRNAs and, to a lesser extent, rRNAs undergo several different base modifications. This complex gene expression system relies on the coordinated action of mitochondrial and nuclear genes and undergoes rapid evolution, contributing to speciation events. Moving beyond the classical model yeast Saccharomyces cerevisiae to other budding yeasts should provide important insights into the coevolution of both genomes that constitute the eukaryotic genetic system.

The eukaryotic translation initiation factor eIF4E unexpectedly acts in splicing thereby coupling mRNA processing with translation: eIF4E induces widescale splicing reprogramming providing system-wide connectivity between splicing, nuclear mRNA export and translation

Recent findings position the eukaryotic translation initiation factor eIF4E as a novel modulator of mRNA splicing, a process that impacts the form and function of resultant proteins. eIF4E physically interacts with the spliceosome and with some intron-containing transcripts implying a direct role in some splicing events. Moreover, eIF4E drives the production of key components of the splicing machinery underpinning larger scale impacts on splicing. These drive eIF4E-dependent reprogramming of the splicing signature. This work completes a series of studies demonstrating eIF4E acts in all the major mRNA maturation steps whereby eIF4E drives production of the RNA processing machinery and escorts some transcripts through various maturation steps. In this way, eIF4E couples the mRNA processing-export-translation axis linking nuclear mRNA processing to cytoplasmic translation. eIF4E elevation is linked to worse outcomes in acute myeloid leukemia patients where these activities are dysregulated. Understanding these effects provides new insight into post-transcriptional control and eIF4E-driven cancers.

Functional Characterization of Two Novel Intron 4 SERPING1 Gene Splice Site Pathogenic Variants in Families with Hereditary Angioedema

Variants that affect splice sites comprise 14.3% of all pathogenic variants in the SERPING1 gene; more than half of them are located outside the canonical sites. To make a clinical decision concerning patients with such variants, it is essential to know the exact way in which the effect of the variant would be realized. The optimal approach to determine the consequences is considered to be mRNA analysis. In the current study, we present the results of functional analysis of two previously non-described variants in the SERPING1 gene (NM_000062.3) affecting intron 4: c.686-1G>A and c.685+4dup, which were detected in members of two Russian families with autosomal dominant inheritance of angioedema type 1. Analysis of the patients' mRNA (extracted from whole blood) showed that the SERPING1(NM_000062.3):c.685+4dup variant leads to the loss of the donor splice site and the activation of the cryptic site in exon 4: r.710_745del (p.Gly217_Pro228del), while the SERPING1(NM_000062.3):c.686-1G>A variant leads to the skipping of exon 5: r.746_949del (p.Asp229_Ser296del).

Pathogenic Variants in USH1G/SANS Alter Protein Interaction with Pre-RNA Processing Factors PRPF6 and PRPF31 of the Spliceosome

Pre-mRNA splicing is an essential process orchestrated by the spliceosome, a dynamic complex assembled stepwise on pre-mRNA. We have previously identified that USH1G protein SANS regulates pre-mRNA splicing by mediating the intranuclear transfer of the spliceosomal U4/U6.U5 tri-snRNP complex. During this process, SANS interacts with the U4/U6 and U5 snRNP-specific proteins PRPF31 and PRPF6 and regulates splicing, which is disturbed by variants of USH1G/SANS causative for human Usher syndrome (USH), the most common form of hereditary deaf-blindness. Here, we aim to gain further insights into the molecular interaction of the splicing molecules PRPF31 and PRPF6 to the CENTn domain of SANS using fluorescence resonance energy transfer assays in cells and in silico deep learning-based protein structure predictions. This demonstrates that SANS directly binds via two distinct conserved regions of its CENTn to the two PRPFs. In addition, we provide evidence that these interactions occur sequentially and a conformational change of an intrinsically disordered region to a short α-helix of SANS CENTn2 is triggered by the binding of PRPF6. Furthermore, we find that pathogenic variants of USH1G/SANS perturb the binding of SANS to both PRPFs, implying a significance for the USH1G pathophysiology.

Trans-Regulation of Alternative PD-L1 mRNA Processing by CDK12 in Non-Small-Cell Lung Cancer Cells

Immunotherapy using checkpoint inhibitors targeting the interaction between PD-1 on T cells and PD-L1 on cancer cells has shown significant results in non-small-cell lung cancer (NSCLC). Not all patients respond to the therapy, and PD-L1 expression heterogeneity is proposed to be one determinant for this. The alternative processing of PD-L1 RNA, which depends on an alternative poly-A site in intron 4, generates a shorter mRNA variant (PD-L1v4) encoding soluble PD-L1 (sPD-L1), relative to the canonical PD-L1v1 mRNA encoding membrane-associated PD-L1 (mPD-L1). This study aimed to identify factors influencing the ratio between these two PD-L1 mRNAs in NSCLC cells. First, we verified the existence of the alternative PD-L1 RNA processing in NSCLC cells, and from in silico analyses, we identified a candidate list of regulatory factors. Examining selected candidates showed that CRISPR/Cas9-generated loss-of-function mutations in CDK12 increased the PD-L1v4/PD-L1v1 mRNA ratio and, accordingly, the sPD-L1/mPD-L1 balance. The CDK12/13 inhibitor THZ531 could also increase the PD-L1v4/PD-L1v1 mRNA ratio and impact the PD-L1 transcriptional response to IFN-γ stimulation. The fact that CDK12 regulates PD-L1 transcript variant formation in NSCLC cells is consistent with CDK12's role in promoting transcriptional elongation over intron-located poly-A sites. This study lays the groundwork for clinical investigations to delineate the implications of the CDK12-mediated balancing of sPD-L1 relative to mPD-L1 for immunotherapeutic responses in NSCLC.

Involvement of a Rarely Used Splicing SD2b Site in the Regulation of HIV-1 vif mRNA Production as Revealed by a Growth-Adaptive Mutation

We have previously reported an HIV-1 mutant designated NL-Y226tac that expresses Vif at an ultra-low level, being replication-defective in high-APOBEC3G cells, such as H9. It carries a synonymous mutation within the splicing SA1 site relative to its parental clone. In order to determine whether a certain mutant(s) emerges during multi-infection cycles, we maintained H9 cells infected with a relatively low or high input of NL-Y226tac for extended time periods. Unexpectedly, we reproducibly identified a g5061a mutation in the SD2b site in the two independent long-term culture experiments that partially increases Vif expression and replication ability. Importantly, the adaptive mutation g5061a was demonstrated to enhance vif mRNA production by activation of the SA1 site mediated through increasing usage of a rarely used SD2b site. In the long-term culture initiated by a high virus input, we additionally found a Y226Fttc mutation at the original Y226tac site in SA1 that fully restores Vif expression and replication ability. As expected, the adaptive mutation Y226Fttc enhances vif mRNA production through increasing the splicing site usage of SA1. Our results here revealed the importance of the SD2b nucleotide sequence in producing vif mRNA involved in the HIV-1 adaptation and of mutual antagonism between Vif and APOBEC3 proteins in HIV-1 adaptation/evolution and survival.

SPLASH: A statistical, reference-free genomic algorithm unifies biological discovery

Today's genomics workflows typically require alignment to a reference sequence, which limits discovery. We introduce a unifying paradigm, SPLASH (Statistically Primary aLignment Agnostic Sequence Homing), which directly analyzes raw sequencing data, using a statistical test to detect a signature of regulation: sample-specific sequence variation. SPLASH detects many types of variation and can be efficiently run at scale. We show that SPLASH identifies complex mutation patterns in SARS-CoV-2, discovers regulated RNA isoforms at the single-cell level, detects the vast sequence diversity of adaptive immune receptors, and uncovers biology in non-model organisms undocumented in their reference genomes: geographic and seasonal variation and diatom association in eelgrass, an oceanic plant impacted by climate change, and tissue-specific transcripts in octopus. SPLASH is a unifying approach to genomic analysis that enables expansive discovery without metadata or references.

Oligomeric, phosphorylated, and truncated tau and spliceosome pathology within the entorhinal-hippocampal connectome across stages of Alzheimer's disease

Neurofibrillary tangles (NFTs) contain abnormally phosphorylated tau proteins, which spread within components of the medial temporal lobe (MTL) memory circuit in Alzheimer's disease (AD). Here, we used quantitative immunohistochemistry to determine the density of posttranslational oligomeric (TOC1 and TNT1), phosphorylated (AT8), and late truncated (TauC3) tau epitopes within the MTL subfields including entorhinal cortex (EC) layer II, subiculum, Cornu Ammonis (CA) subfields, and dentate gyrus (DG) in subjects who died with a clinical diagnosis of no cognitive impairment (NCI), mild cognitive impairment (MCI), and AD. We also examined whether alterations of the nuclear alternative splicing protein, SRSF2, are associated with tau pathology. Although a significant increase in TOC1, TNT1, and AT8 neuron density occurred in the EC in MCI and AD, subicular, DG granule cell, and CA1 and CA3 densities were only significantly higher in AD. TauC3 counts were not different between connectome regions and clinical groups. SRSF2 intensity in AT8-positive cells decreased significantly in all regions independent of the clinical groups examined. CA1 and subicular AT8, TauC3, and oligomeric densities correlated across clinical groups. EC AT8 counts correlated with CA subfields and subicular and DG values across clinical groups. Oligomeric and AT8 CA1, EC, and subicular density correlated with Braak stage. Decreased nuclear SRSF2 in the presence of cytoplasmic phosphorylated tau suggests a dual-hit process in NFT formation within the entorhinal hippocampal connectome during the onset of AD. Although oligomeric and phosphorylated tau follow a stereotypical pattern, clinical disease stage determined density of tau deposition and not anatomic location within the entorhinal-hippocampal connectome.

Calorie restriction modulates the transcription of genes related to stress response and longevity in human muscle: The CALERIE study

The lifespan extension induced by 40% caloric restriction (CR) in rodents is accompanied by postponement of disease, preservation of function, and increased stress resistance. Whether CR elicits the same physiological and molecular responses in humans remains mostly unexplored. In the CALERIE study, 12% CR for 2 years in healthy humans induced minor losses of muscle mass (leg lean mass) without changes of muscle strength, but mechanisms for muscle quality preservation remained unclear. We performed high-depth RNA-Seq (387-618 million paired reads) on human vastus lateralis muscle biopsies collected from the CALERIE participants at baseline, 12- and 24-month follow-up from the 90 CALERIE participants randomized to CR and "ad libitum" control. Using linear mixed effect model, we identified protein-coding genes and splicing variants whose expression was significantly changed in the CR group compared to controls, including genes related to proteostasis, circadian rhythm regulation, DNA repair, mitochondrial biogenesis, mRNA processing/splicing, FOXO3 metabolism, apoptosis, and inflammation. Changes in some of these biological pathways mediated part of the positive effect of CR on muscle quality. Differentially expressed splicing variants were associated with change in pathways shown to be affected by CR in model organisms. Two years of sustained CR in humans positively affected skeletal muscle quality, and impacted gene expression and splicing profiles of biological pathways affected by CR in model organisms, suggesting that attainable levels of CR in a lifestyle intervention can benefit muscle health in humans.

Expression of circular RNAs in myelodysplastic neoplasms and their association with mutations in the splicing factor gene SF3B1

Mutations in the splicing factor 3b subunit 1 (SF3B1) gene are frequent in myelodysplastic neoplasms (MDS). Because the splicing process is involved in the production of circular RNAs (circRNAs), we investigated the impact of SF3B1 mutations on circRNA processing. Using RNA sequencing, we measured circRNA expression in CD34+ bone marrow MDS cells. We defined circRNAs deregulated in a heterogeneous group of MDS patients and described increased circRNA formation in higher-risk MDS. We showed that the presence of SF3B1 mutations did not affect the global production of circRNAs; however, deregulation of specific circRNAs was observed. Particularly, we demonstrated that strong upregulation of circRNAs processed from the zinc finger E-box binding homeobox 1 (ZEB1) transcription factor; this upregulation was exclusive to SF3B1-mutated patients and was not observed in those with mutations in other splicing factors or other recurrently mutated genes, or with other clinical variables. Furthermore, we focused on the most upregulated ZEB1-circRNA, hsa_circ_0000228, and, by its knockdown, we demonstrated that its expression is related to mitochondrial activity. Using microRNA analyses, we proposed miR-1248 as a direct target of hsa_circ_0000228. To conclude, we demonstrated that mutated SF3B1 leads to deregulation of ZEB1-circRNAs, potentially contributing to the defects in mitochondrial metabolism observed in SF3B1-mutated MDS.

Targeting PRMT1-mediated SRSF1 methylation to suppress oncogenic exon inclusion events and breast tumorigenesis

PRMT1 plays a vital role in breast tumorigenesis; however, the underlying molecular mechanisms remain incompletely understood. Herein, we show that PRMT1 plays a critical role in RNA alternative splicing, with a preference for exon inclusion. PRMT1 methylome profiling identifies that PRMT1 methylates the splicing factor SRSF1, which is critical for SRSF1 phosphorylation, SRSF1 binding with RNA, and exon inclusion. In breast tumors, PRMT1 overexpression is associated with increased SRSF1 arginine methylation and aberrant exon inclusion, which are critical for breast cancer cell growth. In addition, we identify a selective PRMT1 inhibitor, iPRMT1, which potently inhibits PRMT1-mediated SRSF1 methylation, exon inclusion, and breast cancer cell growth. Combination treatment with iPRMT1 and inhibitors targeting SRSF1 phosphorylation exhibits an additive effect of suppressing breast cancer cell growth. In conclusion, our study dissects a mechanism underlying PRMT1-mediated RNA alternative splicing. Thus, PRMT1 has great potential as a therapeutic target in breast cancer treatment.

Tridimensional Structural Analysis of Tau Isoforms Generated by Intronic Retention

Background

Tauopathies are a subset of neurodegenerative diseases characterized by abnormal tau inclusions. Recently, we have discovered a new, human specific, tau isoform termed W-tau that originates by intron 12 retention. Our preliminary data suggests this newly discovered W-tau isoform might prevent aberrant aggregation of other tau isoforms but is significantly downregulated in tauopathies such as Alzheimer's disease.

Objective

To accurately predict, examine, and understand tau protein structure and the conformational basis for the neuroprotective role of W-tau.

Methods

A tridimensional deep learning-based approach and in vitro polymerization assay was included to accurately predict, analyze, and understand tau protein structure and the conformational basis for the neuroprotective role of W-tau.

Results

Our findings demonstrate: a) the predicted protein tridimensionality structure of the tau isoforms raised by intron retention and their comparison with the other tau isoforms; b) the interaction of W-tau peptide (from W-tau isoform) with other tau isoforms; c) the effect of W-tau peptide in the polymerization of those tau isoforms.

Conclusions

This study supports the importance of the structure-function relationship on the neuroprotective behavior of W-tau inhibiting tau fibrillization in vitro.

Patterns of Gene Expression, Splicing, and Allele-Specific Expression Vary among Macular Tissues and Clinical Stages of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, and elucidating its underlying disease mechanisms is vital to the development of appropriate therapeutics. We identified differentially expressed genes (DEGs) and differentially spliced genes (DSGs) across the clinical stages of AMD in disease-affected tissue, the macular retina pigment epithelium (RPE)/choroid and the macular neural retina within the same eye. We utilized 27 deeply phenotyped donor eyes (recovered within a 6 h postmortem interval time) from Caucasian donors (60-94 years) using a standardized published protocol. Significant findings were then validated in an independent set of well-characterized donor eyes (n = 85). There was limited overlap between DEGs and DSGs, suggesting distinct mechanisms at play in AMD pathophysiology. A greater number of previously reported AMD loci overlapped with DSGs compared to DEGs between disease states, and no DEG overlap with previously reported loci was found in the macular retina between disease states. Additionally, we explored allele-specific expression (ASE) in coding regions of previously reported AMD risk loci, uncovering a significant imbalance in C3 rs2230199 and CFH rs1061170 in the macular RPE/choroid for normal eyes and intermediate AMD (iAMD), and for CFH rs1061147 in the macular RPE/choroid for normal eyes and iAMD, and separately neovascular AMD (NEO). Only significant DEGs/DSGs from the macular RPE/choroid were found to overlap between disease states. STAT1, validated between the iAMD vs. normal comparison, and AGTPBP1, BBS5, CERKL, FGFBP2, KIFC3, RORα, and ZNF292, validated between the NEO vs. normal comparison, revealed an intricate regulatory network with transcription factors and miRNAs identifying potential upstream and downstream regulators. Findings regarding the complement genes C3 and CFH suggest that coding variants at these loci may influence AMD development via an imbalance of gene expression in a tissue-specific manner. Our study provides crucial insights into the multifaceted genomic underpinnings of AMD (i.e., tissue-specific gene expression changes, potential splice variation, and allelic imbalance), which may open new avenues for AMD diagnostics and therapies specific to iAMD and NEO.

The Diverse Genomic Landscape of Diamond-Blackfan Anemia: Two Novel Variants and a Mini-Review

Diamond-Blackfan anemia (DBA) is a ribosomopathy characterized by bone marrow erythroid hypoplasia, which typically presents with severe anemia within the first months of life. DBA is typically attributed to a heterozygous mutation in a ribosomal protein (RP) gene along with a defect in the ribosomal RNA (rRNA) maturation or levels. Besides classic DBA, DBA-like disease has been described with variations in 16 genes (primarily in GATA1, followed by ADA2 alias CECR1, HEATR3, and TSR2). To date, more than a thousand variants have been reported in RP genes. Splice variants represent 6% of identifiable genetic defects in DBA, while their prevalence is 14.3% when focusing on pathogenic and likely pathogenic (P/LP) variants, thus highlighting the impact of such alterations in RP translation and, subsequently, in ribosome levels. We hereby present two cases with novel pathogenic splice variants in RPS17 and RPS26. Associations of DBA-related variants with specific phenotypic features and malignancies and the molecular consequences of pathogenic variations for each DBA-related gene are discussed. The determinants of the spontaneous remission, cancer development, variable expression of the same variants between families, and selectivity of RP defects towards the erythroid lineage remain to be elucidated.

Genome-wide detection of human intronic AG-gain variants located between splicing branchpoints and canonical splice acceptor sites

Human genetic variants that introduce an AG into the intronic region between the branchpoint (BP) and the canonical splice acceptor site (ACC) of protein-coding genes can disrupt pre-mRNA splicing. Using our genome-wide BP database, we delineated the BP-ACC segments of all human introns and found extreme depletion of AG/YAG in the [BP+8, ACC-4] high-risk region. We developed AGAIN as a genome-wide computational approach to systematically and precisely pinpoint intronic AG-gain variants within the BP-ACC regions. AGAIN identified 350 AG-gain variants from the Human Gene Mutation Database, all of which alter splicing and cause disease. Among them, 74% created new acceptor sites, whereas 31% resulted in complete exon skipping. AGAIN also predicts the protein-level products resulting from these two consequences. We performed AGAIN on our exome/genomes database of patients with severe infectious diseases but without known genetic etiology and identified a private homozygous intronic AG-gain variant in the antimycobacterial gene SPPL2A in a patient with mycobacterial disease. AGAIN also predicts a retention of six intronic nucleotides that encode an in-frame stop codon, turning AG-gain into stop-gain. This allele was then confirmed experimentally to lead to loss of function by disrupting splicing. We further showed that AG-gain variants inside the high-risk region led to misspliced products, while those outside the region did not, by two case studies in genes STAT1 and IRF7. We finally evaluated AGAIN on our 14 paired exome-RNAseq samples and found that 82% of AG-gain variants in high-risk regions showed evidence of missplicing. AGAIN is publicly available from https://hgidsoft.rockefeller.edu/AGAIN and https://github.com/casanova-lab/AGAIN.

Exploring the interplay between PARP1 and circRNA biogenesis and function

PARP1 (poly-ADP-ribose polymerase 1) is a multidomain protein with a flexible and self-folding structure that allows it to interact with a wide range of biomolecules, including nucleic acids and target proteins. PARP1 interacts with its target molecules either covalently via PARylation or non-covalently through its PAR moieties induced by auto-PARylation. These diverse interactions allow PARP1 to participate in complex regulatory circuits and cellular functions. Although the most studied PARP1-mediated functions are associated with DNA repair and cellular stress response, subsequent discoveries have revealed additional biological functions. Based on these findings, PARP1 is now recognized as a major modulator of gene expression. Several discoveries show that this multifunctional protein has been intimately connected to several steps of mRNA biogenesis, from transcription initiation to mRNA splicing, polyadenylation, export, and translation of mRNA to proteins. Nevertheless, our understanding of PARP1's involvement in the biogenesis of both coding and noncoding RNA, notably circular RNA (circRNA), remains restricted. In this review, we outline the possible roles of PARP1 in circRNA biogenesis. A full examination of the regulatory roles of PARP1 in nuclear processes with an emphasis on circRNA may reveal new avenues to control dysregulation implicated in the pathogenesis of several diseases such as neurodegenerative disorders and cancers. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Processing > Splicing Regulation/Alternative Splicing.

An overview of RB1 transcript alterations detected during retinoblastoma genetic screening

Identification of pathogenic RB1 variants aids in the clinical management of families with retinoblastoma. We routinely screen DNA for RB1 variants, but transcript analysis can also be used for variant screening, and to help decide variant pathogenicity. DNA was screened by conformation analysis followed by Sanger sequencing. Large deletion/insertions were detected by polymorphism analysis, MLPA and quantitative-PCR. Methylation-specific PCR was used to detect hypermethylation. RNA screening was performed when a DNA pathogenic variant was missing, or to determine effects on splicing.Two hundred and thirteen small coding variants were predicted to affect splicing in 207 patients. Splice donor (sd) variants were nearly twice as frequent as splice acceptor (sa) with the most affected positions being sd + 1 and sa-1. Some missense and nonsense codons altered splicing, while some splice consensus variants did not. Large deletion/insertions can disrupt splicing, but RNA analysis showed that some of these are more complex than indicated by DNA testing. RNA screening found pathogenic variants in 53.8% of samples where DNA analysis did not. RB1 splicing is altered by changes at consensus splice sites, some missense and nonsense codons, deep intronic changes and large deletion/insertions. Common alternatively spliced transcripts may complicate analysis. An effective molecular screening strategy would include RNA analysis to help determine pathogenicity.

RNA Editing in Cancer Progression

Coding and noncoding RNA molecules play their roles in ensuring cell function and tissue homeostasis in an ordered and systematic fashion. RNA chemical modifications can occur both at bases and ribose sugar, and, similarly to DNA and histone modifications, can be written, erased, and recognized by the corresponding enzymes, thus modulating RNA activities and fine-tuning gene expression programs. RNA editing is one of the most prevalent and abundant forms of post-transcriptional RNA modification in normal physiological processes. By altering the sequences of mRNAs, it makes them different from the corresponding genomic template. Hence, edited mRNAs can produce protein isoforms that are functionally different from the corresponding genome-encoded variants. Abnormalities in regulatory enzymes and changes in RNA-modification patterns are closely associated with the occurrence and development of various human diseases, including cancer. To date, the roles played by RNA modifications in cancer are gathering increasing interest. In this review, we focus on the role of RNA editing in cancer transformation and provide a new perspective on its impact on tumorigenesis, by regulating cell proliferation, differentiation, invasion, migration, stemness, metabolism, and drug resistance.

Contribution and therapeutic implications of retroelement insertions in ataxia telangiectasia

Certain classes of genetic variation still escape detection in clinical sequencing analysis. One such class is retroelement insertion, which has been reported as a cause of Mendelian diseases and may offer unique therapeutic implications. Here, we conducted retroelement profiling on whole-genome sequencing data from a cohort of 237 individuals with ataxia telangiectasia (A-T). We found 15 individuals carrying retroelement insertions in ATM, all but one of which integrated in noncoding regions. Systematic functional characterization via RNA sequencing, RT-PCR, and/or minigene splicing assays showed that 12 out of 14 intronic insertions led or contributed to ATM loss of function by exon skipping or activating cryptic splice sites. We also present proof-of-concept antisense oligonucleotides that suppress cryptic exonization caused by a deep intronic retroelement insertion. These results provide an initial systematic estimate of the contribution of retroelements to the genetic architecture of recessive Mendelian disorders as ∼2.1%-5.5%. Our study highlights the importance of retroelement insertions as causal variants and therapeutic targets in genetic diseases.

Spliceosomal introns in the diplomonad parasite Giardia duodenalis revisited

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain 'WB' genome and further elucidating Giardia's gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.