Splice site recognition - deciphering Exon-Intron transitions for genetic insights using Enhanced integrated Block-Level gated LSTM model

Bioinformatics is a contemporary interdisciplinary area focused on analyzing the growing number of genome sequences. Gene variants are differences in DNA sequences among individuals within a population. Splice site recognition is a crucial step in the process of gene expression, where the coding sequences of genes are joined together to form mature messenger RNA (mRNA). These genetic variants that disrupt genes are believed to be the primary reason for neuro-developmental disorders like ASD (Autism Spectrum Disorder) is a neuro-developmental disorder that is diagnosed in individuals, families, and society and occurs as the developmental delay in one among the hundred genes that are associated with these disorders. Missense variants, premature stop codons, or deletions alter both the quality and quantity of encoded proteins. Predicting genes within exons and introns presents main challenges, such as dealing with sequencing errors, short reads, incomplete genes, overlapping, and more. Although many traditional techniques have been utilized in creating an exon prediction system, the primary challenge lies in accurately identifying the length and spliced strand location classification of exons in conjunction with introns. From now on, the suggested approach utilizes a Deep Learning algorithm to analyze intricate and extensive genomic datasets. M-LSTM is utilized to categorize three binary combinations (EI as 1, IE as 2, and none as 3) using spliced DNA strands. The M-LSTM system is able to sequence extensive datasets, ensuring that long information can be stored without any impact on the current input or output. This enables it to recognize and address long-term connections and problems with rapidly increasing gradients. The proposed model is compared internally with Naïve Bayes and Random Forest to assess its efficacy. Additionally, the proposed model's performance is forecasted by utilizing probabilistic parameters like recall, F1-score, precision, and accuracy to assess the effectiveness of the proposed system.

Root Primordium Defective 1 Encodes an Essential PORR Protein Required for the Splicing of Mitochondria-Encoded Group II Introns and for Respiratory Complex I Biogenesis

Cellular respiration involves complex organellar metabolic activities that are pivotal for plant growth and development. Mitochondria contain their own genetic system (mitogenome, mtDNA), which encodes key elements of the respiratory machinery. Plant mtDNAs are notably larger than their counterparts in Animalia, with complex genome organization and gene expression characteristics. The maturation of the plant mitochondrial transcripts involves extensive RNA editing, trimming and splicing events. These essential processing steps rely on the activities of numerous nuclear-encoded cofactors, which may also play key regulatory roles in mitochondrial biogenesis and function and hence in plant physiology. Proteins that harbor the plant organelle RNA recognition (PORR) domain are represented in a small gene family in plants. Several PORR members, including WTF1, WTF9 and LEFKOTHEA, are known to act in the splicing of organellar group II introns in angiosperms. The AT4G33495 gene locus encodes an essential PORR protein in Arabidopsis, termed ROOT PRIMORDIUM DEFECTIVE 1 (RPD1). A null mutation of At.RPD1 causes arrest in early embryogenesis, while the missense mutant lines, rpd1.1 and rpd1.2, exhibit a strong impairment in root development and retarded growth phenotypes, especially under high-temperature conditions. Here, we further show that RPD1 functions in the splicing of introns that reside in the coding regions of various complex I (CI) subunits (i.e. nad2, nad4, nad5 and nad7), as well as in the maturation of the ribosomal rps3 pre-RNA in Arabidopsis mitochondria. The altered growth and developmental phenotypes and modified respiration activities are tightly correlated with respiratory chain CI defects in rpd1 mutants.

2-Desaza-annomontine (C81) impedes angiogenesis through reduced VEGFR2 expression derived from inhibition of CDC2-like kinases

Angiogenesis is a crucial process in the progression of various pathologies, like solid tumors, wet age-related macular degeneration, and chronic inflammation. Current anti-angiogenic treatments still have major drawbacks like limited efficacy in diseases that also rely on inflammation. Therefore, new anti-angiogenic approaches are sorely needed, and simultaneous inhibition of angiogenesis and inflammation is desirable. Here, we show that 2-desaza-annomontine (C81), a derivative of the plant alkaloid annomontine previously shown to inhibit endothelial inflammation, impedes angiogenesis by inhibiting CDC2-like kinases (CLKs) and WNT/β-catenin signaling. C81 reduced choroidal neovascularization in a laser-induced murine in vivo model, inhibited sprouting from vascular endothelial growth factor A (VEGF-A)-activated murine aortic rings ex vivo, and reduced angiogenesis-related activities of endothelial cells in multiple functional assays. This was largely phenocopied by CLK inhibitors and knockdowns, but not by inhibitors of the other known targets of C81. Mechanistically, CLK inhibition reduced VEGF receptor 2 (VEGFR2) mRNA and protein expression as well as downstream signaling. This was partly caused by a reduction of WNT/β-catenin pathway activity, as activating the pathway induced, while β-catenin knockdown impeded VEGFR2 expression. Surprisingly, alternative splicing of VEGFR2 was not detected. In summary, C81 and other CLK inhibitors could be promising compounds in the treatment of diseases that depend on angiogenesis and inflammation due to their impairment of both processes.

Assessing the reliability of point mutation as data augmentation for deep learning with genomic data

BACKGROUND

Deep neural networks (DNNs) have the potential to revolutionize our understanding and treatment of genetic diseases. An inherent limitation of deep neural networks, however, is their high demand for data during training. To overcome this challenge, other fields, such as computer vision, use various data augmentation techniques to artificially increase the available training data for DNNs. Unfortunately, most data augmentation techniques used in other domains do not transfer well to genomic data.

RESULTS

Most genomic data possesses peculiar properties and data augmentations may significantly alter the intrinsic properties of the data. In this work, we propose a novel data augmentation technique for genomic data inspired by biology: point mutations. By employing point mutations as substitutes for codons, we demonstrate that our newly proposed data augmentation technique enhances the performance of DNNs across various genomic tasks that involve coding regions, such as translation initiation and splice site detection.

CONCLUSION

Silent and missense mutations are found to positively influence effectiveness, while nonsense mutations and random mutations in non-coding regions generally lead to degradation. Overall, point mutation-based augmentations in genomic datasets present valuable opportunities for improving the accuracy and reliability of predictive models for DNA sequences.

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Gene knockout of RNA binding motif 5 in the brain alters RIMS2 protein homeostasis in the cerebellum and Hippocampus and exacerbates behavioral deficits after a TBI in mice

RNA binding motif 5 (RBM5) is a tumor suppressor in cancer but its role in the brain is unclear. We used conditional gene knockout (KO) mice to test if RBM5 inhibition in the brain affects chronic cortical brain tissue survival or function after a controlled cortical impact (CCI) traumatic brain injury (TBI). RBM5 KO decreased baseline contralateral hemispheric volume (p < 0.0001) and exacerbated ipsilateral tissue loss at 21 d after CCI in male mice vs. wild type (WT) (p = 0.0019). CCI injury, but not RBM5 KO, impaired beam balance performance (0-5d post-injury) and swim speed on the Morris Water Maze (MWM) (19-20d) (p < 0.0001). RBM5 KO was associated with mild learning impairment in female mice (p = 0.0426), reflected as a modest increase in escape latency early in training (14-18d post-injury). However, KO did not affect spatial memory at 19d post-injury in male or in female mice but it was impaired by CCI in females (p = 0.0061). RBM5 KO was associated with impaired visual function in male mice on the visible platform test at 20d post-injury (p = 0.0256). To explore signaling disturbances in KOs related to behavior, we first cross-referenced known brain-specific RBM5-regulated gene targets with genes in the curated RetNet database that impact vision. We then performed a secondary literature search on RBM5-regulated genes with a putative role in hippocampal function. Regulating synaptic membrane exocytosis 2 (RIMS) 2 was identified as a gene of interest because it regulates both vision and hippocampal function. Immunoprecipitation and western blot confirmed protein expression of a novel ~170 kDa RIMS2 variant in the cerebellum, and in the hippocampus, it was significantly increased in KO vs WT (p < 0.0001), and in a sex-dependent manner (p = 0.0390). Furthermore, male KOs had decreased total canonical RIMS2 levels in the cerebellum (p = 0.0027) and hippocampus (p < 0.0001), whereas female KOs had increased total RIMS1 levels in the cerebellum (p = 0.0389). In summary, RBM5 modulates brain function in mammals. Future work is needed to test if RBM5 dependent regulation of RIMS2 splicing effects vision and cognition, and to verify potential sex differences on behavior in a larger cohort of mice.

The molecular complexity of COL2A1 splicing variants and their significance in phenotype severity

Pathogenic single nucleotide variants (SNVs) found in the COL2A1 gene are associated with a broad range of skeletal dysplasias due to their impact on the structure and function of the Col2a1 protein. However, the molecular mechanisms of some nucleotide variants detected during diagnostic testing remain unclear. The interpretation of missense and splicing variants caused by SNVs poses a significant challenge for clinicians. In this work, we analyzed 22 splicing variants in the COL2A1 gene which have been found in patients with COL2A1-associated skeletal dysplasias. Using a minigene system, we investigated the impact of these SNVs on splicing and gained insights into their molecular mechanisms and genotype-phenotype correlations for each patient. The results of our study are very useful for improving the accuracy of diagnosis and the management of patients with skeletal dysplasias caused by SNVs in the COL2A1 gene.

Expression- and splicing-based multi-tissue transcriptome-wide association studies identified multiple genes for breast cancer by estrogen-receptor status

BACKGROUND

Although several transcriptome-wide association studies (TWASs) have been performed to identify genes associated with overall breast cancer (BC) risk, only a few TWAS have explored the differences in estrogen receptor-positive (ER+) and estrogen receptor-negative (ER-) breast cancer. Additionally, these studies were based on gene expression prediction models trained primarily in breast tissue, and they did not account for alternative splicing of genes.

METHODS

In this study, we utilized two approaches to perform multi-tissue TWASs of breast cancer by ER subtype: (1) an expression-based TWAS that combined TWAS signals for each gene across multiple tissues and (2) a splicing-based TWAS that combined TWAS signals of all excised introns for each gene across tissues. To perform this TWAS, we utilized summary statistics for ER + BC from the Breast Cancer Association Consortium (BCAC) and for ER- BC from a meta-analysis of BCAC and the Consortium of Investigators of Modifiers of BRCA1 and BRCA2 (CIMBA).

RESULTS

In total, we identified 230 genes in 86 loci that were associated with ER + BC and 66 genes in 29 loci that were associated with ER- BC at a Bonferroni threshold of significance. Of these genes, 2 genes associated with ER + BC at the 1q21.1 locus were located at least 1 Mb from published GWAS hits. For several well-studied tumor suppressor genes such as TP53 and CHEK2 which have historically been thought to impact BC risk through rare, penetrant mutations, we discovered that common variants, which modulate gene expression, may additionally contribute to ER + or ER- etiology.

CONCLUSIONS

Our study comprehensively examined how differences in common variation contribute to molecular differences between ER + and ER- BC and introduces a novel, splicing-based framework that can be used in future TWAS studies.

Tunnel ultrasound can guide the use of peritoneal dialysis catheter exit site relocation by external splicing and cuff removal in refractory tunnel infection

BACKGROUND

Peritoneal dialysis (PD) catheter related infections continue to be a major cause of morbidity and transfer to hemodialysis (HD) in PD patients. The treatment of tunnel infection (TI) could be challenging, especially when the infection involves the superficial cuff requiring the removal of the catheter. To spare the patient the loss of the catheter and the transfer to HD, several mini-invasive surgical techniques have been proposed as rescue therapy. Furthermore, nowadays, the rapid growth of digital technology has enormously increased the diagnostic sensibility of the echo signal allowing to accurately defines the extent of the infectious process along the PD catheter tunnel.

METHODS

Between 1st January 2020 and 31st December 2021 seven patients who underwent exit-site relocation by external splicing and cuff removal at our institution due to refractory TI were included in the study. All patients were followed until 12 months after the procedure. As soon as TI was defined refractory to the medical therapy, an ultrasonographic examination of the catheter tunnel was performed to define the extent of the infectious episode.

RESULTS

Among the 7 infectious episodes, 4 were caused by P. aeruginosa, and 3 by S. aureus. Around the superficial cuff the hypo/anechoic collections detected by ultrasounds showed a mean diameter of 3.05 ± 0.79 mm. The exit-site relocation by external splicing and cuff removal was successful in all cases (7/7, 100%).

CONCLUSIONS

In our experience the use of exit site relocation by external splicing and cuff removal as rescue therapy for TI with positive ultrasounds for TI limited to superficial cuff involvement and without secondary peritonitis, yielded to promising results with a success rate of 100%. This preliminary experience underlines the paramount usefulness of tunnel echography in accurately defining the extent of TI and, consequently, guiding the choice of the therapeutical approach in refractory TI.

Alternative splicing expands the clinical spectrum of NDUFS6-related mitochondrial disorders

PURPOSE

We describe 3 families with Charcot-Marie-Tooth neuropathy (CMT), harboring a homozygous NDUFS6 NM_004553.6:c.309+5G>A variant previously linked to fatal Leigh syndrome. We aimed to characterize clinically and molecularly the newly identified patients and understand the mechanism underlying their milder phenotype.

METHODS

The patients underwent extensive clinical examinations. Exome sequencing was done in 4 affected individuals. The functional effect of the c.309+5G>A variant was investigated in patient-derived EBV-transformed lymphoblasts at the complementary DNA, protein, and mitochondrial level. Alternative splicing was evaluated using complementary DNA long-read sequencing.

RESULTS

All patients presented with early-onset, slowly progressive axonal CMT, and nystagmus; some exhibited additional central nervous system symptoms. The c.309+5G>A substitution caused the expression of aberrantly spliced transcripts and negligible levels of the canonical transcript. Immunoblotting showed reduced levels of mutant isoforms. No detectable defects in mitochondrial complex stability or bioenergetics were found.

CONCLUSION

We expand the clinical spectrum of NDUFS6-related mitochondrial disorders to include axonal CMT, emphasizing the clinical and pathophysiologic overlap between these 2 clinical entities. This work demonstrates the critical role that alternative splicing may play in modulating the severity of a genetic disorder, emphasizing the need for careful consideration when interpreting splice variants and their implications on disease prognosis.

A transient in planta editing assay identifies specific binding of the splicing regulator PTB as a prerequisite for cassette exon inclusion

The dynamic interaction of RNA-binding proteins (RBPs) with their target RNAs contributes to the diversity of ribonucleoprotein (RNP) complexes that are involved in a myriad of biological processes. Identifying the RNP components at high resolution and defining their interactions are key to understanding their regulation and function. Expressing fusions between an RBP of interest and an RNA editing enzyme can result in nucleobase changes in target RNAs, representing a recent addition to experimental approaches for profiling RBP/RNA interactions. Here, we have used the MS2 protein/RNA interaction to test four RNA editing proteins for their suitability to detect target RNAs of RBPs in planta. We have established a transient test system for fast and simple quantification of editing events and identified the hyperactive version of the catalytic domain of an adenosine deaminase (hADARcd) as the most suitable editing enzyme. Examining fusions between homologs of polypyrimidine tract binding proteins (PTBs) from Arabidopsis thaliana and hADARcd allowed determining target RNAs with high sensitivity and specificity. Moreover, almost complete editing of a splicing intermediate provided insight into the order of splicing reactions and PTB dependency of this particular splicing event. Addition of sequences for nuclear localisation of the fusion protein increased the editing efficiency, highlighting this approach's potential to identify RBP targets in a compartment-specific manner. Our studies have established the editing-based analysis of interactions between RBPs and their RNA targets in a fast and straightforward assay, offering a new system to study the intricate composition and functions of plant RNPs in vivo.

Splicing variant of WDR37 in a case of Neurooculocardiogenitourinary syndrome

BACKGROUND

Neurooculocardiogenitourinary syndrome (NOCGUS), a multisystemic syndrome characterized by motor disorder, intellectual disability, seizures, abnormal brain structure, ocular diseases, and cardiac diseases, has been reported with missense variant of WD repeat-containing protein 37 (WDR37) in humans. This report aimed to identify the cause of NOCGUS in an affected patient.

CASE PRESENTATION

We identified a de novo intronic 4-bp deletion of WDR37, c.727-27_727-24del, which were predicted to cause abnormal splicing by SpliceAI, in the patient with NOCGUS. Reverse transcription polymerase chain reaction (RT-PCR) revealed intron retention of 63 base pairs before exon 10 in messenger RNA, which was predicted to insert 21 additional aberrant amino acids (p.S242_I243insLCQKKLKISRKCLFWPSLWQQ). The patient had novel phenotypes, anal atresia, and polycystic kidney, in addition to intellectual disability, seizures, cerebellar vermian anomaly, and coloboma, which are typical in NOCGUS. We did not observe motor impairments or cardiovascular anomalies.

CONCLUSION

This is the first reported case of NOCGUS with the splicing variant of WDR37, which manifests with distinctive but variable features. Our findings may expand a possible phenotypic expression of NOCGUS.

Therapeutic potential of alternative splicing in cardiovascular diseases

RNA splicing is an important RNA processing step required by multiexon protein-coding mRNAs and some noncoding RNAs. Precise RNA splicing is required for maintaining gene and cell function; however, mis-spliced RNA transcripts can lead to loss- or gain-of-function effects in human diseases. Mis-spliced RNAs induced by gene mutations or the dysregulation of splicing regulators may result in frameshifts, nonsense-mediated decay (NMD), or inclusion/exclusion of exons. Genetic animal models have characterised multiple splicing factors required for cardiac development or function. Moreover, sarcomeric and ion channel genes, which are closely associated with cardiovascular function and disease, are hotspots for AS. Here, we summarise splicing factors and their targets that are associated with cardiovascular diseases, introduce some therapies potentially related to pathological AS targets, and raise outstanding questions and future directions in this field.

cerevisiae Splicing Factor Luc7

Identification of splice sites is a critical step in pre-mRNA splicing since definition of the exon/intron boundaries controls what nucleotides are incorporated into mature mRNAs. The intron boundary with the upstream exon is initially identified through interactions with the U1 snRNP. This involves both base pairing between the U1 snRNA and the pre-mRNA as well as snRNP proteins interacting with the 5' splice site/snRNA duplex. In yeast, this duplex is buttressed by two conserved protein factors, Yhc1 and Luc7. Luc7 has three human paralogs (LUC7L, LUC7L2, and LUC7L3) which play roles in alternative splicing. What domains of these paralogs promote splicing at particular sites is not yet clear. Here, we humanized the zinc finger domains of the yeast Luc7 protein in order to understand their roles in splice site selection using reporter assays, transcriptome analysis, and genetic interactions. While we were unable to determine a function for the first zinc finger domain, humanization of the second zinc finger domain to mirror that found in LUC7L or LUC7L2 resulted in altered usage of nonconsensus 5' splice sites. In contrast, the corresponding zinc finger domain of LUC7L3 could not support yeast viability. Further, humanization of Luc7 can suppress mutation of the ATPase Prp28, which is involved in U1 release and exchange for U6 at the 5' splice site. Our work reveals a role for the second zinc finger of Luc7 in splice site selection and suggests that different zinc finger domains may have different ATPase requirements for release by Prp28.

A common cellular response to broad splicing perturbations is characterized by metabolic transcript downregulation driven by the Mdm2-p53 axis

Disruptions in core cellular processes elicit stress responses that drive cell-state changes leading to organismal phenotypes. Perturbations in the splicing machinery cause widespread mis-splicing, resulting in p53-dependent cell-state changes that give rise to cell-type-specific phenotypes and disease. However, a unified framework for how cells respond to splicing perturbations, and how this response manifests itself in nuanced disease phenotypes, has yet to be established. Here, we show that a p53-stabilizing Mdm2 alternative splicing event and the resulting widespread downregulation of metabolic transcripts are common events that arise in response to various splicing perturbations in both cellular and organismal models. Together, our results classify a common cellular response to splicing perturbations, put forth a new mechanism behind the cell-type-specific phenotypes that arise when splicing is broadly disrupted, and lend insight into the pleiotropic nature of the effects of p53 stabilization in disease.

The splicing factor SF3B1 is involved in brown adipocyte thermogenic activation

The ability of alternative splicing mechanisms to control gene expression is increasingly being recognized as relevant for adipose tissue function. The expression of SF3B1, a key component of the SF3B complex directly involved in spliceosome formation, was previously reported to be significantly induced in brown adipose tissue under cold-induced thermogenic activation. Here, we identify that noradrenergic cAMP-mediated thermogenic stimulation increases SF3B1 expression in brown and beige adipocytes. We further show that pladienolide-B, a drug that binds SF3B1 to inhibit pre-mRNA splicing by targeting the SF3B complex, down-regulates key components of the thermogenic machinery (e.g., UCP1 gene expression), differentially alters the expression of alternative splicing-regulated transcripts encoding molecular actors involved in the oxidative metabolism of brown adipocytes (e.g., peroxisome proliferator-activated receptor-gamma co-activator-alpha [PGC-1α] and cytochrome oxidase subunit 7a genes), and impairs the respiratory activity of brown adipocytes. Similar alterations were found in brown adipocytes with siRNA-mediated knockdown of SF3B1 protein levels. Our findings collectively indicate that SF3B1 is a key factor in the appropriate thermogenic activation of differentiated brown adipocytes. This work exemplifies the importance of splicing processes in adaptive thermogenesis and suggests that pharmacological tools, such as pladienolide-B, may be used to modulate brown adipocyte thermogenic activity.

Deep intronic variant causes aberrant splicing of ATP7A in a family with a variable occipital horn syndrome phenotype

Genetic variants in ATP7A are associated with a spectrum of X-linked disorders. In descending order of severity, these are Menkes disease, occipital horn syndrome, and X-linked distal spinal muscular atrophy. After 30 years of diagnostic investigation, we identified a deep intronic ATP7A variant in four males from a family affected to variable degrees by a predominantly skeletal phenotype, featuring bowing of long bones, elbow joints with restricted mobility which dislocate frequently, coarse curly hair, chronic diarrhoea, and motor coordination difficulties. Analysis of whole genome sequencing data from the Genomics England 100,000 Genomes Project following clinical re-evaluation identified a deep intronic ATP7A variant, which was predicted by SpliceAI to have a modest splicing effect. Using a mini-gene splicing assay, we determined that the intronic variant results in aberrant splicing. Sanger sequencing of patient cDNA revealed ATP7A transcripts with exon 5 skipping, or inclusion of a novel intron 4 pseudoexon. In both instances, frameshift leading to premature termination are predicted. Quantification of ATP7A mRNA transcripts using a qPCR assay indicated that the majority of transcripts (86.1 %) have non-canonical splicing, with 68.0 % featuring exon 5 skipping, and 18.1 % featuring the novel pseudoexon. We suggest that the variability of the phenotypes within the affected males results from the stochastic effects of splicing. This deep intronic variant, resulting in aberrant ATP7A splicing, expands the understanding of intronic variation on the ATP7A-related disease spectrum.

Erythromycin for myotonic dystrophy type 1: a multicentre, randomised, double-blind, placebo-controlled, phase 2 trial

Background

Myotonic dystrophy type 1 (DM1) is a devastating multisystemic disorder caused by a CTG repeat expansion in the DMPK gene, which subsequently triggers toxic RNA expression and dysregulated splicing. In a preclinical study, we demonstrated that erythromycin reduces the toxicity of abnormal RNA and ameliorates the aberrant splicing and motor phenotype in DM1 model mice.

Methods

This multicentre, randomised, double-blind, placebo-controlled, phase 2 trial was conducted at three centres in Japan to translate preclinical findings into practical applications in patients with DM1 by evaluating the safety and efficacy of erythromycin. Between Nov 29, 2019, and Jan 20, 2022, a total of 30 adult patients with DM1 were enrolled and randomly assigned in a 1:2:2 ratio to receive either placebo or erythromycin at two daily doses (500 mg or 800 mg) for 24 weeks. The primary outcome included the safety and tolerability of erythromycin. The secondary efficacy measures included splicing biomarkers, 6-min walk test results, muscle strength, and serum creatinine kinase (CK) values. This trial is registered with the Japan Registry of Clinical Trials, jRCT2051190069.

Findings

Treatment-related gastrointestinal symptoms occurred more frequently in the erythromycin group, but all adverse events were mild to moderate and resolved spontaneously. No serious safety concerns were identified. The CK levels from baseline to week 24 decreased in the overall erythromycin group compared with the placebo group (mean change of -6.4 U/L [SD 149] vs +182.8 [SD 228]), although this difference was not statistically significant (p = 0.070). Statistically significant improvements in the overall erythromycin treated groups compared to placebo were seen for two of the eleven splicing biomarkers that were each evaluated in half of the trial sample. These were MBNL1 (p = 0.048) and CACNA1S (p = 0.042).

Interpretation

Erythromycin demonstrated favourable safety and tolerability profiles in patients with DM1. A well-powered phase 3 trial is needed to evaluate efficacy, building on the preliminary findings from this study.

Funding

Japan Agency for Medical Research and Development.

Molecular diagnosis of cystic fibrosis by RNA obtained from nasal epithelial cells

BACKGROUND

The diagnosis of cystic fibrosis (CF) is established when characteristic clinical signs are coupled with biallelic CFTR pathogenic variants. No previously reported non-canonical splice site variants have to be considered as variants of uncertain significance unless their effect on splicing has been validated.

METHODS

Two variants identified by next-generation sequencing were evaluated. We assayed their effects on splicing employing RNA analysis and real-time expression quantification from RNA obtained from the nasal epithelial cells of a patient with clinically suspected CF and of two patients with milder phenotypes (CFTR-related disorders).

RESULTS

The variant c.164+2dup causes skipping of exon 2 (p.(Ser18_Glu54del)) and exon 2 plus 3 (p.(Ser18Argfs*16)) in CFTR mRNA. Exon 2 expression in the patient heterozygous for c.164+2dup was decreased to 7 % of the exon 2 expression in the controls. The synonymous variant c.1584G>A causes a partial skipping of exon 11. The exon 11 expression in the two patients heterozygous for this variant was 22 % and 42 % of that of the controls, respectively.

CONCLUSION

We conclude that variant c.164+2dup affects mRNA processing and can be considered a CF-causing variant. The results of the functional assay also showed that the p.(Glu528=) variant, usually categorized as a neutral variant based on epidemiological data, partially affects mRNA processing in our patients. This finding would allow us to reclassify the variant as a CFTR-related variant with incomplete penetrance. RNA obtained from nasal epithelial cells is an easy and accurate tool for CFTR functional studies in patients with unclassified splice variants.

Combined Bioinformatic and Splicing Analysis of Likely Benign Intronic and Synonymous Variants Reveals Evidence for Pathogenicity

Background

Current clinical variant analysis pipelines focus on coding variants and intronic variants within 10-20 bases of an exon-intron boundary that may affect splicing. The impact of newer splicing prediction algorithms combined with in vitro splicing assays on rare variants currently considered Benign/Likely Benign (B/LB) is unknown.

Methods

Exome sequencing data from 576 pediatric cancer patients enrolled in the Texas KidsCanSeq study were filtered for intronic or synonymous variants absent from population databases, predicted to alter splicing via SpliceAI (>0.20), and scored as potentially deleterious by CADD (>10.0). Total cellular RNA was extracted from monocytes and RT-PCR products analyzed. Subsequently, rare synonymous or intronic B/LB variants in a subset of genes submitted to ClinVar were similarly evaluated. Variants predicted to lead to a frameshifted splicing product were functionally assessed using an in vitro splicing reporter assay in HEK-293T cells.

Results

KidsCanSeq exome data analysis revealed a rare, heterozygous, intronic variant (NM_177438.3(DICER1):c.574-26A>G) predicted by SpliceAI to result in gain of a secondary splice acceptor site. The proband had a personal and family history of pleuropulmonary blastoma consistent with DICER1 syndrome but negative clinical sequencing reports. Proband RNA analysis revealed alternative DICER1 transcripts including the SpliceAI-predicted transcript.Similar bioinformatic analysis of synonymous or intronic B/LB variants (n=31,715) in ClinVar from 61 Mendelian disease genes yielded 18 variants, none of which could be scored by MaxEntScan. Eight of these variants were assessed (DICER1 n=4, CDH1 n=2, PALB2 n=2) using in vitro splice reporter assay and demonstrated abnormal splice products (mean 66%; range 6% to 100%). Available phenotypic information from submitting laboratories demonstrated DICER1 phenotypes in 2 families (1 variant) and breast cancer phenotypes for PALB2 in 3 families (2 variants).

Conclusions

Our results demonstrate the power of newer predictive splicing algorithms to highlight rare variants previously considered B/LB in patients with features of hereditary conditions. Incorporation of SpliceAI annotation of existing variant data combined with either direct RNA analysis or in vitro assays has the potential to identify disease-associated variants in patients without a molecular diagnosis.

Spliceosome mutations are associated with clinical response in a phase 1b/2 study of the PLK1 inhibitor onvansertib in combination with decitabine in relapsed or refractory acute myeloid leukemia

PLK1 is overexpressed in acute myeloid leukemia (AML). A phase 1b trial of the PLK1 inhibitor onvansertib (ONV) combined with decitabine (DAC) demonstrated initial safety and efficacy in patients with relapsed/refractory (R/R) AML. The current study aimed to identify molecular predictors of response to ONV + DAC in R/R AML patients. A total of 44 R/R AML patients were treated with ONV + DAC and considered evaluable for efficacy. Bone marrow (BM) samples were collected at baseline for genomic and transcriptomic analysis (n = 32). A 10-gene expression signature, predictive of response to ONV + DAC, was derived from the leading-edge genes of gene set enrichment analyses (GSEA). The gene signature was evaluated in independent datasets and used to identify associated mutated genes. Twenty percent of the patients achieved complete remission, with or without hematologic count recovery (CR/CRi), and 32% exhibited a ≥50% reduction in bone marrow blasts. Patients who responded to treatment had elevated mitochondrial function and OXPHOS. The gene signature was not associated with response to DAC alone in an independent dataset. By applying the signature to the BeatAML cohort (n = 399), we identified a positive association between predicted ONV + DAC response and mutations in splicing factors (SF). In the phase 1b/2 trial, patients with SF mutations (SRSF2, SF3B1) had a higher CR/CRi rate (50%) compared to those without SF mutations (9%). PLK1 inhibition with ONV in combination with DAC could be a potential therapy in R/R AML patients, particularly those with high OXPHOS gene expression and SF mutations.

[Pathophysiology of myelodysplastic syndromes]

During aging, the onset of mutations at low frequency in hematopoietic cells or clonal hematopoiesis of indeterminate significance favors the evolution towards hemopathies such as myelodysplastic syndromes or acute leukemias, but also cardiovascular diseases and other pathologies. Acute or chronic inflammation related to age influences the clonal evolution and the immune response. Conversely, mutated hematopoietic cells create an inflammatory bone marrow environment facilitating their expansion. Various pathophysiological mechanisms depending on the type of mutation produce the diversity of phenotypes. Identifying factors affecting clonal selection is mandatory to improve patient care.

Splicing alterations in pancreatic ductal adenocarcinoma: a new molecular landscape with translational potential

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most lethal cancers worldwide, mainly due to its late diagnosis and lack of effective therapies, translating into a low 5-year 12% survival rate, despite extensive clinical efforts to improve outcomes. International cooperative studies have provided informative multiomic landscapes of PDAC, but translation of these discoveries into clinical advances are lagging. Likewise, early diagnosis biomarkers and new therapeutic tools are sorely needed to tackle this cancer. The study of poorly explored molecular processes, such as splicing, can provide new tools in this regard. Alternative splicing of pre-RNA allows the generation of multiple RNA variants from a single gene and thereby contributes to fundamental biological processes by finely tuning gene expression. However, alterations in alternative splicing are linked to many diseases, and particularly to cancer, where it can contribute to tumor initiation, progression, metastasis and drug resistance. Splicing defects are increasingly being associated with PDAC, including both mutations or dysregulation of components of the splicing machinery and associated factors, and altered expression of specific relevant gene variants. Such disruptions can be a key element enhancing pancreatic tumor progression or metastasis, while they can also provide suitable tools to identify potential candidate biomarkers and discover new actionable targets. In this review, we aimed to summarize the current information about dysregulation of splicing-related elements and aberrant splicing isoforms in PDAC, and to describe their relationship with the development, progression and/or aggressiveness of this dismal cancer, as well as their potential as therapeutic tools and targets.

miR-6893-3p is a bonafide negative regulator of splicing activator, RNPS1

RNA-binding protein with serine-rich domain 1, RNPS1, is a global guardian of splicing fidelity and has implications in cervical cancer cell progression. We previously observed elevated RNPS1 expression in cervical cancer cells compared to normal cells. To understand the mechanisms that lead to the dysregulation of RNPS1 expression in cervical cancer cells, we focused on microRNAs. Using an in silico approach, we predicted potential miRNA candidates targeting RNPS1. Among the candidate miRNAs, we found miR-6893-3p as a potential regulator of RNPS1 expression. Interestingly, the expression of miR-6893-3p is downregulated in cervical cancer cells compared to normal cells and its level is negatively correlated with the expression of RNPS1. Further, qPCR, Western blot analysis, and luciferase reporter assay confirmed that miR-6893-3p negatively regulates RNPS1 in HeLa cells. In this regard, overexpression of miR-6893-3p suppresses the endogenous mRNA and protein levels of RNPS1 in HeLa cells. Further investigation revealed that miR-6893-3p mediated regulation of RNPS1 is dependent on the binding of miR-6893-3p to a microRNA response element in the 3'UTR of RNPS1 mRNA. Furthermore, mechanistic analysis showed that targeted negative regulation of RNPS1 by miR-6893-3p occurs via enhanced mRNA degradation. Collectively, our findings establish miR-6893-3p as an important player in the post-transcriptional regulation of RNPS1 in HeLa cells.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03761-2.

Prmt5 deficient mouse B cells display RNA processing complexity and slower colorectal tumor progression

Protein arginine methyltransferase 5 (Prmt5) is essential for normal B-cell development; however, the roles of Prmt5 in tumor-infiltrating B cells in tumor therapy have not been well elucidated. Here, we revealed that CD19-cre-Prmt5fl/fl (Prmt5cko) mice showed smaller tumor weights and volumes in the colorectal cancer mouse model; B cells expressed higher levels of Ccl22 and Il12a, which attracted T cells to the tumor site. Furthermore, we used direct RNA sequencing to comprehensively profile RNA processes in Prmt5 deletion B cells to explore underline mechanisms. We found significantly differentially expressed isoforms, mRNA splicing, poly(A) tail lengths, and m6A modification changes between the Prmt5cko and control groups. Cd74 isoform expressions might be regulated by mRNA splicing; the expression of two novel Cd74 isoforms was decreased, while one isoform was elevated in the Prmt5cko group, but the Cd74 gene expression showed no changes. We observed Ccl22, Ighg1, and Il12a expression was significantly increased in the Prmt5cko group, whereas Jak3 and Stat5b expression was decreased. Ccl22 and Ighg1 expression might be associated with poly(A) tail length, Jak3, Stat5b, and Il12a expression might be modulated by m6A modification. Our study demonstrated that Prmt5 regulates B-cell function through different mechanisms and supported the development of Prmt5-targeted antitumor treatments.

SGC-CLK-1: A chemical probe for the Cdc2-like kinases CLK1, CLK2, and CLK4

Small molecule modulators are important tools to study both basic biology and the complex signaling of protein kinases. The cdc2-like kinases (CLK) are a family of four kinases that have garnered recent interest for their involvement in a diverse set of diseases such as neurodegeneration, autoimmunity, and many cancers. Targeted medicinal chemistry around a CLK inhibitor hit identified through screening of a kinase inhibitor set against a large panel of kinases allowed us to identify a potent and selective inhibitor of CLK1, 2, and 4. Here, we present the synthesis, selectivity, and preliminary biological characterization of this compound - SGC-CLK-1 (CAF-170). We further show CLK2 has the highest binding affinity, and high CLK2 expression correlates with a lower IC50 in a screen of multiple cancer cell lines. Finally, we show that SGC-CLK-1 not only reduces serine arginine-rich (SR) protein phosphorylation but also alters SR protein and CLK2 subcellular localization in a reversible way. Therefore, we anticipate that this compound will be a valuable tool for increasing our understanding of CLKs and their targets, SR proteins, at the level of phosphorylation and subcellular localization.

ParSE-seq: A Calibrated Multiplexed Assay to Facilitate the Clinical Classification of Putative Splice-altering Variants

Background

Interpreting the clinical significance of putative splice-altering variants outside 2-base pair canonical splice sites remains difficult without functional studies.

Methods

We developed Parallel Splice Effect Sequencing (ParSE-seq), a multiplexed minigene-based assay, to test variant effects on RNA splicing quantified by high-throughput sequencing. We studied variants in SCN5A, an arrhythmia-associated gene which encodes the major cardiac voltage-gated sodium channel. We used the computational tool SpliceAI to prioritize exonic and intronic candidate splice variants, and ClinVar to select benign and pathogenic control variants. We generated a pool of 284 barcoded minigene plasmids, transfected them into Human Embryonic Kidney (HEK293) cells and induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), sequenced the resulting pools of splicing products, and calibrated the assay to the American College of Medical Genetics and Genomics scheme. Variants were interpreted using the calibrated functional data, and experimental data were compared to SpliceAI predictions. We further studied some splice-altering missense variants by cDNA-based automated patch clamping (APC) in HEK cells and assessed splicing and sodium channel function in CRISPR-edited iPSC-CMs.

Results

ParSE-seq revealed the splicing effect of 224 SCN5A variants in iPSC-CMs and 244 variants in HEK293 cells. The scores between the cell types were highly correlated (R2=0.84). In iPSCs, the assay had concordant scores for 21/22 benign/likely benign and 24/25 pathogenic/likely pathogenic control variants from ClinVar. 43/112 exonic variants and 35/70 intronic variants with determinate scores disrupted splicing. 11 of 42 variants of uncertain significance were reclassified, and 29 of 34 variants with conflicting interpretations were reclassified using the functional data. SpliceAI computational predictions correlated well with experimental data (AUC = 0.96). We identified 20 unique SCN5A missense variants that disrupted splicing, and 2 clinically observed splice-altering missense variants of uncertain significance had normal function when tested with the cDNA-based APC assay. A splice-altering intronic variant detected by ParSE-seq, c.1891-5C>G, also disrupted splicing and sodium current when introduced into iPSC-CMs at the endogenous locus by CRISPR editing.

Conclusions

ParSE-seq is a calibrated, multiplexed, high-throughput assay to facilitate the classification of candidate splice-altering variants.

The proper interplay between the expression of Spo11 splice isoforms and the structure of the pseudoautosomal region promotes XY chromosomes recombination

XY chromosome missegregation is relatively common in humans and can lead to sterility or the generation of aneuploid spermatozoa. A leading cause of XY missegregation in mammals is the lack of formation of double-strand breaks (DSBs) in the pseudoautosomal region (PAR), a defect that may occur in mice due to faulty expression of Spo11 splice isoforms. Using a knock-in (ki) mouse that expresses only the single Spo11β splice isoform, here we demonstrate that by varying the genetic background of mice, the length of chromatin loops extending from the PAR axis and the XY recombination proficiency varies. In spermatocytes of C57Spo11βki/- mice, in which loops are relatively short, recombination/synapsis between XY is fairly normal. In contrast, in cells of C57/129Spo11βki/- males where PAR loops are relatively long, formation of DSBs in the PAR (more frequently the Y-PAR) and XY synapsis fails at a high rate, and mice produce sperm with sex-chromosomal aneuploidy. However, if the entire set of Spo11 splicing isoforms is expressed by a wild type allele in the C57/129 background, XY recombination and synapsis is recovered. By generating a Spo11αki mouse model, we prove that concomitant expression of SPO11β and SPO11α isoforms, boosts DSB formation in the PAR. Based on these findings, we propose that SPO11 splice isoforms cooperate functionally in promoting recombination in the PAR, constraining XY asynapsis defects that may arise due to differences in the conformation of the PAR between mouse strains.

A panel of TDP-43-regulated splicing events verifies loss of TDP-43 function in amyotrophic lateral sclerosis brain tissue

TDP-43 dysfunction is a molecular hallmark of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). A major hypothesis of TDP-43 dysfunction in disease is the loss of normal nuclear function, resulting in impaired RNA regulation and the emergence of cryptic exons. Cryptic exons and differential exon usage are emerging as promising markers of lost TDP-43 function in addition to revealing biological pathways involved in neurodegeneration in ALS/FTD. In this brief report, we identified markers of TDP-43 loss of function by depleting TARDBP from post-mortem human brain pericytes, a manipulable in vitro primary human brain cell model, and identifying differential exon usage events with bulk RNA-sequencing analysis. We present these data in an interactive database (https://www.scotterlab.auckland.ac.nz/research-themes/tdp43-lof-db-v2/) together with seven other TDP-43-depletion datasets we meta-analysed previously, for user analysis of differential expression and splicing signatures. Differential exon usage events that were validated by qPCR were then compiled into a 'differential exon usage panel' with other well-established TDP-43 loss-of-function exon markers. This differential exon usage panel was investigated in ALS and control motor cortex tissue to verify whether, and to what extent, TDP-43 loss of function occurs in ALS. We find that profiles of TDP-43-regulated cryptic exons, changed exon usage and changed 3' UTR usage discriminate ALS brain tissue from controls, verifying that TDP-43 loss of function occurs in ALS. We propose that TDP-43-regulated splicing events that occur in brain tissue will have promise as predictors of disease.

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

Motivation

Although transcriptomics data is typically used to analyse mature spliced mRNA, recent attention has focused on jointly investigating spliced and unspliced (or precursor-) mRNA, which can be used to study gene regulation and changes in gene expression production. Nonetheless, most methods for spliced/unspliced inference (such as RNA velocity tools) focus on individual samples, and rarely allow comparisons between groups of samples (e.g., healthy vs. diseased). Furthermore, this kind of inference is challenging, because spliced and unspliced mRNA abundance is characterized by a high degree of quantification uncertainty, due to the prevalence of multi-mapping reads, i.e., reads compatible with multiple transcripts (or genes), and/or with both their spliced and unspliced versions.

Results

Here, we present DifferentialRegulation, a Bayesian hierarchical method to discover changes between experimental conditions with respect to the relative abundance of unspliced mRNA (over the total mRNA). We model the quantification uncertainty via a latent variable approach, where reads are allocated to their gene/transcript of origin, and to the respective splice version. We designed several benchmarks where our approach shows good performance, in terms of sensitivity and error control, versus state-of-the-art competitors. Importantly, our tool is flexible, and works with both bulk and single-cell RNA-sequencing data.

Availability and implementation

DifferentialRegulation is distributed as a Bioconductor R package.

Integration of photoperiodic and temperature cues by the circadian clock to regulate insect seasonal adaptations

Organisms adapt to unfavorable seasonal conditions to survive. These seasonal adaptations rely on the correct interpretation of environmental cues such as photoperiod, and temperature. Genetic studies in several organisms, including the genetic powerhouse Drosophila melanogaster, indicate that circadian clock components, such as period and timeless, are involved in photoperiodic-dependent seasonal adaptations, but our understanding of this process is far from complete. In particular, the role of temperature as a key factor to complement photoperiodic response is not well understood. The development of new sequencing technologies has proven extremely useful in understanding the plastic changes that the clock and other cellular components undergo in different environmental conditions, including changes in gene expression and alternative splicing. This article discusses the integration of photoperiod and temperature for seasonal biology as well as downstream molecular and cellular pathways involved in the regulation of physiological adaptations that occur with changing seasons. We focus our discussion on the current understanding of the involvement of the molecular clock and the circadian clock neuronal circuits in these adaptations in D. melanogaster.

Pathogenic Mis-splicing of CPEB4 in Schizophrenia

BACKGROUND

Schizophrenia (SCZ) is caused by an interplay of polygenic risk and environmental factors, which may alter regulators of gene expression leading to pathogenic misexpression of SCZ risk genes. The CPEB family of RNA-binding proteins (CPEB1-4) regulates translation of target RNAs (approximately 40% of overall genes). We previously identified CPEB4 as a key dysregulated translational regulator in autism spectrum disorder (ASD) because its neuronal-specific microexon (exon 4) is mis-spliced in ASD brains, causing underexpression of numerous ASD risk genes. The genetic factors and pathogenic mechanisms shared between SCZ and ASD led us to hypothesize CPEB4 mis-splicing in SCZ leading to underexpression of multiple SCZ-related genes.

METHODS

We performed MAGMA-enrichment analysis on Psychiatric Genomics Consortium genome-wide association study data and analyzed RNA sequencing data from the PsychENCODE Consortium. Reverse transcriptase polymerase chain reaction and Western blot were performed on postmortem brain tissue, and the presence/absence of antipsychotics was assessed through toxicological analysis. Finally, mice with mild overexpression of exon 4-lacking CPEB4 (CPEB4Δ4) were generated and analyzed biochemically and behaviorally.

RESULTS

First, we found enrichment of SCZ-associated genes for CPEB4-binder transcripts. We also found decreased usage of CPEB4 microexon in SCZ probands, which was correlated with decreased protein levels of CPEB4-target SCZ-associated genes only in antipsychotic-free individuals. Interestingly, differentially expressed genes fit those reported for SCZ, specifically in the SCZ probands with decreased CPEB4-microexon inclusion. Finally, we demonstrated that mice with mild overexpression of CPEB4Δ4 showed decreased protein levels of CPEB4-target SCZ genes and SCZ-linked behaviors.

CONCLUSIONS

We identified aberrant CPEB4 splicing and downstream misexpression of SCZ risk genes as a novel etiological mechanism in SCZ.

A heterozygous splicing variant IVS9-7A > T in intron 9 of the MAPT gene in a patient with right-temporal variant frontotemporal dementia with atypical 4 repeat tauopathy

Right temporal variant frontotemporal dementia, also called right-predominant semantic dementia, often has an unclear position within the framework of the updated diagnostic criteria for behavioral variant frontotemporal dementia or primary progressive aphasia. Recent studies have suggested that this population may be clinically, neuropathologically, and genetically distinct from those with behavioral variant frontotemporal dementia or left-predominant typical semantic variant primary progressive aphasia. Here we describe a Japanese case of right temporal variant frontotemporal dementia with novel heterozygous MAPT mutation Adenine to Thymidine in intervening sequence (IVS) 9 at position -7 from 3' splicing site of intron 9/exon 10 boundary (MAPT IVS9-7A > T). Postmortem neuropathological analysis revealed a predominant accumulation of 4 repeat tau, especially in the temporal lobe, amygdala, and substantia nigra, but lacked astrocytic plaques or tufted astrocytes. Immunoelectron microscopy of the tau filaments extracted from the brain revealed a ribbon-like structure. Moreover, a cellular MAPT splicing assay confirmed that this novel variant promoted the inclusion of exon 10, resulting in the predominant production of 4 repeat tau. These data strongly suggest that the MAPT IVS9-7 A > T variant found in our case is a novel mutation that stimulates the inclusion of exon 10 through alternative splicing of MAPT transcript and causes predominant 4 repeat tauopathy which clinically presents as right temporal variant frontotemporal dementia.

Multicenter clinical and functional evidence reclassifies a recurrent noncanonical filamin C splice-altering variant

BACKGROUND

Truncating variants in filamin C (FLNC) can cause arrhythmogenic cardiomyopathy (ACM) through haploinsufficiency. Noncanonical splice-altering variants may contribute to this phenotype.

OBJECTIVE

The purpose of this study was to investigate the clinical and functional consequences of a recurrent FLNC intronic variant of uncertain significance (VUS), c.970-4A>G.

METHODS

Clinical data in 9 variant heterozygotes from 4 kindreds were obtained from 5 tertiary health care centers. We used in silico predictors and functional studies with peripheral blood and patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Isolated RNA was studied by reverse transcription polymerase chain reaction. iPSC-CMs were further characterized at baseline and after nonsense-mediated decay (NMD) inhibition, using quantitative polymerase chain reaction (qPCR), RNA-sequencing, and cellular electrophysiology. American College of Medical Genetics and Genomics (ACMG) criteria were used to adjudicate variant pathogenicity.

RESULTS

Variant heterozygotes displayed a spectrum of disease phenotypes, spanning from mild ventricular dysfunction with palpitations to severe ventricular arrhythmias requiring device shocks or progressive cardiomyopathy requiring heart transplantation. Consistent with in silico predictors, the c.970-4A>G FLNC variant activated a cryptic splice acceptor site, introducing a 3-bp insertion containing a premature termination codon. NMD inhibition upregulated aberrantly spliced transcripts by qPCR and RNA-sequencing. Patch clamp studies revealed irregular spontaneous action potentials, increased action potential duration, and increased sodium late current in proband-derived iPSC-CMs. These findings fulfilled multiple ACMG criteria for pathogenicity.

CONCLUSION

Clinical, in silico, and functional evidence support the prediction that the intronic c.970-4A>G VUS disrupts splicing and drives ACM, enabling reclassification from VUS to pathogenic.

fruitless is sex-differentially spliced and is important for the courtship behavior and development of silkmoth Bombyx mori

Sexual dimorphisms of the brain play essential roles in successful reproduction. Silkmoth Bombyx mori exhibits extensive sexual differences in sexual behavior, as well as their morphology. Although the neural circuits that transmit information about sex pheromone in the male brain are extensively analyzed, the molecular mechanisms that regulate their development are still elusive. In the present study, we focused on the silkmoth ortholog of fruitless (fru) as a candidate gene that regulates sexual dimorphisms of the brain. fru transcripts were expressed from multiple promoters in various tissues, and brain-specific transcripts were sex-specifically spliced, in a manner similar to Drosophila. Interestingly, fru was highly expressed in the adult female brain and the male larval testis. Analysis of CRISPR/Cas9-mediated fru knockout strains revealed that fru plays important roles in survival during late larval and pupal stages, testis development, and adult sexual behavior. fru mutant males exhibited highly reduced levels of courtship and low copulation rate, indicating that fru plays significant roles in the sexual behavior of silkmoths, although it is not absolutely necessary for copulation. In the fru mutant males, sexually dimorphic pattern of the odorant receptor expression was impaired, possibly causing the defects in courtship behavior. These results provide important clues to elucidate the development of sexual dimorphisms of silkmoth brains, as well as the evolution of fruitless gene in insects.

The new world of RNA diagnostics and therapeutics

The 5th Workshop IRE on Translational Oncology was held in Rome (Italy) on 27-28 March at the IRCCS Regina Elena National Cancer Institute. This meeting entitled "The New World of RNA diagnostics and therapeutics" highlightes the significant progress in the RNA field made over the last years. Research moved from pure discovery towards the development of diagnostic biomarkers or RNA-base targeted therapies seeking validation in several clinical trials. Non-coding RNAs in particular have been the focus of this workshop due to their unique properties that make them attractive tools for the diagnosis and therapy of cancer.This report collected the presentations of many scientists from different institutions that discussed recent oncology research providing an excellent overview and representative examples for each possible application of RNA as biomarker, for therapy or to increase the number of patients that can benefit from precision oncology treatment.In particular, the meeting specifically emphasized two key features of RNA applications: RNA diagnostic (Blandino, Palcau, Sestito, Díaz Méndez, Cappelletto, Pulito, Monteonofrio, Calin, Sozzi, Cheong) and RNA therapeutics (Dinami, Marcia, Anastasiadou, Ryan, Fattore, Regazzo, Loria, Aharonov).

Transformer 2β regulates the alternative splicing of cell cycle regulatory genes to promote the malignant phenotype of ovarian cancer

Late-stage ovarian cancer (OC) has a poor prognosis and a high metastasis rate, but the underlying molecular mechanism is unclear. RNA binding proteins (RBPs) play important roles in posttranscriptional regulation in the contexts of neoplasia and tumor metastasis. In this study, we explored the molecular functions of a canonical RBP, Transformer 2β homolog (TRA2B), in cancer cells. TRA2B knockdown in HeLa cells and subsequent whole-transcriptome RNA sequencing (RNA-seq) analysis revealed the TRA2B-regulated alternative splicing (AS) profile. We disrupted TRA2B expression in epithelial OC cells and performed a series of experiments to confirm the resulting effects on OC cell proliferation, apoptosis and invasion. TRA2B-regulated AS was tightly associated with the mitotic cell cycle, apoptosis and several cancer pathways. Moreover, the expression of hundreds of genes was regulated by TRA2B, and these genes were enriched in the functions of cell proliferation, cell adhesion and angiogenesis, which are related to the malignant phenotype of OC. By integrating the alternatively spliced and differentially expressed genes, we found that AS events and gene expression were regulated independently. We then explored and validated the oncogenic functions of TRA2B by knocking down its expression in OC cells. The high TRA2B expression was associated with poor prognosis in patients with OC. In ovarian tissue, TRA2B expression showed a gradual increasing trend with increasing malignancy. We demonstrated the important roles of TRA2B in ovarian neoplasia and aggressive OC behaviors and identified the underlying molecular mechanisms, facilitating the targeted treatment of OC.

Post-transcriptional control of haemostatic genes: mechanisms and emerging therapeutic concepts in thrombo-inflammatory disorders

The haemostatic system is pivotal to maintaining vascular integrity. Multiple components involved in blood coagulation have central functions in inflammation and immunity. A derailed haemostasis is common in prevalent pathologies such as sepsis, cardiovascular disorders, and lately, COVID-19. Physiological mechanisms limit the deleterious consequences of a hyperactivated haemostatic system through adaptive changes in gene expression. While this is mainly regulated at the level of transcription, co- and posttranscriptional mechanisms are increasingly perceived as central hubs governing multiple facets of the haemostatic system. This layer of regulation modulates the biogenesis of haemostatic components, for example in situations of increased turnover and demand. However, they can also be 'hijacked' in disease processes, thereby perpetuating and even causally entertaining associated pathologies. This review summarizes examples and emerging concepts that illustrate the importance of posttranscriptional mechanisms in haemostatic control and crosstalk with the immune system. It also discusses how such regulatory principles can be used to usher in new therapeutic concepts to combat global medical threats such as sepsis or cardiovascular disorders.

A transposable element into the human long noncoding RNA CARMEN is a switch for cardiac precursor cell specification

AIMS

The major cardiac cell types composing the adult heart arise from common multipotent precursor cells. Cardiac lineage decisions are guided by extrinsic and cell-autonomous factors, including recently discovered long noncoding RNAs (lncRNAs). The human lncRNA CARMEN, which is known to dictate specification toward the cardiomyocyte (CM) and the smooth muscle cell (SMC) fates, generates a diversity of alternatively spliced isoforms.

METHODS AND RESULTS

The CARMEN locus can be manipulated to direct human primary cardiac precursor cells (CPCs) into specific cardiovascular fates. Investigating CARMEN isoform usage in differentiating CPCs represents therefore a unique opportunity to uncover isoform-specific functions in lncRNAs. Here, we identify one CARMEN isoform, CARMEN-201, to be crucial for SMC commitment. CARMEN-201 activity is encoded within an alternatively spliced exon containing a MIRc short interspersed nuclear element. This element binds the transcriptional repressor REST (RE1 Silencing Transcription Factor), targets it to cardiogenic loci, including ISL1, IRX1, IRX5, and SFRP1, and thereby blocks the CM gene program. In turn, genes regulating SMC differentiation are induced.

CONCLUSIONS

These data show how a critical physiological switch is wired by alternative splicing and functional transposable elements in a long noncoding RNA. They further demonstrated the crucial importance of the lncRNA isoform CARMEN-201 in SMC specification during heart development.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

Development of MAPT S305 mutation models exhibiting elevated 4R tau expression, resulting in altered neuronal and astrocytic function

Due to the importance of 4R tau in the pathogenicity of primary tauopathies, it has been challenging to model these diseases in iPSC-derived neurons, which express very low levels of 4R tau. To address this problem we have developed a panel of isogenic iPSC lines carrying the MAPT splice-site mutations S305S, S305I or S305N, derived from four different donors. All three mutations significantly increased the proportion of 4R tau expression in iPSC-neurons and astrocytes, with up to 80% 4R transcripts in S305N neurons from as early as 4 weeks of differentiation. Transcriptomic and functional analyses of S305 mutant neurons revealed shared disruption in glutamate signaling and synaptic maturity, but divergent effects on mitochondrial bioenergetics. In iPSC-astrocytes, S305 mutations induced lysosomal disruption and inflammation and exacerbated internalization of exogenous tau that may be a precursor to the glial pathologies observed in many tauopathies. In conclusion, we present a novel panel of human iPSC lines that express unprecedented levels of 4R tau in neurons and astrocytes. These lines recapitulate previously characterized tauopathy-relevant phenotypes, but also highlight functional differences between the wild type 4R and mutant 4R proteins. We also highlight the functional importance of MAPT expression in astrocytes. These lines will be highly beneficial to tauopathy researchers enabling a more complete understanding of the pathogenic mechanisms underlying 4R tauopathies across different cell types.

Circular RNA regulation and function in drug seeking phenotypes

Drug overdoses have increased dramatically in the United States over the last decade where they are now the leading cause of accidental death. To develop efficient therapeutic options for decreasing drug consumption and overdose risk, it is critical to understand the neurobiological changes induced by drug exposure. Chronic systemic exposure to all drug classes, including opioids, psychostimulants, nicotine, cannabis, and alcohol, induces profound molecular neuroadaptations within the central nervous system that may reveal crucial information about the lasting effects that these substances impart on brain cells. Transcriptome analyses of messenger RNAs (mRNAs) have identified gene patterns in the brain that result from exposure to various classes of drugs. However, mRNAs represent only a small fraction of the RNA within the cell, and drug exposure also impacts other classes of RNA that are largely understudied, especially circular RNAs. Circular RNAs (circRNAs) are a naturally occurring RNA species formed from back-splicing events during mRNA processing and are enriched in the nervous system. circRNAs are a pleiotropic class of RNAs and have a diverse impact on cellular function, with putative functions including regulation of mRNA transcription, protein translation, microRNA sponging, and sequestration of RNA-binding proteins. Recent studies have demonstrated that circRNAs can modulate cognition and are regulated in the brain in response to drug exposure, yet very few studies have explored the contribution of circRNAs to drug seeking phenotypes. In this review, we will provide an overview of the mechanisms of circRNA function in the cell to highlight how drug-induced circRNA dysregulation may impact the molecular substrates that mediate drug seeking behavior and the current studies that have reported drug-induced dysregulation of circRNAs in the brain. Furthermore, we will discuss how principles of circRNA biology can be adapted to study circRNAs in models of drug exposure and seek to provide further insight into the neurobiology of addiction.

Introme accurately predicts the impact of coding and noncoding variants on gene splicing, with clinical applications

Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .

Autosomal recessive congenital cataract is associated with a novel 4-bp splicing deletion mutation in a novel C10orf71 human gene

Congenital cataract is one of the most genetically heterogeneous ocular conditions with different genes involved in its etiology. Here, we describe the analysis of a new candidate gene of a congenital bilateral cataract associated with polymalformative syndrome, moderate global developmental delay, microcephaly, axial hypotonia, intrauterine growth restriction and facial dysmorphism for two affected siblings. Molecular analysis included exome sequencing and genome wide homozygosity mapping revealed a region of homozygosity shared by the two affected siblings at 10q11.23. The new C10orf71 gene was included in this interval and direct sequencing of this gene revealed an already described homozygous c. 2123T > G mutation (p. L708R) for the two affected subjects. Interestingly, we revealed in contrast a 4-bp deletion on the 3'-splicing acceptor site of intron 3-exon 4, namely defined as IVS3-5delGCAA. The C10Orf71 gene expression analysis using RT-PCR showed an expression pattern in different fetal organs and tissues as well as in leukocytes and confirmed that the IVS3-5delGCAA deletion of the C10orf71 gene is a splicing mutation responsible for the shortening of the C10orf71 protein in the two related patients. The C10orf71 gene has not been described to date as associated to the autosomal recessive phenotype.

Sf3b4 regulates chromatin remodeler splicing and Hox expression

SF3B proteins form a heptameric complex in the U2 small nuclear ribonucleoprotein, essential for pre-mRNA splicing. Heterozygous pathogenic variants in human SF3B4 are associated with head, face, limb, and vertebrae defects. Using the CRISPR/Cas9 system, we generated mice with constitutive heterozygous deletion of Sf3b4 and showed that mutant embryos have abnormal vertebral development. Vertebrae abnormalities were accompanied by changes in levels and expression pattern of Hox genes in the somites. RNA sequencing analysis of whole embryos and somites of Sf3b4 mutant and control litter mates revealed increased expression of other Sf3b4 genes. However, the mutants exhibited few differentially expressed genes and a large number of transcripts with differential splicing events (DSE), predominantly increased exon skipping and intron retention. Transcripts with increased DSE included several genes involved in chromatin remodeling that are known to regulate Hox expression. Our study confirms that Sf3b4 is required for normal vertebrae development and shows, for the first time, that like Sf3b1, Sf3b4 also regulates Hox expression. We propose that abnormal splicing of chromatin remodelers is primarily responsible for vertebral defects found in Sf3b4 heterozygous mutant embryos.

Differential splicing of neuronal genes in a Trem2*R47H mouse model mimics alterations associated with Alzheimer's disease

BACKGROUND

Molecular characterization of late-onset Alzheimer's disease (LOAD), the leading cause of age-related dementia, has revealed transcripts, proteins, and pathway alterations associated with disease. Assessing these postmortem signatures of LOAD in experimental model systems can further elucidate their relevance to disease origins and progression. Model organisms engineered with human genetic factors further link these signatures to disease-associated variants, especially when studies are designed to leverage homology across species. Here we assess differential gene splicing patterns in aging mouse models carrying humanized APOE4 and/or the Trem2*R47H variant on a C57BL/6J background. We performed a differential expression of gene (DEG) and differential splicing analyses on whole brain transcriptomes at multiple ages. To better understand the difference between differentially expressed and differentially spliced genes, we evaluated enrichment of KEGG pathways and cell-type specific gene signatures of the adult brain from each alteration type. To determine LOAD relevance, we compared differential splicing results from mouse models with multiple human AD splicing studies.

RESULTS

We found that differentially expressed genes in Trem2*R47H mice were significantly enriched in multiple AD-related pathways, including immune response, osteoclast differentiation, and metabolism, whereas differentially spliced genes were enriched for neuronal related functions, including GABAergic synapse and glutamatergic synapse. These results were reinforced by the enrichment of microglial genes in DEGs and neuronal genes in differentially spliced genes in Trem2*R47H mice. We observed significant overlap between differentially spliced genes in Trem2*R47H mice and brains from human AD subjects. These effects were absent in APOE4 mice and suppressed in APOE4.Trem2*R47H double mutant mice relative to Trem2*R47H mice.

CONCLUSIONS

The cross-species observation that alternative splicing observed in LOAD are present in Trem2*R47H mouse models suggests a novel link between this candidate risk gene and molecular signatures of LOAD in neurons and demonstrates how deep molecular analysis of new genetic models links molecular disease outcomes to a human candidate gene.

The uprise of RNA biology in neuroendocrine neoplasms: altered splicing and RNA species unveil translational opportunities

Neuroendocrine neoplasms (NENs) comprise a highly heterogeneous group of tumors arising from the diffuse neuroendocrine system. NENs mainly originate in gastrointestinal, pancreatic, and pulmonary tissues, and despite being rare, show rising incidence. The molecular mechanisms underlying NEN development are still poorly understood, although recent studies are unveiling their genomic, epigenomic and transcriptomic landscapes. RNA was originally considered as an intermediary between DNA and protein. Today, compelling evidence underscores the regulatory relevance of RNA processing, while new RNA molecules emerge with key functional roles in core cell processes. Indeed, correct functioning of the interrelated complementary processes comprising RNA biology, its processing, transport, and surveillance, is essential to ensure adequate cell homeostasis, and its misfunction is related to cancer at multiple levels. This review is focused on the dysregulation of RNA biology in NENs. In particular, we survey alterations in the splicing process and available information implicating the main RNA species and processes in NENs pathology, including their role as biomarkers, and their functionality and targetability. Understanding how NENs precisely (mis)behave requires a profound knowledge at every layer of their heterogeneity, to help improve NEN management. RNA biology provides a wide spectrum of previously unexplored processes and molecules that open new avenues for NEN detection, classification and treatment. The current molecular biology era is rapidly evolving to facilitate a detailed comprehension of cancer biology and is enabling the arrival of personalized, predictive and precision medicine to rare tumors like NENs.

Clinical and functional characterization of a long survivor congenital titinopathy patient with a novel metatranscript-only titin variant

Congenital titinopathies are an emerging group of a potentially severe form of congenital myopathies caused by biallelic mutations in titin, encoding the largest existing human protein involved in the formation and stability of sarcomeres. In this study we describe a patient with a congenital myopathy characterized by multiple contractures, a rigid spine, non progressive muscular weakness, and a novel homozygous TTN pathogenic variant in a metatranscript-only exon: the c.36400A > T, p.Lys12134*. Muscle biopsies showed increased internalized nuclei, variability in fiber size, mild fibrosis, type 1 fiber predominance, and a slight increase in the number of satellite cells. RNA studies revealed the retention of intron 170 and 171 in the open reading frame, and immunoflourescence and western blot studies, a normal titin content. Single fiber functional studies showed a slight decrease in absolute maximal force and a cross-sectional area with no decreases in tension, suggesting that weakness is not sarcomere-based but due to hypotrophy. Passive properties of single fibers were not affected, but the observed increased calcium sensitivity of force generation might contribute to the contractural phenotype and rigid spine of the patient. Our findings provide evidence for a pathogenic, causative role of a metatranscript-only titin variant in a long survivor congenital titinopathy patient with distal arthrogryposis and rigid spine.

The era of cryptic exons: implications for ALS-FTD

TDP-43 is an RNA-binding protein with a crucial nuclear role in splicing, and mislocalises from the nucleus to the cytoplasm in a range of neurodegenerative disorders. TDP-43 proteinopathy spans a spectrum of incurable, heterogeneous, and increasingly prevalent neurodegenerative diseases, including the amyotrophic lateral sclerosis and frontotemporal dementia disease spectrum and a significant fraction of Alzheimer's disease. There are currently no directed disease-modifying therapies for TDP-43 proteinopathies, and no way to distinguish who is affected before death. It is now clear that TDP-43 proteinopathy leads to a number of molecular changes, including the de-repression and inclusion of cryptic exons. Importantly, some of these cryptic exons lead to the loss of crucial neuronal proteins and have been shown to be key pathogenic players in disease pathogenesis (e.g., STMN2), as well as being able to modify disease progression (e.g., UNC13A). Thus, these aberrant splicing events make promising novel therapeutic targets to restore functional gene expression. Moreover, presence of these cryptic exons is highly specific to patients and areas of the brain affected by TDP-43 proteinopathy, offering the potential to develop biomarkers for early detection and stratification of patients. In summary, the discovery of cryptic exons gives hope for novel diagnostics and therapeutics on the horizon for TDP-43 proteinopathies.

Metformin and simvastatin exert additive antitumour effects in glioblastoma via senescence-state: clinical and translational evidence

BACKGROUND

Glioblastoma is one of the most devastating and incurable cancers due to its aggressive behaviour and lack of available therapies, being its overall-survival from diagnosis ∼14-months. Thus, identification of new therapeutic tools is urgently needed. Interestingly, metabolism-related drugs (e.g., metformin/statins) are emerging as efficient antitumour agents for several cancers. Herein, we evaluated the in vitro/in vivo effects of metformin and/or statins on key clinical/functional/molecular/signalling parameters in glioblastoma patients/cells.

METHODS

An exploratory-observational-randomized retrospective glioblastoma patient cohort (n = 85), human glioblastoma/non-tumour brain human cells (cell lines/patient-derived cell cultures), mouse astrocytes progenitor cell cultures, and a preclinical xenograft glioblastoma mouse model were used to measure key functional parameters, signalling-pathways and/or antitumour progression in response to metformin and/or simvastatin.

FINDINGS

Metformin and simvastatin exerted strong antitumour actions in glioblastoma cell cultures (i.e., proliferation/migration/tumoursphere/colony-formation/VEGF-secretion inhibition and apoptosis/senescence induction). Notably, their combination additively altered these functional parameters vs. individual treatments. These actions were mediated by the modulation of key oncogenic signalling-pathways (i.e., AKT/JAK-STAT/NF-κB/TGFβ-pathways). Interestingly, an enrichment analysis uncovered a TGFβ-pathway activation, together with AKT inactivation, in response to metformin + simvastatin combination, which might be linked to an induction of the senescence-state, the associated secretory-phenotype, and to the dysregulation of spliceosome components. Remarkably, the antitumour actions of metformin + simvastatin combination were also observed in vivo [i.e., association with longer overall-survival in human, and reduction in tumour-progression in a mouse model (reduced tumour-size/weight/mitosis-number, and increased apoptosis)].

INTERPRETATION

Altogether, metformin and simvastatin reduce aggressiveness features in glioblastomas, being this effect significantly more effective (in vitro/in vivo) when both drugs are combined, offering a clinically relevant opportunity that should be tested for their use in humans.

FUNDING

Spanish Ministry of Science, Innovation and Universities; Junta de Andalucía; CIBERobn (CIBER is an initiative of Instituto de Salud Carlos III, Spanish Ministry of Health, Social Services and Equality).