Age-dependent regulation of ELP1 exon 20 splicing in Familial Dysautonomia by RNA Polymerase II kinetics and chromatin structure.

Familial Dysautonomia (FD) is a rare disease caused by ELP1 exon 20 skipping. Here we clarify the role of RNA Polymerase II (RNAPII) and chromatin on this splicing event. A slow RNAPII mutant and chromatin-modifying chemicals that reduce the rate of RNAPII elongation induce exon skipping whereas chemicals that create a more relaxed chromatin exon inclusion. In the brain of a mouse transgenic for the human FD-ELP1 we observed on this gene an age-dependent decrease in the RNAPII density profile that was most pronounced on the alternative exon, a robust increase in the repressive marks H3K27me3 and H3K9me3 and a decrease of H3K27Ac, together with a progressive reduction in ELP1 exon 20 inclusion level. In HEK 293T cells, selective drug-induced demethylation of H3K27 increased RNAPII elongation on ELP1 and SMN2, promoted the inclusion of the corresponding alternative exons, and, by RNA-sequencing analysis, induced changes in several alternative splicing events. These data suggest a co-transcriptional model of splicing regulation in which age-dependent changes in H3K27me3/Ac modify the rate of RNAPII elongation and affect processing of ELP1 alternative exon 20.

An interpretable model of pre-mRNA splicing for animal and plant genes.

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Disentangling genetic effects on transcriptional and post-transcriptional gene regulation through integrating exon and intron expression QTLs.

Expression quantitative trait loci (eQTL) studies typically consider exon expression of genes and discard intronic RNA sequencing reads despite their information on RNA metabolism. Here, we quantify genetic effects on exon and intron levels of genes and their ratio in lymphoblastoid cell lines, revealing thousands of cis-QTLs of each type. While genetic effects are often shared between cis-QTL types, 7814 (47%) are not detected as top cis-QTLs at exon levels. We show that exon levels preferentially capture genetic effects on transcriptional regulation, while exon-intron-ratios better detect those on co- and post-transcriptional processes. Considering all cis-QTL types substantially increases (by 71%) the number of colocalizing variants identified by genome-wide association studies (GWAS). It further allows dissecting the potential gene regulatory processes underlying GWAS associations, suggesting comparable contributions by transcriptional (50%) and co- and post-transcriptional regulation (46%) to complex traits. Overall, integrating intronic RNA sequencing reads in eQTL studies expands our understanding of genetic effects on gene regulatory processes.

In vivo genome editing via CRISPR/Cas9-mediated homology-independent targeted integration for Bietti crystalline corneoretinal dystrophy treatment.

Bietti crystalline corneoretinal dystrophy (BCD) is an autosomal recessive chorioretinal degenerative disease without approved therapeutic drugs. It is caused by mutations in CYP4V2 gene, and about 80% of BCD patients carry mutations in exon 7 to 11. Here, we apply CRISPR/Cas9 mediated homology-independent targeted integration (HITI)-based gene editing therapy in HEK293T cells, BCD patient derived iPSCs, and humanized Cyp4v3 mouse model (h-Cyp4v3mut/mut) using two rAAV2/8 vectors via sub-retinal administration. We find that sgRNA-guided Cas9 generates double-strand cleavage on intron 6 of the CYP4V2 gene, and the HITI donor inserts the carried sequence, part of intron 6, exon 7-11, and a stop codon into the DNA break, achieving precise integration, effective transcription and translation both in vitro and in vivo. HITI-based editing restores the viability of iPSC-RPE cells from BCD patient, improves the morphology, number and metabolism of RPE and photoreceptors in h-Cyp4v3mut/mut mice. These results suggest that HITI-based editing could be a promising therapeutic strategy for those BCD patients carrying mutations in exon 7 to 11, and one injection will achieve lifelong effectiveness.

Current status and breakthroughs in treating advanced non-small cell lung cancer with EGFR exon 20 insertion mutations.

Recently, targeted therapy and immunotherapy have emerged as effective treatment options for non-small cell lung cancer (NSCLC). This progress has been facilitated by the rapid development of diagnostic and therapeutic technologies and the continuous research and development of new drugs, leading to a new era in precision medicine for NSCLC. This is a breakthrough for patients with common mutations in the human epidermal growth factor receptor (EGFR) gene in NSCLC. Consequently, the use of targeted drugs has significantly improved survival. Nevertheless, certain rare genetic mutations are referred to as EGFR exon 20 insertion (ex20ins) mutations, which differ in structure from conventional EGFR gene mutations, namely, exon 19 deletion mutations (19-Del) and exon 21 point mutations. Owing to their distinct structural characteristics, patients harboring these EGFR ex20ins mutations are unresponsive to traditional tyrosine kinase inhibitor (TKI) therapy. This particular group of patients did not fall within the scope of their applicability. However, the activating A763_Y764insFQEA mutation elicits a more pronounced response than mutations in the near and far regions of the C-helix immediately following it and should, therefore, be treated differently. Currently, there is a lack of effective treatments for EGFR ex20ins mutations NSCLC. The efficacy of chemotherapy has been relatively favorable, whereas the effectiveness of immunotherapy remains ambiguous owing to inadequate clinical data. In addition, the efficacy of the first- and second-generation targeted drugs remains limited. However, third-generation and novel targeted drugs have proven to be effective. Although novel EGFR-TKIs are expected to treat EGFR ex20ins mutations in patients with NSCLC, they face many challenges. The main focus of this review is on emerging therapies that target NSCLC with EGFR ex20ins and highlight major ongoing clinical trials while also providing an overview of the associated challenges and research advancements in this area.

Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation.

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Assessing the Impact of Novel BRCA1 Exon 11 Variants on Pre-mRNA Splicing.

Our study focused on assessing the effects of three newly identified BRCA1 exon 11 variants (c.1019T>C, c.2363T>G, and c.3192T>C) on breast cancer susceptibility. Using computational predictions and experimental splicing assays, we evaluated their potential as pathogenic mutations. Our in silico analyses suggested that the c.2363T>G and c.3192T>C variants could impact both splicing and protein function, resulting in the V340A and V788G mutations, respectively. We further examined their splicing effects using minigene assays in MCF7 and SKBR3 breast cancer cell lines. Interestingly, we found that the c.2363T>G variant significantly altered splicing patterns in MCF7 cells but not in SKBR3 cells. This finding suggests a potential influence of cellular context on the variant's effects. While attempts to correlate in silico predictions with RNA binding factors were inconclusive, this observation underscores the complexity of splicing regulation. Splicing is governed by various factors, including cellular contexts and protein interactions, making it challenging to predict outcomes accurately. Further research is needed to fully understand the functional consequences of the c.2363T>G variant in breast cancer pathogenesis. Integrating computational predictions with experimental data will provide valuable insights into the role of alternative splicing regulation in different breast cancer types and stages.

Discovery of immunotherapy targets for pediatric solid and brain tumors by exon-level expression.

Immunotherapy with chimeric antigen receptor T cells for pediatric solid and brain tumors is constrained by available targetable antigens. Cancer-specific exons present a promising reservoir of targets; however, these have not been explored and validated systematically in a pan-cancer fashion. To identify cancer specific exon targets, here we analyze 1532 RNA-seq datasets from 16 types of pediatric solid and brain tumors for comparison with normal tissues using a newly developed workflow. We find 2933 exons in 157 genes encoding proteins of the surfaceome or matrisome with high cancer specificity either at the gene (n = 148) or the alternatively spliced isoform (n = 9) level. Expression of selected alternatively spliced targets, including the EDB domain of fibronectin 1, and gene targets, such as COL11A1, are validated in pediatric patient derived xenograft tumors. We generate T cells expressing chimeric antigen receptors specific for the EDB domain or COL11A1 and demonstrate that these have antitumor activity. The full target list, explorable via an interactive web portal ( https://cseminer.stjude.org/ ), provides a rich resource for developing immunotherapy of pediatric solid and brain tumors using gene or AS targets with high expression specificity in cancer.

A splice site variant in MADD affects hormone expression in pancreatic ß cells and pituitary gonadotropes.

MAPK activating death domain (MADD) is a multifunctional protein regulating small GTPases RAB3 and RAB27, MAPK signaling, and cell survival. Polymorphisms in the MADD locus are associated with glycemic traits, but patients with biallelic variants in MADD manifest a complex syndrome affecting nervous, endocrine, exocrine, and hematological systems. We identified a homozygous splice site variant in MADD in 2 siblings with developmental delay, diabetes, congenital hypogonadotropic hypogonadism, and growth hormone deficiency. This variant led to skipping of exon 30 and in-frame deletion of 36 amino acids. To elucidate how this mutation causes pleiotropic endocrine phenotypes, we generated relevant cellular models with deletion of MADD exon 30 (dex30). We observed reduced numbers of ß cells, decreased insulin content, and increased proinsulin-to-insulin ratio in dex30 human embryonic stem cell-derived pancreatic islets. Concordantly, dex30 led to decreased insulin expression in human ß cell line EndoC-ßH1. Furthermore, dex30 resulted in decreased luteinizing hormone expression in mouse pituitary gonadotrope cell line LßT2 but did not affect ontogeny of stem cell-derived GnRH neurons. Protein-protein interactions of wild-type and dex30 MADD revealed changes affecting multiple signaling pathways, while the GDP/GTP exchange activity of dex30 MADD remained intact. Our results suggest MADD-specific processes regulate hormone expression in pancreatic ß cells and pituitary gonadotropes.

CRISPR-dCas13d-based deep screening of proximal and distal splicing-regulatory elements.

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is hampered by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically inactive CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identifies not only known SREs but also a previously unknown distal intronic SRE, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

Transcriptome- and proteome-wide effects of a circular RNA encompassing four early exons of the spinal muscular atrophy genes.

Spinal muscular atrophy (SMA) genes, SMN1 and SMN2 (hereinafter referred to as SMN1/2), produce multiple circular RNAs (circRNAs), including C2A-2B-3-4 that encompasses early exons 2A, 2B, 3 and 4. C2A-2B-3-4 is a universally and abundantly expressed circRNA of SMN1/2. Here we report the transcriptome- and proteome-wide effects of overexpression of C2A-2B-3-4 in inducible HEK293 cells. Our RNA-Seq analysis revealed altered expression of ~ 15% genes (4172 genes) by C2A-2B-3-4. About half of the affected genes by C2A-2B-3-4 remained unaffected by L2A-2B-3-4, a linear transcript encompassing exons 2A, 2B, 3 and 4 of SMN1/2. These findings underscore the unique role of the structural context of C2A-2B-3-4 in gene regulation. A surprisingly high number of upregulated genes by C2A-2B-3-4 were located on chromosomes 4 and 7, whereas many of the downregulated genes were located on chromosomes 10 and X. Supporting a cross-regulation of SMN1/2 transcripts, C2A-2B-3-4 and L2A-2B-3-4 upregulated and downregulated SMN1/2 mRNAs, respectively. Proteome analysis revealed 61 upregulated and 57 downregulated proteins by C2A-2B-3-4 with very limited overlap with those affected by L2A-2B-3-4. Independent validations confirmed the effect of C2A-2B-3-4 on expression of genes associated with chromatin remodeling, transcription, spliceosome function, ribosome biogenesis, lipid metabolism, cytoskeletal formation, cell proliferation and neuromuscular junction formation. Our findings reveal a broad role of C2A-2B-3-4, and expands our understanding of functions of SMN1/2 genes.

Exon-junction complex association with stalled ribosomes and slow translation-independent disassembly.

Exon junction complexes are deposited at exon-exon junctions during splicing. They are primarily known to activate non-sense mediated degradation of transcripts harbouring premature stop codons before the last intron. According to a popular model, exon-junction complexes accompany mRNAs to the cytoplasm where the first translating ribosome pushes them out. However, they are also removed by uncharacterized, translation-independent mechanisms. Little is known about kinetic and transcript specificity of these processes. Here we tag core subunits of exon-junction complexes with complementary split nanoluciferase fragments to obtain sensitive and quantitative assays for complex formation. Unexpectedly, exon-junction complexes form large stable mRNPs containing stalled ribosomes. Complex assembly and disassembly rates are determined after an arrest in transcription and/or translation. 85% of newly deposited exon-junction complexes are disassembled by a translation-dependent mechanism. However as this process is much faster than the translation-independent one, only 30% of the exon-junction complexes present in cells at steady state require translation for disassembly. Deep RNA sequencing shows a bias of exon-junction complex bound transcripts towards microtubule and centrosome coding ones and demonstrate that the lifetimes of exon-junction complexes are transcript-specific. This study provides a dynamic vision of exon-junction complexes and uncovers their unexpected stable association with ribosomes.

Copper-binding proteins and exonic splicing enhancers and silencers.

Eukaryotic DNA codes not only for proteins but contains a wealth of information required for accurate splicing of messenger RNA precursors and inclusion of constitutively or alternatively spliced exons in mature transcripts. This "auxiliary" splicing code has been characterized as exonic splicing enhancers and silencers (ESE and ESS). The exact interplay between protein and splicing codes is, however, poorly understood. Here, we show that exons encoding copper-coordinating amino acids in human cuproproteins lack ESEs and/or have an excess of ESSs, yet RNA sequencing and expressed sequence tags data show that they are more efficiently included in mature transcripts by the splicing machinery than average exons. Their largely constitutive inclusion in messenger RNA is facilitated by stronger splice sites, including polypyrimidine tracts, consistent with an important role of the surrounding intron architecture in ensuring high expression of metal-binding residues during evolution. ESE/ESS profiles of codons and entire exons that code for copper-coordinating residues were very similar to those encoding residues that coordinate zinc but markedly different from those that coordinate calcium. Together, these results reveal how the traditional and auxiliary splicing motifs responded to constraints of metal coordination in proteins.

Neuregulin 2 Is a Candidate Gene for Autism Spectrum Disorder.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with heterogeneous and complex genetic underpinnings. Our previous microarray gene expression profiling identified significantly different neuregulin-2 gene (NRG2) expression between ASD patients and controls. Thus, we aimed to clarify whether NRG2 is a candidate gene associated with ASD. The study consisted of two stages. First, we used real-time quantitative PCR in 20 ASDs and 20 controls to confirm the microarray gene expression profiling results. The average NRG2 gene expression level in patients with ASD (3.23 ± 2.80) was significantly lower than that in the controls (9.27 ± 4.78, p < 0.001). Next, we conducted resequencing of all the exons of NRG2 in a sample of 349 individuals with ASD, aiming to identify variants of the NRG2 associated with ASD. We identified three variants, including two single nucleotide variants (SNVs), IVS3 + 13A > G (rs889022) and IVS10 + 32T > A (rs182642591), and one small deletion at exon 11 of NRG2 (delGCCCGG, rs933769137). Using data from the Taiwan Biobank as the controls, we found no significant differences in allele frequencies of rs889022 and rs182642591 between two groups. However, there is a significant difference in the genotype and allele frequency distribution of rs933769137 between ASDs and controls (p < 0.0001). The small deletion is located in the EGF-like domain at the C-terminal of the NRG2 precursor protein. Our findings suggest that NRG2 might be a susceptibility gene for ASD.

Comparative Genome-Wide Alternative Splicing Analysis between Preadipocytes and Adipocytes.

Alternative splicing is a ubiquitous regulatory mechanism in gene expression that allows a single gene to generate multiple messenger RNAs (mRNAs). Adipocyte development is regulated by many processes, and recent studies have found that splicing factors also play an important role in adipogenic development. In the present study, we further investigated the differences in selective shearing during different periods of adipocyte differentiation. We identified five alternative splicing types including skipped exon, mutually exclusive exon, Alternative 5' splice site, Alternative 3' splice site, and Retained intron, with skipped exons being the most abundant type of selective shearing. In total, 641 differentially expressed selective shearing genes were obtained, enriched in 279 pathways, from which we selected and verified the accuracy of the sequencing results. Overall, RNA-seq revealed changes in the splicing and expression levels of these new candidate genes between precursor adipocytes and adipocytes, suggesting that they may be involved in adipocyte generation and differentiation.

The debranching enzyme Dbr1 regulates lariat turnover and intron splicing.

The majority of genic transcription is intronic. Introns are removed by splicing as branched lariat RNAs which require rapid recycling. The branch site is recognized during splicing catalysis and later debranched by Dbr1 in the rate-limiting step of lariat turnover. Through generation of a viable DBR1 knockout cell line, we find the predominantly nuclear Dbr1 enzyme to encode the sole debranching activity in human cells. Dbr1 preferentially debranches substrates that contain canonical U2 binding motifs, suggesting that branchsites discovered through sequencing do not necessarily represent those favored by the spliceosome. We find that Dbr1 also exhibits specificity for particular 5' splice site sequences. We identify Dbr1 interactors through co-immunoprecipitation mass spectrometry. We present a mechanistic model for Dbr1 recruitment to the branchpoint through the intron-binding protein AQR. In addition to a 20-fold increase in lariats, Dbr1 depletion increases exon skipping. Using ADAR fusions to timestamp lariats, we demonstrate a defect in spliceosome recycling. In the absence of Dbr1, spliceosomal components remain associated with the lariat for a longer period of time. As splicing is co-transcriptional, slower recycling increases the likelihood that downstream exons will be available for exon skipping.

On a kneading theory for gene-splicing.

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

Alternative splicing for tuneable expression of protein subunits at desired ratios.

The controlled expression of two or more proteins at a defined and stable ratio remains a substantial challenge, particularly in the bi- and multispecific antibody field. Achieving an optimal ratio of protein subunits can facilitate the assembly of multimeric proteins with high efficiency and minimize the production of by-products. In this study, we propose a solution based on alternative splicing, enabling the expression of a tunable and predefined ratio of two distinct polypeptide chains from the same pre-mRNA under the control of a single promoter. The pre-mRNA used in this study contains two open reading frames situated on separate exons. The first exon is flanked by two copies of the chicken troponin intron 4 (cTNT-I4) and is susceptible to excision from the pre-mRNA by means of alternative splicing. This specific design enables the modulation of the splice ratio by adjusting the strength of the splice acceptor. To illustrate this approach, we developed constructs expressing varying ratios of GFP and dsRED and extended their application to multimeric proteins such as monoclonal antibodies, achieving industrially relevant expression levels (>1 g/L) in a 14-day fed-batch process. The stability of the splice ratio was confirmed by droplet digital PCR in a stable pool cultivated over a 28-day period, while product quality was assessed via intact mass analysis, demonstrating absence of product-related impurities resulting from undesired splice events. Furthermore, we showcased the versatility of the construct by expressing two subunits of a bispecific antibody of the BEAT® type, which contains three distinct subunits in total.