Search

Loss of nuclear TDP-43 is a hallmark of neurodegeneration in TDP-43 proteinopathies, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 mislocalization results in cryptic splicing and polyadenylation of pre-messenger RNAs (pre-mRNAs) encoding stathmin-2 (also known as SCG10), a protein that is required for axonal regeneration. We found that TDP-43 binding to a GU-rich region sterically blocked recognition of the cryptic 3 splice site in STMN2 pre-mRNA. Targeting dCasRx or antisense oligonucleotides (ASOs) suppressed cryptic splicing, which restored axonal regeneration and stathmin-2-dependent lysosome trafficking in TDP-43-deficient human motor neurons. In mice that were gene-edited to contain human STMN2 cryptic splice-polyadenylation sequences, ASO injection into cerebral spinal fluid successfully corrected Stmn2 pre-mRNA misprocessing and restored stathmin-2 expression levels independently of TDP-43 binding.

The major spliceosome mediates pre-mRNA splicing by recognizing the highly conserved sequences at the 5' and 3' splice sites and the branch point. More than 150 proteins participate in the splicing process and are organized in the spliceosomal A, B, and C complexes. FRA10AC1 is a peripheral protein of the spliceosomal C complex and its ortholog in the green alga facilitates recognition or interaction with splice sites. We identified biallelic pathogenic variants in FRA10AC1 in five individuals from three consanguineous families. The two unrelated Patients 1 and 2 with loss-of-function variants showed developmental delay, intellectual disability, and no speech, while three siblings with the c.494_496delAAG (p.Glu165del) variant had borderline to mild intellectual disability. All patients had microcephaly, hypoplasia or agenesis of the corpus callosum, growth retardation, and craniofacial dysmorphism. FRA10AC1 transcripts and proteins were drastically reduced or absent in fibroblasts of Patients 1 and 2. In a heterologous expression system, the p.Glu165del variant impacts intrinsic stability of FRA10AC1 but does not affect its nuclear localization. By co-immunoprecipitation, we found ectopically expressed HA-FRA10AC1 in complex with endogenous DGCR14, another component of the spliceosomal C complex, while the splice factors CHERP, NKAP, RED, and SF3B2 could not be co-immunoprecipitated. Using an in vitro splicing reporter assay, we did not obtain evidence for FRA10AC1 deficiency to suppress missplicing events caused by mutations in the highly conserved dinucleotides of 5' and 3' splice sites in an in vitro splicing assay in patient-derived fibroblasts. Our data highlight the importance of specific peripheral spliceosomal C complex proteins for neurodevelopment. It remains possible that FRA10AC1 may have other and/or additional cellular functions, such as coupling of transcription and splicing reactions.

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

The spliceosome undergoes extensive rearrangements as it assembles onto precursor messenger RNAs. In the earliest assembly step, U1snRNA identifies the 5 splice site. However, U1snRNA leaves the spliceosome relatively early in assembly, and 5 splice site identity is subsequently maintained through interactions with U6snRNA, protein factor PRP8, and other components during the rearrangements that build the catalytic site. Using a forward genetic screen in Caenorhabditis elegans, we have identified suppressors of a locomotion defect caused by a 5ss mutation. Here we report three new suppressor alleles from this screen, two in PRP8 and one in SNRNP200/BRR2. mRNASeq studies of these suppressor strains indicate that they also affect specific native alternative 5ss, especially for suppressor PRP8 D1549N. A strong suppressor at the unstructured N-terminus of SNRNP200, N18K, indicates a novel role for this region. By examining distinct changes in the splicing of native genes, examining double mutants between suppressors, comparing these new suppressors to previously identified splicing suppressors from yeast, and mapping conserved suppressor residues onto cryoEM structural models of assembling human spliceosomes, we conclude that there are multiple interactions at multiple stages in spliceosome assembly responsible for maintaining the initial 5ss identified by U1snRNA for entry into the catalytic core.

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Hereditary Diffuse Gastric Cancer (HDGC) syndrome is associated with CDH1 germline likely pathogenic/pathogenic variants. Carriers of CDH1 germline likely pathogenic/pathogenic variants are predisposed to diffuse gastric cancer and lobular breast cancer. This study aims to classify the CDH1 c.[715G>A] missense variant identified in a diffuse gastric cancer prone family by performing splicing studies. RT-PCR and subsequent cloning experiments were performed to investigate whether this variant completely disrupts normal splicing. This variant preferentially abolishes normal splicing through activation of a cryptic 3' acceptor splice site within exon 6 of CDH1, presumably leading to a premature protein truncation within first extracellular domain repeat of E-cadherin protein. Our results contributed to evidence necessary to resolve pathogenicity classification of this variant, indicating that this variant is to be classified as pathogenic.

Trafficking protein particle (TRAPP) complexes, which include the TRAPPC4 protein, regulate membrane trafficking between lipid organelles in a process termed vesicular tethering. TRAPPC4 was recently implicated in a recessive neurodevelopmental condition in four unrelated families due to a shared c.454+3A>G splice variant. Here, we report 23 patients from 17 independent families with an early-infantile-onset neurodegenerative presentation, where we also identified the homozygous variant hg38:11:119020256 A>G (NM_016146.5:c.454+3A>G) in TRAPPC4 through exome or genome sequencing. No other clinically relevant TRAPPC4 variants were identified among any of over 10,000 patients with neurodevelopmental conditions. We found the carrier frequency of TRAPPC4 c.454+3A>G was 2.4-5.4 per 10,000 healthy individuals. Affected individuals with the homozygous TRAPPC4 c.454+3A>G variant showed profound psychomotor delay, developmental regression, early-onset epilepsy, microcephaly and progressive spastic tetraplegia. Based upon RNA sequencing, the variant resulted in partial exon 3 skipping and generation of an aberrant transcript owing to use of a downstream cryptic splice donor site, predicting a premature stop codon and nonsense mediated decay. These data confirm the pathogenicity of the TRAPPC4 c.454+3A>G variant, and refine the clinical presentation of TRAPPC4-related encephalopathy.

MotivationCircular RNA (circRNA) is a novel class of long non-coding RNAs that have been broadly discovered in the eukaryotic transcriptome. The circular structure arises from a non-canonical splicing process, where the donor site backspliced to an upstream acceptor site. These circRNA sequences are conserved across species. More importantly, rising evidence suggests their vital roles in gene regulation and association with diseases. As the fundamental effort toward elucidating their functions and mechanisms, several computational methods have been proposed to predict the circular structure from the primary sequence. Recently, advanced computational methods leverage deep learning to capture the relevant patterns from RNA sequences and model their interactions to facilitate the prediction. However, these methods fail to fully explore positional information of splice junctions and their deep interaction.ResultsWe present a robust end-to-end framework, Junction Encoder with Deep Interaction (JEDI), for circRNA prediction using only nucleotide sequences. JEDI first leverages the attention mechanism to encode each junction site based on deep bidirectional recurrent neural networks and then presents the novel cross-attention layer to model deep interaction among these sites for backsplicing. Finally, JEDI can not only predict circRNAs but also interpret relationships among splice sites to discover backsplicing hotspots within a gene region. Experiments demonstrate JEDI significantly outperforms state-of-the-art approaches in circRNA prediction on both isoform level and gene level. Moreover, JEDI also shows promising results on zero-shot backsplicing discovery, where none of the existing approaches can achieve.Availability and implementationThe implementation of our framework is available at https://github.com/hallogameboy/JEDI.Supplementary informationSupplementary data are available at Bioinformatics online.

MotivationCircular RNA (circRNA) is a novel class of long non-coding RNAs that have been broadly discovered in the eukaryotic transcriptome. The circular structure arises from a non-canonical splicing process, where the donor site backspliced to an upstream acceptor site. These circRNA sequences are conserved across species. More importantly, rising evidence suggests their vital roles in gene regulation and association with diseases. As the fundamental effort toward elucidating their functions and mechanisms, several computational methods have been proposed to predict the circular structure from the primary sequence. Recently, advanced computational methods leverage deep learning to capture the relevant patterns from RNA sequences and model their interactions to facilitate the prediction. However, these methods fail to fully explore positional information of splice junctions and their deep interaction.ResultsWe present a robust end-to-end framework, Junction Encoder with Deep Interaction (JEDI), for circRNA prediction using only nucleotide sequences. JEDI first leverages the attention mechanism to encode each junction site based on deep bidirectional recurrent neural networks and then presents the novel cross-attention layer to model deep interaction among these sites for backsplicing. Finally, JEDI can not only predict circRNAs but also interpret relationships among splice sites to discover backsplicing hotspots within a gene region. Experiments demonstrate JEDI significantly outperforms state-of-the-art approaches in circRNA prediction on both isoform level and gene level. Moreover, JEDI also shows promising results on zero-shot backsplicing discovery, where none of the existing approaches can achieve.Availability and implementationThe implementation of our framework is available at https://github.com/hallogameboy/JEDI.Supplementary informationSupplementary data are available at Bioinformatics online.

RNA splicing and spliceosome assembly in eukaryotes occur mainly during transcription. However, co-transcriptional splicing has not yet been explored in plants. Here, we built transcriptomes of nascent chromatin RNAs in Arabidopsis thaliana and showed that nearly all introns undergo co-transcriptional splicing, which occurs with higher efficiency for introns in protein-coding genes than for those in noncoding RNAs. Total intron number and intron position are two predominant features that correlate with co-transcriptional splicing efficiency, and introns with alternative 5' or 3' splice sites are less efficiently spliced. Furthermore, we found that mutations in genes encoding trans-acting proteins lead to more introns with increased splicing defects in nascent RNAs than in mature RNAs, and that introns with increased splicing defects in mature RNAs are inefficiently spliced at the co-transcriptional level. Collectively, our results not only uncovered widespread co-transcriptional splicing in Arabidopsis but also identified features that may affect or be affected by co-transcriptional splicing efficiency.

Although immune checkpoint blockade (ICB) therapy has revolutionized cancer treatment, many patients do not respond or develop resistance to ICB. N6 -methylation of adenosine (m6A) in RNA regulates many pathophysiological processes. Here, we show that deletion of the m6A demethylase Alkbh5 sensitized tumors to cancer immunotherapy. Alkbh5 has effects on m6A density and splicing events in tumors during ICB. Alkbh5 modulates Mct4/Slc16a3 expression and lactate content of the tumor microenvironment and the composition of tumor-infiltrating Treg and myeloid-derived suppressor cells. Importantly, a small-molecule Alkbh5 inhibitor enhanced the efficacy of cancer immunotherapy. Notably, the ALKBH5 gene mutation and expression status of melanoma patients correlate with their response to immunotherapy. Our results suggest that m6A demethylases in tumor cells contribute to the efficacy of immunotherapy and identify ALKBH5 as a potential therapeutic target to enhance immunotherapy outcome in melanoma, colorectal, and potentially other cancers.

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

The decoy exon model has been proposed to regulate a subset of intron retention (IR) events involving predominantly larger introns (>1 kb). Splicing reporter studies have shown that decoy splice sites are essential for activity, suggesting that decoys act by engaging intron-terminal splice sites and competing with cross-intron interactions required for intron excision. The decoy model predicts that antisense oligonucleotides may be able to block decoy splice sites in endogenous pre-mRNA, thereby reducing IR and increasing productive gene expression. Indeed, we now demonstrate that targeting a decoy 5' splice site in the O-GlcNAc transferase (OGT) gene reduced IR from ∼80% to ∼20% in primary human erythroblasts, accompanied by increases in spliced OGT RNA and OGT protein expression. The remaining OGT IR was refractory to antisense treatment and might be mediated by independent mechanism(s). In contrast, other retained introns were strongly dependent on decoy function, since antisense targeting of decoy 5' splice sites greatly reduced (SNRNP70) or nearly eliminated (SF3B1) IR in two widely expressed splicing factors, and also greatly reduced IR in transcripts encoding the erythroid-specific structural protein, α-spectrin (SPTA1). These results show that modulating decoy exon function can dramatically alter IR and suggest that dynamic regulation of decoy exons could be a mechanism to fine-tune gene expression post-transcriptionally in many cell types.

Although immune checkpoint blockade (ICB) therapy has revolutionized cancer treatment, many patients do not respond or develop resistance to ICB. N6 -methylation of adenosine (m6A) in RNA regulates many pathophysiological processes. Here, we show that deletion of the m6A demethylase Alkbh5 sensitized tumors to cancer immunotherapy. Alkbh5 has effects on m6A density and splicing events in tumors during ICB. Alkbh5 modulates Mct4/Slc16a3 expression and lactate content of the tumor microenvironment and the composition of tumor-infiltrating Treg and myeloid-derived suppressor cells. Importantly, a small-molecule Alkbh5 inhibitor enhanced the efficacy of cancer immunotherapy. Notably, the ALKBH5 gene mutation and expression status of melanoma patients correlate with their response to immunotherapy. Our results suggest that m6A demethylases in tumor cells contribute to the efficacy of immunotherapy and identify ALKBH5 as a potential therapeutic target to enhance immunotherapy outcome in melanoma, colorectal, and potentially other cancers.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

Pre-mRNA splicing must occur with extremely high fidelity. Spliceosomes assemble onto pre-mRNA guided by specific sequences (5' splice site, 3' splice site, and branchpoint). When splice sites are mutated, as in many hereditary diseases, the spliceosome can aberrantly select nearby pseudo- or "cryptic" splice sites, often resulting in nonfunctional protein. How the spliceosome distinguishes authentic splice sites from cryptic splice sites is poorly understood. We performed a Caenorhabditis elegans genetic screen to find cellular factors that affect the frequency with which the spliceosome uses cryptic splice sites and identified two alleles in core spliceosome component Prp8 that alter cryptic splicing frequency. Subsequent complementary genetic and structural analyses in yeast implicate these alleles in the stability of the spliceosome's catalytic core. However, despite a clear effect on cryptic splicing, high-throughput mRNA sequencing of these prp-8 mutant C. elegans reveals that overall alternative splicing patterns are relatively unchanged. Our data suggest the spliceosome evolved intrinsic mechanisms to reduce the occurrence of cryptic splicing and that these mechanisms are distinct from those that impact alternative splicing.

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5' splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5' splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

The epithelial cell adhesion molecule gene (EPCAM, previously known as TACSTD1 or TROP1) encodes a membrane-bound protein that is localized to the basolateral membrane of epithelial cells and is overexpressed in some tumors. Biallelic mutations in EPCAM cause congenital tufting enteropathy (CTE), which is a rare chronic diarrheal disorder presenting in infancy. Monoallelic deletions of the 3' end of EPCAM that silence the downstream gene, MSH2, cause a form of Lynch syndrome, which is a cancer predisposition syndrome associated with loss of DNA mismatch repair. Here, we report 13 novel EPCAM mutations from 17 CTE patients from two separate centers, review EPCAM mutations associated with CTE and Lynch syndrome, and structurally model pathogenic missense mutations. Statistical analyses indicate that the c.499dupC (previously reported as c.498insC) frameshift mutation was associated with more severe treatment regimens and greater mortality in CTE, whereas the c.556-14A>G and c.491+1G>A splice site mutations were not correlated with treatments or outcomes significantly different than random simulation. These findings suggest that genotype-phenotype correlations may be useful in contributing to management decisions of CTE patients. Depending on the type and nature of EPCAM mutation, one of two unrelated diseases may occur, CTE or Lynch syndrome.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

Many viruses shut off host gene expression to inhibit antiviral responses. Viral proteins and host proteins required for viral replication are typically spared in this process, but the mechanisms of target selectivity during host shutoff remain poorly understood. Using transcriptome-wide and targeted reporter experiments, we demonstrate that the influenza A virus endoribonuclease PA-X usurps RNA splicing to selectively target host RNAs for destruction. Proximity-labeling proteomics reveals that PA-X interacts with cellular RNA processing proteins, some of which are partially required for host shutoff. Thus, PA-X taps into host nuclear pre-mRNA processing mechanisms to destroy nascent mRNAs shortly after their synthesis. This mechanism sets PA-X apart from other viral host shutoff proteins that target actively translating mRNAs in the cytoplasm. Our study reveals a unique mechanism of host shutoff that helps us understand how influenza viruses suppress host gene expression.

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

The epithelial cell adhesion molecule gene (EPCAM, previously known as TACSTD1 or TROP1) encodes a membrane-bound protein that is localized to the basolateral membrane of epithelial cells and is overexpressed in some tumors. Biallelic mutations in EPCAM cause congenital tufting enteropathy (CTE), which is a rare chronic diarrheal disorder presenting in infancy. Monoallelic deletions of the 3' end of EPCAM that silence the downstream gene, MSH2, cause a form of Lynch syndrome, which is a cancer predisposition syndrome associated with loss of DNA mismatch repair. Here, we report 13 novel EPCAM mutations from 17 CTE patients from two separate centers, review EPCAM mutations associated with CTE and Lynch syndrome, and structurally model pathogenic missense mutations. Statistical analyses indicate that the c.499dupC (previously reported as c.498insC) frameshift mutation was associated with more severe treatment regimens and greater mortality in CTE, whereas the c.556-14A>G and c.491+1G>A splice site mutations were not correlated with treatments or outcomes significantly different than random simulation. These findings suggest that genotype-phenotype correlations may be useful in contributing to management decisions of CTE patients. Depending on the type and nature of EPCAM mutation, one of two unrelated diseases may occur, CTE or Lynch syndrome.

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5' splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5' splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

Many viruses shut off host gene expression to inhibit antiviral responses. Viral proteins and host proteins required for viral replication are typically spared in this process, but the mechanisms of target selectivity during host shutoff remain poorly understood. Using transcriptome-wide and targeted reporter experiments, we demonstrate that the influenza A virus endoribonuclease PA-X usurps RNA splicing to selectively target host RNAs for destruction. Proximity-labeling proteomics reveals that PA-X interacts with cellular RNA processing proteins, some of which are partially required for host shutoff. Thus, PA-X taps into host nuclear pre-mRNA processing mechanisms to destroy nascent mRNAs shortly after their synthesis. This mechanism sets PA-X apart from other viral host shutoff proteins that target actively translating mRNAs in the cytoplasm. Our study reveals a unique mechanism of host shutoff that helps us understand how influenza viruses suppress host gene expression.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

In a recent issue of Molecular Cell, Boehm et al. (2018), Blazquez et al. (2018), and Gonatopoulos-Pournatzis et al. (2018) uncover novel mechanisms by which the cell regulates splicing of cryptic splice sites and microexons.

Alternative pre-mRNA splicing plays a major role in expanding the transcript output of human genes. This process is regulated, in part, by the interplay of trans-acting RNA binding proteins (RBPs) with myriad cis-regulatory elements scattered throughout pre-mRNAs. These molecular recognition events are critical for defining the protein-coding sequences (exons) within pre-mRNAs and directing spliceosome assembly on noncoding regions (introns). One of the earliest events in this process is recognition of the 3' splice site (3'ss) by U2 small nuclear RNA auxiliary factor 2 (U2AF2). Splicing regulators, such as the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), influence spliceosome assembly both in vitro and in vivo, but their mechanisms of action remain poorly described on a global scale. HNRNPA1 also promotes proofreading of 3'ss sequences though a direct interaction with the U2AF heterodimer. To determine how HNRNPA1 regulates U2AF-RNA interactions in vivo, we analyzed U2AF2 RNA binding specificity using individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) in control and HNRNPA1 overexpression cells. We observed changes in the distribution of U2AF2 crosslinking sites relative to the 3'ss of alternative cassette exons but not constitutive exons upon HNRNPA1 overexpression. A subset of these events shows a concomitant increase of U2AF2 crosslinking at distal intronic regions, suggesting a shift of U2AF2 to "decoy" binding sites. Of the many noncanonical U2AF2 binding sites, Alu-derived RNA sequences represented one of the most abundant classes of HNRNPA1-dependent decoys. We propose that one way HNRNPA1 regulates exon definition is to modulate the interaction of U2AF2 with decoy or bona fide 3'ss.

The tri-snRNP 27K protein is a component of the human U4/U6-U5 tri-snRNP and contains an N-terminal phosphorylated RS domain. In a forward genetic screen in C. elegans, we previously identified a dominant mutation, M141T, in the highly-conserved C-terminal region of this protein. The mutant allele promotes changes in cryptic 5 splice site choice. To better understand the function of this poorly characterized splicing factor, we performed high-throughput mRNA sequencing analysis on worms containing this dominant mutation. Comparison of alternative splice site usage between the mutant and wild-type strains led to the identification of 26 native genes whose splicing changes in the presence of the snrp-27 mutation. The changes in splicing are specific to alternative 5 splice sites. Analysis of new alleles suggests that snrp-27 is an essential gene for worm viability. We performed a novel directed-mutation experiment in which we used the CRISPR-cas9 system to randomly generate mutations specifically at M141 of SNRP-27. We identified eight amino acid substitutions at this position that are viable, and three that are homozygous lethal. All viable substitutions at M141 led to varying degrees of changes in alternative 5 splicing of native targets. We hypothesize a role for this SR-related factor in maintaining the position of the 5 splice site as U1snRNA trades interactions at the 5 end of the intron with U6snRNA and PRP8 as the catalytic site is assembled.

During terminal erythropoiesis, the splicing machinery in differentiating erythroblasts executes a robust intron retention (IR) program that impacts expression of hundreds of genes. We studied IR mechanisms in the SF3B1 splicing factor gene, which expresses ∼50% of its transcripts in late erythroblasts as a nuclear isoform that retains intron 4. RNA-seq analysis of nonsense-mediated decay (NMD)-inhibited cells revealed previously undescribed splice junctions, rare or not detected in normal cells, that connect constitutive exons 4 and 5 to highly conserved cryptic cassette exons within the intron. Minigene splicing reporter assays showed that these cassettes promote IR. Genome-wide analysis of splice junction reads demonstrated that cryptic noncoding cassettes are much more common in large (>1 kb) retained introns than they are in small retained introns or in nonretained introns. Functional assays showed that heterologous cassettes can promote retention of intron 4 in the SF3B1 splicing reporter. Although many of these cryptic exons were spliced inefficiently, they exhibited substantial binding of U2AF1 and U2AF2 adjacent to their splice acceptor sites. We propose that these exons function as decoys that engage the intron-terminal splice sites, thereby blocking cross-intron interactions required for excision. Developmental regulation of decoy function underlies a major component of the erythroblast IR program.

Prp8 is an essential protein that regulates spliceosome assembly and conformation during pre-mRNA splicing. Recent cryo-EM structures of the spliceosome model Prp8 as a scaffold for the spliceosome's catalytic U snRNA components. Using a new amino acid probing strategy, we identified a dynamic region in human Prp8 that is positioned to stabilize the pre-mRNA in the spliceosome active site through interactions with U5 snRNA. Mutagenesis of the identified Prp8 residues in yeast indicates a role in 5' splice site recognition. Genetic interactions with spliceosome proteins Isy1, which buttresses the intron branch point, and Snu114, a regulatory GTPase that directly contacts Prp8, further corroborate a role for the same Prp8 residues in substrate positioning and activation. Together the data suggest that adjustments in interactions between Prp8 and U5 snRNA help establish proper positioning of the pre-mRNA into the active site to enhance 5' splice site fidelity.

We previously investigated the transcriptome and proteome profiles of the murine ocular lens at six developmental time points including two embryonic (E15 and E18) and four postnatal time points (P0, P3, P6, and P9). Here, we extend our analyses to identify novel transcripts and peptides in developing  mouse lens. We identified a total of 9,707 novel transcripts and 325 novel fusion genes in developing mouse lens. Additionally, we identified 13,281 novel alternative splicing (AS) events in mouse lens including 6,990 exon skipping (ES), 2,447 alternative 3' splice site (A3SS), 1,900 alternative 5' splice site (A5SS), 1,771 mutually exclusive exons (MXE), and 173 intron retention (IR). Finally, we integrated our OMIC (Transcriptome and Proteome) datasets identifying 20 novel peptides in mouse lens. All 20 peptides were validated through matching MS/MS spectra of synthetic peptides. To the best of our knowledge, this is the first report integrating OMIC datasets to identify novel peptides in developing murine lens.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

The group II intron and the spliceosome share a common active site architecture and are thought to be evolutionarily related. Here we report the 3.7 Å crystal structure of a eukaryotic group II intron in the lariat-3' exon form, immediately preceding the second step of splicing, analogous to the spliceosomal P complex. This structure reveals the location of the intact 3' splice site within the catalytic core of the group II intron. The 3'-OH of the 5' exon is positioned in close proximity to the 3' splice site for nucleophilic attack and exon ligation. The active site undergoes conformational rearrangements with the catalytic triplex having different configurations before and after the second step of splicing. We describe a complete model for the second step of group II intron splicing that incorporates a dynamic catalytic triplex being responsible for creating the binding pocket for 3' splice site capture.

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

Prp8 is an essential protein that regulates spliceosome assembly and conformation during pre-mRNA splicing. Recent cryo-EM structures of the spliceosome model Prp8 as a scaffold for the spliceosome's catalytic U snRNA components. Using a new amino acid probing strategy, we identified a dynamic region in human Prp8 that is positioned to stabilize the pre-mRNA in the spliceosome active site through interactions with U5 snRNA. Mutagenesis of the identified Prp8 residues in yeast indicates a role in 5' splice site recognition. Genetic interactions with spliceosome proteins Isy1, which buttresses the intron branch point, and Snu114, a regulatory GTPase that directly contacts Prp8, further corroborate a role for the same Prp8 residues in substrate positioning and activation. Together the data suggest that adjustments in interactions between Prp8 and U5 snRNA help establish proper positioning of the pre-mRNA into the active site to enhance 5' splice site fidelity.

The group II intron and the spliceosome share a common active site architecture and are thought to be evolutionarily related. Here we report the 3.7 Å crystal structure of a eukaryotic group II intron in the lariat-3' exon form, immediately preceding the second step of splicing, analogous to the spliceosomal P complex. This structure reveals the location of the intact 3' splice site within the catalytic core of the group II intron. The 3'-OH of the 5' exon is positioned in close proximity to the 3' splice site for nucleophilic attack and exon ligation. The active site undergoes conformational rearrangements with the catalytic triplex having different configurations before and after the second step of splicing. We describe a complete model for the second step of group II intron splicing that incorporates a dynamic catalytic triplex being responsible for creating the binding pocket for 3' splice site capture.

We previously investigated the transcriptome and proteome profiles of the murine ocular lens at six developmental time points including two embryonic (E15 and E18) and four postnatal time points (P0, P3, P6, and P9). Here, we extend our analyses to identify novel transcripts and peptides in developing mouse lens. We identified a total of 9,707 novel transcripts and 325 novel fusion genes in developing mouse lens. Additionally, we identified 13,281 novel alternative splicing (AS) events in mouse lens including 6,990 exon skipping (ES), 2,447 alternative 3' splice site (A3SS), 1,900 alternative 5' splice site (A5SS), 1,771 mutually exclusive exons (MXE), and 173 intron retention (IR). Finally, we integrated our OMIC (Transcriptome and Proteome) datasets identifying 20 novel peptides in mouse lens. All 20 peptides were validated through matching MS/MS spectra of synthetic peptides. To the best of our knowledge, this is the first report integrating OMIC datasets to identify novel peptides in developing murine lens.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

During terminal erythropoiesis, the splicing machinery in differentiating erythroblasts executes a robust intron retention (IR) program that impacts expression of hundreds of genes. We studied IR mechanisms in the SF3B1 splicing factor gene, which expresses ∼50% of its transcripts in late erythroblasts as a nuclear isoform that retains intron 4. RNA-seq analysis of nonsense-mediated decay (NMD)-inhibited cells revealed previously undescribed splice junctions, rare or not detected in normal cells, that connect constitutive exons 4 and 5 to highly conserved cryptic cassette exons within the intron. Minigene splicing reporter assays showed that these cassettes promote IR. Genome-wide analysis of splice junction reads demonstrated that cryptic noncoding cassettes are much more common in large (>1 kb) retained introns than they are in small retained introns or in nonretained introns. Functional assays showed that heterologous cassettes can promote retention of intron 4 in the SF3B1 splicing reporter. Although many of these cryptic exons were spliced inefficiently, they exhibited substantial binding of U2AF1 and U2AF2 adjacent to their splice acceptor sites. We propose that these exons function as decoys that engage the intron-terminal splice sites, thereby blocking cross-intron interactions required for excision. Developmental regulation of decoy function underlies a major component of the erythroblast IR program.

The spliceosome undergoes dramatic changes in a splicing cycle. Structures of B, Bact, C, C*, and intron lariat spliceosome complexes revealed mechanisms of 5'-splice site (ss) recognition, branching, and intron release, but lacked information on 3'-ss recognition, exon ligation, and exon release. Here we report a cryo-electron microscopy structure of the postcatalytic P complex at 3.3-angstrom resolution, revealing that the 3' ss is mainly recognized through non-Watson-Crick base pairing with the 5' ss and branch point. Furthermore, one or more unidentified proteins become stably associated with the P complex, securing the 3' exon and potentially regulating activity of the helicase Prp22. Prp22 binds nucleotides 15 to 21 in the 3' exon, enabling it to pull the intron-exon or ligated exons in a 3' to 5' direction to achieve 3'-ss proofreading or exon release, respectively.

We report the clinical, biochemical, and molecular findings in two brothers with encephalopathy and multi-systemic disease. Abnormal transferrin glycoforms were suggestive of a type I congenital disorder of glycosylation (CDG). While exome sequencing was negative for CDG related candidate genes, the testing revealed compound heterozygous mutations in the mitochondrial elongation factor G gene (GFM1). One of the mutations had been reported previously while the second, novel variant was found deep in intron 6, activating a cryptic splice site. Functional studies demonstrated decreased GFM1 protein levels, suggested disrupted assembly of mitochondrial complexes III and V and decreased activities of mitochondrial complexes I and IV, all indicating combined OXPHOS deficiency.

U1 snRNP plays a critical role in 5'-splice site recognition and is a frequent target of alternative splicing factors. These factors transiently associate with human U1 snRNP and are not amenable for structural studies, while their Saccharomyces cerevisiae (yeast) homologs are stable components of U1 snRNP. Here, we report the cryoEM structure of yeast U1 snRNP at 3.6 Å resolution with atomic models for ten core proteins, nearly all essential domains of its RNA, and five stably associated auxiliary proteins. The foot-shaped yeast U1 snRNP contains a core in the "ball-and-toes" region architecturally similar to the human U1 snRNP. All auxiliary proteins are in the "arch-and-heel" region and connected to the core through the Prp42/Prp39 paralogs. Our demonstration that homodimeric human PrpF39 directly interacts with U1C-CTD, mirroring yeast Prp42/Prp39, supports yeast U1 snRNP as a model for understanding how transiently associated auxiliary proteins recruit human U1 snRNP in alternative splicing.

A major goal of biomedicine is to understand the function of every gene in the human genome. Loss-of-function mutations can disrupt both copies of a given gene in humans and phenotypic analysis of such 'human knockouts' can provide insight into gene function. Consanguineous unions are more likely to result in offspring carrying homozygous loss-of-function mutations. In Pakistan, consanguinity rates are notably high. Here we sequence the protein-coding regions of 10,503 adult participants in the Pakistan Risk of Myocardial Infarction Study (PROMIS), designed to understand the determinants of cardiometabolic diseases in individuals from South Asia. We identified individuals carrying homozygous predicted loss-of-function (pLoF) mutations, and performed phenotypic analysis involving more than 200 biochemical and disease traits. We enumerated 49,138 rare (<1% minor allele frequency) pLoF mutations. These pLoF mutations are estimated to knock out 1,317 genes, each in at least one participant. Homozygosity for pLoF mutations at PLA2G7 was associated with absent enzymatic activity of soluble lipoprotein-associated phospholipase A2; at CYP2F1, with higher plasma interleukin-8 concentrations; at TREH, with lower concentrations of apoB-containing lipoprotein subfractions; at either A3GALT2 or NRG4, with markedly reduced plasma insulin C-peptide concentrations; and at SLC9A3R1, with mediators of calcium and phosphate signalling. Heterozygous deficiency of APOC3 has been shown to protect against coronary heart disease; we identified APOC3 homozygous pLoF carriers in our cohort. We recruited these human knockouts and challenged them with an oral fat load. Compared with family members lacking the mutation, individuals with APOC3 knocked out displayed marked blunting of the usual post-prandial rise in plasma triglycerides. Overall, these observations provide a roadmap for a 'human knockout project', a systematic effort to understand the phenotypic consequences of complete disrup

U1 snRNP plays a critical role in 5'-splice site recognition and is a frequent target of alternative splicing factors. These factors transiently associate with human U1 snRNP and are not amenable for structural studies, while their Saccharomyces cerevisiae (yeast) homologs are stable components of U1 snRNP. Here, we report the cryoEM structure of yeast U1 snRNP at 3.6 Å resolution with atomic models for ten core proteins, nearly all essential domains of its RNA, and five stably associated auxiliary proteins. The foot-shaped yeast U1 snRNP contains a core in the "ball-and-toes" region architecturally similar to the human U1 snRNP. All auxiliary proteins are in the "arch-and-heel" region and connected to the core through the Prp42/Prp39 paralogs. Our demonstration that homodimeric human PrpF39 directly interacts with U1C-CTD, mirroring yeast Prp42/Prp39, supports yeast U1 snRNP as a model for understanding how transiently associated auxiliary proteins recruit human U1 snRNP in alternative splicing.