Genome-scale exon perturbation screens uncover exons critical for cell fitness.

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

Meet the authors: Mei-Sheng Xiao, Arun Prasath Damodaran, and Thomas Gonatopoulos-Pournatzis.

We talk to corresponding author Thomas Gonatopoulos-Pournatzis and co-first authors Arun Prasath Damodaran and Mei-Sheng Xiao about their paper "Genome-scale exon perturbation screens uncover exons critical for cell fitness" (in this issue of Molecular Cell) and get insights into their findings, career trajectories, and future directions in the pre-mRNA processing field.

CCAR1 promotes DNA repair via alternative splicing.

DNA repair is directly performed by hundreds of core factors and indirectly regulated by thousands of others. We massively expanded a CRISPR inhibition and Cas9-editing screening system to discover factors indirectly modulating homology-directed repair (HDR) in the context of ∼18,000 individual gene knockdowns. We focused on CCAR1, a poorly understood gene that we found the depletion of reduced both HDR and interstrand crosslink repair, phenocopying the loss of the Fanconi anemia pathway. CCAR1 loss abrogated FANCA protein without substantial reduction in the level of its mRNA or that of other FA genes. We instead found that CCAR1 prevents inclusion of a poison exon in FANCA. Transcriptomic analysis revealed that the CCAR1 splicing modulatory activity is not limited to FANCA, and it instead regulates widespread changes in alternative splicing that would damage coding sequences in mouse and human cells. CCAR1 therefore has an unanticipated function as a splicing fidelity factor.

hnRNPM protects against the dsRNA-mediated interferon response by repressing LINE-associated cryptic splicing.

RNA splicing is pivotal in post-transcriptional gene regulation, yet the exponential expansion of intron length in humans poses a challenge for accurate splicing. Here, we identify hnRNPM as an essential RNA-binding protein that suppresses cryptic splicing through binding to deep introns, maintaining human transcriptome integrity. Long interspersed nuclear elements (LINEs) in introns harbor numerous pseudo splice sites. hnRNPM preferentially binds at intronic LINEs to repress pseudo splice site usage for cryptic splicing. Remarkably, cryptic exons can generate long dsRNAs through base-pairing of inverted ALU transposable elements interspersed among LINEs and consequently trigger an interferon response, a well-known antiviral defense mechanism. Significantly, hnRNPM-deficient tumors show upregulated interferon-associated pathways and elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity by repressing cryptic splicing and suggest that targeting hnRNPM in tumors may be used to trigger an inflammatory immune response, thereby boosting cancer surveillance.

Efficient, specific, and combinatorial control of endogenous exon splicing with dCasRx-RBM25.

Efficient targeted control of splicing is a major goal of functional genomics and therapeutic applications. Guide (g)RNA-directed, deactivated (d)Cas CRISPR enzymes fused to splicing effectors represent a promising strategy due to the flexibility of these systems. However, efficient, specific, and generalizable activation of endogenous exons using this approach has not been previously reported. By screening over 300 dCasRx-splicing factor fusion proteins tethered to splicing reporters, we identify dCasRx-RBM25 as a potent activator of exons. Moreover, dCasRx-RBM25 efficiently activates the splicing of ∼90% of targeted endogenous alternative exons and displays high on-target specificity. Using gRNA arrays for combinatorial targeting, we demonstrate that dCasRx-RBM25 enables multiplexed activation and repression of exons. Using this feature, the targeting of neural-regulated exons in Ptpb1 and Puf60 in embryonic stem cells reveals combinatorial effects on downstream alternative splicing events controlled by these factors. Collectively, our results enable versatile, combinatorial exon-resolution functional assays and splicing-directed therapeutic applications.

The splicing factor CCAR1 regulates the Fanconi anemia/BRCA pathway.

The twenty-three Fanconi anemia (FA) proteins cooperate in the FA/BRCA pathway to repair DNA interstrand cross-links (ICLs). The cell division cycle and apoptosis regulator 1 (CCAR1) protein is also a regulator of ICL repair, though its possible function in the FA/BRCA pathway remains unknown. Here, we demonstrate that CCAR1 plays a unique upstream role in the FA/BRCA pathway and is required for FANCA protein expression in human cells. Interestingly, CCAR1 co-immunoprecipitates with FANCA pre-mRNA and is required for FANCA mRNA processing. Loss of CCAR1 results in retention of a poison exon in the FANCA transcript, thereby leading to reduced FANCA protein expression. A unique domain of CCAR1, the EF hand domain, is required for interaction with the U2AF heterodimer of the spliceosome and for excision of the poison exon. Taken together, CCAR1 is a splicing modulator required for normal splicing of the FANCA mRNA and other mRNAs involved in various cellular pathways.

4.5SH RNA counteracts deleterious exonization of SINE B1 in mice.

4.5SH RNA is a highly abundant, small rodent-specific noncoding RNA that localizes to nuclear speckles enriched in pre-mRNA-splicing regulators. To investigate the physiological functions of 4.5SH RNA, we have created mutant mice that lack the expression of 4.5SH RNA. The mutant mice exhibited embryonic lethality, suggesting that 4.5SH RNA is an essential species-specific noncoding RNA in mice. RNA-sequencing analyses revealed that 4.5SH RNA protects the transcriptome from abnormal exonizations of the antisense insertions of the retrotransposon SINE B1 (asB1), which would otherwise introduce deleterious premature stop codons or frameshift mutations. Mechanistically, 4.5SH RNA base pairs with complementary asB1-containing exons via the target recognition region and recruits effector proteins including Hnrnpm via its 5' stem loop region. The modular organization of 4.5SH RNA allows us to engineer a programmable splicing regulator to induce the skipping of target exons of interest. Our results also suggest the general existence of splicing regulatory noncoding RNAs.

Structural insights into the cross-exon to cross-intron spliceosome switch.

Early spliceosome assembly can occur through an intron-defined pathway, whereby U1 and U2 small nuclear ribonucleoprotein particles (snRNPs) assemble across the intron1. Alternatively, it can occur through an exon-defined pathway2-5, whereby U2 binds the branch site located upstream of the defined exon and U1 snRNP interacts with the 5' splice site located directly downstream of it. The U4/U6.U5 tri-snRNP subsequently binds to produce a cross-intron (CI) or cross-exon (CE) pre-B complex, which is then converted to the spliceosomal B complex6,7. Exon definition promotes the splicing of upstream introns2,8,9 and plays a key part in alternative splicing regulation10-16. However, the three-dimensional structure of exon-defined spliceosomal complexes and the molecular mechanism of the conversion from a CE-organized to a CI-organized spliceosome, a pre-requisite for splicing catalysis, remain poorly understood. Here cryo-electron microscopy analyses of human CE pre-B complex and B-like complexes reveal extensive structural similarities with their CI counterparts. The results indicate that the CE and CI spliceosome assembly pathways converge already at the pre-B stage. Add-back experiments using purified CE pre-B complexes, coupled with cryo-electron microscopy, elucidate the order of the extensive remodelling events that accompany the formation of B complexes and B-like complexes. The molecular triggers and roles of B-specific proteins in these rearrangements are also identified. We show that CE pre-B complexes can productively bind in trans to a U1 snRNP-bound 5' splice site. Together, our studies provide new mechanistic insights into the CE to CI switch during spliceosome assembly and its effect on pre-mRNA splice site pairing at this stage.

Autonomous transposons tune their sequences to ensure somatic suppression.

Transposable elements (TEs) are a major constituent of human genes, occupying approximately half of the intronic space. During pre-messenger RNA synthesis, intronic TEs are transcribed along with their host genes but rarely contribute to the final mRNA product because they are spliced out together with the intron and rapidly degraded. Paradoxically, TEs are an abundant source of RNA-processing signals through which they can create new introns1, and also functional2 or non-functional chimeric transcripts3. The rarity of these events implies the existence of a resilient splicing code that is able to suppress TE exonization without compromising host pre-mRNA processing. Here we show that SAFB proteins protect genome integrity by preventing retrotransposition of L1 elements while maintaining splicing integrity, via prevention of the exonization of previously integrated TEs. This unique dual role is possible because of L1's conserved adenosine-rich coding sequences that are bound by SAFB proteins. The suppressive activity of SAFB extends to tissue-specific, giant protein-coding cassette exons, nested genes and Tigger DNA transposons. Moreover, SAFB also suppresses LTR/ERV elements in species in which they are still active, such as mice and flies. A significant subset of splicing events suppressed by SAFB in somatic cells are activated in the testis, coinciding with low SAFB expression in postmeiotic spermatids. Reminiscent of the division of labour between innate and adaptive immune systems that fight external pathogens, our results uncover SAFB proteins as an RNA-based, pattern-guided, non-adaptive defence system against TEs in the soma, complementing the RNA-based, adaptive Piwi-interacting RNA pathway of the germline.

On the genetic basis of tail-loss evolution in humans and apes.

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.

Antisense oligonucleotide therapeutic approach for Timothy syndrome.

Timothy syndrome (TS) is a severe, multisystem disorder characterized by autism, epilepsy, long-QT syndrome and other neuropsychiatric conditions1. TS type 1 (TS1) is caused by a gain-of-function variant in the alternatively spliced and developmentally enriched CACNA1C exon 8A, as opposed to its counterpart exon 8. We previously uncovered several phenotypes in neurons derived from patients with TS1, including delayed channel inactivation, prolonged depolarization-induced calcium rise, impaired interneuron migration, activity-dependent dendrite retraction and an unanticipated persistent expression of exon 8A2-6. We reasoned that switching CACNA1C exon utilization from 8A to 8 would represent a potential therapeutic strategy. Here we developed antisense oligonucleotides (ASOs) to effectively decrease the inclusion of exon 8A in human cells both in vitro and, following transplantation, in vivo. We discovered that the ASO-mediated switch from exon 8A to 8 robustly rescued defects in patient-derived cortical organoids and migration in forebrain assembloids. Leveraging a transplantation platform previously developed7, we found that a single intrathecal ASO administration rescued calcium changes and in vivo dendrite retraction of patient neurons, suggesting that suppression of CACNA1C exon 8A expression is a potential treatment for TS1. Broadly, these experiments illustrate how a multilevel, in vivo and in vitro stem cell model-based approach can identify strategies to reverse disease-relevant neural pathophysiology.

Cryo-EM analyses of dimerized spliceosomes provide new insights into the functions of B complex proteins.

The B complex is a key intermediate stage of spliceosome assembly. To improve the structural resolution of monomeric, human spliceosomal B (hB) complexes and thereby generate a more comprehensive hB molecular model, we determined the cryo-EM structure of B complex dimers formed in the presence of ATP γ S. The enhanced resolution of these complexes allows a finer molecular dissection of how the 5' splice site (5'ss) is recognized in hB, and new insights into molecular interactions of FBP21, SNU23 and PRP38 with the U6/5'ss helix and with each other. It also reveals that SMU1 and RED are present as a heterotetrameric complex and are located at the interface of the B dimer protomers. We further show that MFAP1 and UBL5 form a 5' exon binding channel in hB, and elucidate the molecular contacts stabilizing the 5' exon at this stage. Our studies thus yield more accurate models of protein and RNA components of hB complexes. They further allow the localization of additional proteins and protein domains (such as SF3B6, BUD31 and TCERG1) whose position was not previously known, thereby uncovering new functions for B-specific and other hB proteins during pre-mRNA splicing.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia.

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

Functional characterization of 2,832 JAG1 variants supports reclassification for Alagille syndrome and improves guidance for clinical variant interpretation.

Pathogenic variants in the JAG1 gene are a primary cause of the multi-system disorder Alagille syndrome. Although variant detection rates are high for this disease, there is uncertainty associated with the classification of missense variants that leads to reduced diagnostic yield. Consequently, up to 85% of reported JAG1 missense variants have uncertain or conflicting classifications. We generated a library of 2,832 JAG1 nucleotide variants within exons 1-7, a region with a high number of reported missense variants, and designed a high-throughput assay to measure JAG1 membrane expression, a requirement for normal function. After calibration using a set of 175 known or predicted pathogenic and benign variants included within the variant library, 486 variants were characterized as functionally abnormal (n = 277 abnormal and n = 209 likely abnormal), of which 439 (90.3%) were missense. We identified divergent membrane expression occurring at specific residues, indicating that loss of the wild-type residue itself does not drive pathogenicity, a finding supported by structural modeling data and with broad implications for clinical variant classification both for Alagille syndrome and globally across other disease genes. Of 144 uncertain variants reported in patients undergoing clinical or research testing, 27 had functionally abnormal membrane expression, and inclusion of our data resulted in the reclassification of 26 to likely pathogenic. Functional evidence augments the classification of genomic variants, reducing uncertainty and improving diagnostics. Inclusion of this repository of functional evidence during JAG1 variant reclassification will significantly affect resolution of variant pathogenicity, making a critical impact on the molecular diagnosis of Alagille syndrome.

De novo missense variants in exon 9 of SEPHS1 cause a neurodevelopmental condition with developmental delay, poor growth, hypotonia, and dysmorphic features.

Selenophosphate synthetase (SEPHS) plays an essential role in selenium metabolism. Two mammalian SEPHS paralogues, SEPHS1 and SEPHS2, share high sequence identity and structural homology with SEPHS. Here, we report nine individuals from eight families with developmental delay, growth and feeding problems, hypotonia, and dysmorphic features, all with heterozygous missense variants in SEPHS1. Eight of these individuals had a recurrent variant at amino acid position 371 of SEPHS1 (p.Arg371Trp, p.Arg371Gln, and p.Arg371Gly); seven of these variants were known to be de novo. Structural modeling and biochemical assays were used to understand the effect of these variants on SEPHS1 function. We found that a variant at residue Trp352 results in local structural changes of the C-terminal region of SEPHS1 that decrease the overall thermal stability of the enzyme. In contrast, variants of a solvent-exposed residue Arg371 do not impact enzyme stability and folding but could modulate direct protein-protein interactions of SEPSH1 with cellular factors in promoting cell proliferation and development. In neuronal SH-SY5Y cells, we assessed the impact of SEPHS1 variants on cell proliferation and ROS production and investigated the mRNA expression levels of genes encoding stress-related selenoproteins. Our findings provided evidence that the identified SEPHS1 variants enhance cell proliferation by modulating ROS homeostasis. Our study supports the hypothesis that SEPHS1 plays a critical role during human development and provides a basis for further investigation into the molecular mechanisms employed by SEPHS1. Furthermore, our data suggest that variants in SEPHS1 are associated with a neurodevelopmental disorder.

Systematic identification and characterization of exon-intron circRNAs.

Exon-intron circRNAs (EIciRNAs) are a circRNA subclass with retained introns. Global features of EIciRNAs remain largely unexplored, mainly owing to the lack of bioinformatic tools. The regulation of intron retention (IR) in EIciRNAs and the associated functionality also require further investigation. We developed a framework, FEICP, which efficiently detected EIciRNAs from high-throughput sequencing (HTS) data. EIciRNAs are distinct from exonic circRNAs (EcircRNAs) in aspects such as with larger length, localization in the nucleus, high tissue specificity, and enrichment mostly in the brain. Deep learning analyses revealed that compared with regular introns, the retained introns of circRNAs (CIRs) are shorter in length, have weaker splice site strength, and have higher GC content. Compared with retained introns in linear RNAs (LIRs), CIRs are more likely to form secondary structures and show greater sequence conservation. CIRs are closer to the 5'-end, whereas LIRs are closer to the 3'-end of transcripts. EIciRNA-generating genes are more actively transcribed and associated with epigenetic marks of gene activation. Computational analyses and genome-wide CRISPR screening revealed that SRSF1 binds to CIRs and inhibits the biogenesis of most EIciRNAs. SRSF1 regulates the biogenesis of EIciLIMK1, which enhances the expression of LIMK1 in cis to boost neuronal differentiation, exemplifying EIciRNA physiological function. Overall, our study has developed the FEICP pipeline to identify EIciRNAs from HTS data, and reveals multiple features of CIRs and EIciRNAs. SRSF1 has been identified to regulate EIciRNA biogenesis. EIciRNAs and the regulation of EIciRNA biogenesis play critical roles in neuronal differentiation.

Reassessing the exon-foldon correspondence using frustration analysis.

Protein folding and evolution are intimately linked phenomena. Here, we revisit the concept of exons as potential protein folding modules across a set of 38 abundant and conserved protein families. Taking advantage of genomic exon-intron organization and extensive protein sequence data, we explore exon boundary conservation and assess the foldon-like behavior of exons using energy landscape theoretic measurements. We found deviations in the exon size distribution from exponential decay indicating selection in evolution. We show that when taken together there is a pronounced tendency to independent foldability for segments corresponding to the more conserved exons, supporting the idea of exon-foldon correspondence. While 45% of the families follow this general trend when analyzed individually, there are some families for which other stronger functional determinants, such as preserving frustrated active sites, may be acting. We further develop a systematic partitioning of protein domains using exon boundary hotspots, showing that minimal common exons correspond with uninterrupted alpha and/or beta elements for the majority of the families but not for all of them.

Hox gene-specific cellular targeting using split intein Trojan exons.

The Trojan exon method, which makes use of intronically inserted T2A-Gal4 cassettes, has been widely used in Drosophila to create thousands of gene-specific Gal4 driver lines. These dual-purpose lines provide genetic access to specific cell types based on their expression of a native gene while simultaneously mutating one allele of the gene to enable loss-of-function analysis in homozygous animals. While this dual use is often an advantage, the truncation mutations produced by Trojan exons are sometimes deleterious in heterozygotes, perhaps by creating translation products with dominant negative effects. Such mutagenic effects can cause developmental lethality as has been observed with genes encoding essential transcription factors. Given the importance of transcription factors in specifying cell type, alternative techniques for generating specific Gal4 lines that target them are required. Here, we introduce a modified Trojan exon method that retains the targeting fidelity and plug-and-play modularity of the original method but mitigates its mutagenic effects by exploiting the self-splicing capabilities of split inteins. "Split Intein Trojan exons" (siTrojans) ensure that the two truncation products generated from the interrupted allele of the native gene are trans-spliced to create a full-length native protein. We demonstrate the efficacy of siTrojans by generating a comprehensive toolkit of Gal4 and Split Gal4 lines for the segmentally expressed Hox transcription factors and illustrate their use in neural circuit mapping by targeting neurons according to their position along the anterior-posterior axis. Both the method and the Hox gene-specific toolkit introduced here should be broadly useful.

The immunoglobulin J chain is an evolutionarily co-opted chemokine.

The joining (J) chain regulates polymerization of multimeric Immunoglobulin(Ig)M and IgA, forming a disulfide bond to the C termini of their Ig heavy chains, and it controls IgM/IgA transport across mucosal epithelia. Like Ig itself and human-like adaptive immunity, J chain emerged in jawed vertebrates (gnathostomes), but its origin has remained mysterious since its discovery over 50 y ago. Here, we show unexpectedly that J chain is a member of the CXCL chemokine family. The J chain gene (JCHAIN) is linked to clustered CXCL chemokine loci in all gnathostomes except actinopterygians that lost JCHAIN. JCHAIN and most CXCL genes have four exons with the same intron phases, including the same cleavage site for the signal peptide/mature protein. The second exon of both genes encodes a CXC motif at the same position, and the lengths of exons 1 to 3 are similar. No other gene in the human secretome shares all of these characteristics. In contrast, intrachain disulfide bonds of the two proteins are completely different, likely due to modifications in J chain to direct Ig polymerization and mucosal transport. Crystal structures of CXCL8 and J chain share a conserved beta-strand core but diverge otherwise due to different intrachain disulfide bonds and extension of the J chain C terminus. Identification of this ancestral affiliation between J chain and CXCL chemokines addresses an age-old problem in immunology.

Splice modulation strategy applied to deep intronic variants in COL7A1 causing recessive dystrophic epidermolysis bullosa.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare and most often severe genetic disease characterized by recurrent blistering and erosions of the skin and mucous membranes after minor trauma, leading to major local and systemic complications. The disease is caused by loss-of-function variants in COL7A1 encoding type VII collagen (C7), the main component of anchoring fibrils, which form attachment structures stabilizing the cutaneous basement membrane zone. Alterations in C7 protein structure and/or expression lead to abnormal, rare or absent anchoring fibrils resulting in loss of dermal-epidermal adherence and skin blistering. To date, more than 1,200 distinct COL7A1 deleterious variants have been reported and 19% are splice variants. Here, we describe two RDEB patients for whom we identified two pathogenic deep intronic pathogenic variants in COL7A1. One of these variants (c.7795-97C > G) promotes the inclusion of a pseudoexon between exons 104 and 105 in the COL7A1 transcript, while the other causes partial or complete retention of intron 51. We used antisense oligonucleotide (ASO) mediated exon skipping to correct these aberrant splicing events in vitro. This led to increased normal mRNA splicing above 94% and restoration of C7 protein expression at a level (up to 56%) that should be sufficient to reverse the phenotype. This first report of exon skipping applied to counteract deep intronic variants in COL7A1 represents a promising therapeutic strategy for personalized medicine directed at patients with intronic variants at a distance of consensus splice sites.