LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

Quantifying 3'UTR length from scRNA-seq data reveals changes independent of gene expression

Although more than half of all genes generate transcripts that differ in 3'UTR length, current analysis pipelines only quantify the amount but not the length of mRNA transcripts. 3'UTR length is determined by 3' end cleavage sites (CS). We map CS in more than 200 primary human and mouse cell types and increase CS annotations relative to the GENCODE database by 40%. Approximately half of all CS are used in few cell types, revealing that most genes only have one or two major 3' ends. We incorporate the CS annotations into a computational pipeline, called scUTRquant, for rapid, accurate, and simultaneous quantification of gene and 3'UTR isoform expression from single-cell RNA sequencing (scRNA-seq) data. When applying scUTRquant to data from 474 cell types and 2134 perturbations, we discover extensive 3'UTR length changes across cell types that are as widespread and coordinately regulated as gene expression changes but affect mostly different genes. Our data indicate that mRNA abundance and mRNA length are two largely independent axes of gene regulation that together determine the amount and spatial organization of protein synthesis.

Genome-scale requirements for dynein-based transport revealed by a high-content arrayed CRISPR screen

The microtubule motor dynein plays a key role in cellular organization. However, little is known about how dynein's biosynthesis, assembly, and functional diversity are orchestrated. To address this issue, we have conducted an arrayed CRISPR loss-of-function screen in human cells using the distribution of dynein-tethered peroxisomes and early endosomes as readouts. From a genome-wide gRNA library, 195 validated hits were recovered and parsed into those impacting multiple dynein cargoes and those whose effects are restricted to a subset of cargoes. Clustering of high-dimensional phenotypic fingerprints revealed co-functional proteins involved in many cellular processes, including several candidate novel regulators of core dynein functions. Further analysis of one of these factors, the RNA-binding protein SUGP1, indicates that it promotes cargo trafficking by sustaining functional expression of the dynein activator LIS1. Our data represent a rich source of new hypotheses for investigating microtubule-based transport, as well as several other aspects of cellular organization captured by our high-content imaging.

PABPN1 loss-of-function causes APA-shift in oculopharyngeal muscular dystrophy

Alternative polyadenylation (APA) at the 3' UTR of transcripts contributes to the cell transcriptome. APA is suppressed by the nuclear RNA-binding protein PABPN1. Aging-associated reduced PABPN1 levels in skeletal muscles lead to muscle wasting. Muscle weakness in oculopharyngeal muscular dystrophy (OPMD) is caused by short alanine expansion in PABPN1 exon1. The expanded PABPN1 forms nuclear aggregates, an OPMD hallmark. Whether the expanded PABPN1 affects APA and how it contributes to muscle pathology is unresolved. To investigate these questions, we developed a procedure including RNA library preparation and a simple pipeline calculating the APA-shift ratio as a readout for PABPN1 activity. Comparing APA-shift results to previously published PAS utilization and APA-shift results, we validated this procedure. The procedure was then applied on the OPMD cell model and on RNA from OPMD muscles. APA-shift was genome-wide in the mouse OPMD model, primarily affecting muscle transcripts. In OPMD individuals, APA-shift was enriched with muscle transcripts. In an OPMD cell model APA-shift was not significant. APA-shift correlated with reduced expression levels of a subset of PABPN1 isoforms, whereas the expression of the expanded PABPN1 did not correlate with APA-shift. PABPN1 activity is not affected by the expression of expanded PABPN1, but rather by reduced PABPN1 expression levels. In muscles, PABPN1 activity initially affects muscle transcripts. We suggest that muscle weakness in OPMD is caused by PABPN1 loss-of-function leading to APA-shift that primarily affects in muscle transcripts.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

Genome-wide kinetic profiling of pre-mRNA 3' end cleavage

Cleavage and polyadenylation is necessary for the formation of mature mRNA molecules. The rate at which this process occurs can determine the temporal availability of mRNA for subsequent function throughout the cell and is likely tightly regulated. Despite advances in high-throughput approaches for global kinetic profiling of RNA maturation, genome-wide 3' end cleavage rates have never been measured. Here, we describe a novel approach to estimate the rates of cleavage, using metabolic labeling of nascent RNA, high-throughput sequencing, and mathematical modeling. Using in silico simulations of nascent RNA-seq data, we show that our approach can accurately and precisely estimate cleavage half-lives for both constitutive and alternative sites. We find that 3' end cleavage is fast on average, with half-lives under a minute, but highly variable across individual sites. Rapid cleavage is promoted by the presence of canonical sequence elements and an increased density of polyadenylation signals near a cleavage site. Finally, we find that cleavage rates are associated with the localization of RNA polymerase II at the end of a gene, and faster cleavage leads to quicker degradation of downstream readthrough RNA. Our findings shed light on the features important for efficient 3' end cleavage and the regulation of transcription termination.

Multiplexed screening reveals how cancer-specific alternative polyadenylation shapes tumor growth in vivo

Alternative polyadenylation (APA) is strikingly dysregulated in many cancers. Although global APA dysregulation is frequently associated with poor prognosis, the importance of most individual APA events is controversial simply because few have been functionally studied. Here, we address this gap by developing a CRISPR-Cas9-based screen to manipulate endogenous polyadenylation and systematically quantify how APA events contribute to tumor growth in vivo. Our screen reveals individual APA events that control mouse melanoma growth in an immunocompetent host, with concordant associations in clinical human cancer. For example, forced Atg7 3' UTR lengthening in mouse melanoma suppresses ATG7 protein levels, slows tumor growth, and improves host survival; similarly, in clinical human melanoma, a long ATG7 3' UTR is associated with significantly prolonged patient survival. Overall, our study provides an easily adaptable means to functionally dissect APA in physiological systems and directly quantifies the contributions of recurrent APA events to tumorigenic phenotypes.

TDP-43 loss induces extensive cryptic polyadenylation in ALS/FTD

Nuclear depletion and cytoplasmic aggregation of the RNA-binding protein TDP-43 is the hallmark of ALS, occurring in over 97% of cases. A key consequence of TDP-43 nuclear loss is the de-repression of cryptic exons. Whilst TDP-43 regulated cryptic splicing is increasingly well catalogued, cryptic alternative polyadenylation (APA) events, which define the 3' end of last exons, have been largely overlooked, especially when not associated with novel upstream splice junctions. We developed a novel bioinformatic approach to reliably identify distinct APA event types: alternative last exons (ALE), 3'UTR extensions (3'Ext) and intronic polyadenylation (IPA) events. We identified novel neuronal cryptic APA sites induced by TDP-43 loss of function by systematically applying our pipeline to a compendium of publicly available and in house datasets. We find that TDP-43 binding sites and target motifs are enriched at these cryptic events and that TDP-43 can have both repressive and enhancing action on APA. Importantly, all categories of cryptic APA can also be identified in ALS and FTD post mortem brain regions with TDP-43 proteinopathy underlining their potential disease relevance. RNA-seq and Ribo-seq analyses indicate that distinct cryptic APA categories have different downstream effects on transcript and translation. Intriguingly, cryptic 3'Exts occur in multiple transcription factors, such as ELK1, SIX3, and TLX1, and lead to an increase in wild-type protein levels and function. Finally, we show that an increase in RNA stability leading to a higher cytoplasmic localisation underlies these observations. In summary, we demonstrate that TDP-43 nuclear depletion induces a novel category of cryptic RNA processing events and we expand the palette of TDP-43 loss consequences by showing this can also lead to an increase in normal protein translation.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, limitations, and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for continuous extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies, while the containers and reproducible workflows could easily be deployed and extended to evaluate new methods or data sets.

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

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

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, and limitations and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for seamless extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies. Furthermore, the containers and reproducible workflows generated in the course of this project can be seamlessly deployed and extended in the future to evaluate new methods or datasets.

Genome-scale requirements for dynein-based trafficking revealed by a high-content arrayed CRISPR screen

The cytoplasmic dynein-1 (dynein) motor plays a key role in cellular organisation by transporting a wide variety of cellular constituents towards the minus ends of microtubules. However, relatively little is known about how the biosynthesis, assembly and functional diversity of the motor is orchestrated. To address this issue, we have conducted an arrayed CRISPR loss-of-function screen in human cells using the distribution of dynein-tethered peroxisomes and early endosomes as readouts. From a guide RNA library targeting 18,253 genes, 195 validated hits were recovered and parsed into those impacting multiple dynein cargoes and those whose effects are restricted to a subset of cargoes. Clustering of high-dimensional phenotypic fingerprints generated from multiplexed images revealed co-functional genes involved in many cellular processes, including several candidate novel regulators of core dynein functions. Mechanistic analysis of one of these proteins, the RNA-binding protein SUGP1, provides evidence that it promotes cargo trafficking by sustaining functional expression of the dynein activator LIS1. Our dataset represents a rich source of new hypotheses for investigating microtubule-based transport, as well as several other aspects of cellular organisation that were captured by our high-content imaging.

REPAC: analysis of alternative polyadenylation from RNA-sequencing data

Alternative polyadenylation (APA) is an important post-transcriptional mechanism that has major implications in biological processes and diseases. Although specialized sequencing methods for polyadenylation exist, availability of these data are limited compared to RNA-sequencing data. We developed REPAC, a framework for the analysis of APA from RNA-sequencing data. Using REPAC, we investigate the landscape of APA caused by activation of B cells. We also show that REPAC is faster than alternative methods by at least 7-fold and that it scales well to hundreds of samples. Overall, the REPAC method offers an accurate, easy, and convenient solution for the exploration of APA.

A post-transcriptional regulatory landscape of aging in the female mouse hippocampus

Aging is associated with substantial physiological changes and constitutes a major risk factor for neurological disorders including dementia. Alterations in gene expression upon aging have been extensively studied; however, an in-depth characterization of post-transcriptional regulatory events remains elusive. Here, we profiled the age-related changes of the transcriptome and translatome in the female mouse hippocampus by RNA sequencing of total RNA and polysome preparations at four ages (3-, 6-, 12-, 20-month-old); and we implemented a variety of bioinformatics approaches to unravel alterations in transcript abundance, alternative splicing, and polyadenylation site selection. We observed mostly well-coordinated transcriptome and translatome expression signatures across age including upregulation of transcripts related to immune system processes and neuroinflammation, though transcripts encoding ribonucleoproteins or associated with mitochondrial functions, calcium signaling and the cell-cycle displayed substantial discordant profiles, suggesting translational control associated with age-related deficits in hippocampal-dependent behavior. By contrast, alternative splicing was less preserved, increased with age and was associated with distinct functionally-related transcripts encoding proteins acting at synapses/dendrites, RNA-binding proteins; thereby predicting regulatory roles for RBM3 and CIRBP. Only minor changes in polyadenylation site selection were identified, indicating pivotal 3'-end selection in young adults compared to older groups. Overall, our study provides a comprehensive resource of age-associated post-transcriptional regulatory events in the mouse hippocampus, enabling further examination of the molecular features underlying age-associated neurological diseases.

High-throughput identification of RNA localization elements in neuronal cells

Hundreds of RNAs are enriched in the projections of neuronal cells. For the vast majority of them, though, the sequence elements that regulate their localization are unknown. To identify RNA elements capable of directing transcripts to neurites, we deployed a massively parallel reporter assay that tested the localization regulatory ability of thousands of sequence fragments drawn from endogenous mouse 3' UTRs. We identified peaks of regulatory activity within several 3' UTRs and found that sequences derived from these peaks were both necessary and sufficient for RNA localization to neurites in mouse and human neuronal cells. The localization elements were enriched in adenosine and guanosine residues. They were at least tens to hundreds of nucleotides long as shortening of two identified elements led to significantly reduced activity. Using RNA affinity purification and mass spectrometry, we found that the RNA-binding protein Unk was associated with the localization elements. Depletion of Unk in cells reduced the ability of the elements to drive RNAs to neurites, indicating a functional requirement for Unk in their trafficking. These results provide a framework for the unbiased, high-throughput identification of RNA elements and mechanisms that govern transcript localization in neurons.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

Pseudouridine synthases modify human pre-mRNA co-transcriptionally and affect pre-mRNA processing

Pseudouridine is a modified nucleotide that is prevalent in human mRNAs and is dynamically regulated. Here, we investigate when in their life cycle mRNAs become pseudouridylated to illuminate the potential regulatory functions of endogenous mRNA pseudouridylation. Using single-nucleotide resolution pseudouridine profiling on chromatin-associated RNA from human cells, we identified pseudouridines in nascent pre-mRNA at locations associated with alternatively spliced regions, enriched near splice sites, and overlapping hundreds of binding sites for RNA-binding proteins. In vitro splicing assays establish a direct effect of individual endogenous pre-mRNA pseudouridines on splicing efficiency. We validate hundreds of pre-mRNA sites as direct targets of distinct pseudouridine synthases and show that PUS1, PUS7, and RPUSD4-three pre-mRNA-modifying pseudouridine synthases with tissue-specific expression-control widespread changes in alternative pre-mRNA splicing and 3' end processing. Our results establish a vast potential for cotranscriptional pre-mRNA pseudouridylation to regulate human gene expression via alternative pre-mRNA processing.

Widespread occurrence of hybrid internal-terminal exons in human transcriptomes

Messenger RNA isoform differences are predominantly driven by alternative first, internal, and last exons. Despite the importance of classifying exons to understand isoform structure, few tools examine isoform-specific exon usage. We recently observed that alternative transcription start sites often arise near internal exons, often creating “hybrid” first/internal exons. To systematically detect hybrid exons, we built the hybrid-internal-terminal (HIT) pipeline to classify exons depending on their isoform-specific usage. On the basis of splice junction reads in RNA sequencing data and probabilistic modeling, the HIT index identified thousands of previously misclassified hybrid first-internal and internal-last exons. Hybrid exons are enriched in long genes and genes involved in RNA splicing and have longer flanking introns and strong splice sites. Their usage varies considerably across human tissues. By developing the first method to classify exons according to isoform contexts, our findings document the occurrence of hybrid exons, a common quirk of the human transcriptome.

The Role of Alternative Polyadenylation in the Regulation of Subcellular RNA Localization

Alternative polyadenylation (APA) is a widespread and conserved regulatory mechanism that generates diverse 3' ends on mRNA. APA patterns are often tissue specific and play an important role in cellular processes such as cell proliferation, differentiation, and response to stress. Many APA sites are found in 3' UTRs, generating mRNA isoforms with different 3' UTR contents. These alternate 3' UTR isoforms can change how the transcript is regulated, affecting its stability and translation. Since the subcellular localization of a transcript is often regulated by 3' UTR sequences, this implies that APA can also change transcript location. However, this connection between APA and RNA localization has only recently been explored. In this review, we discuss the role of APA in mRNA localization across distinct subcellular compartments. We also discuss current challenges and future advancements that will aid our understanding of how APA affects RNA localization and molecular mechanisms that drive these processes.

The Quaking RNA-binding proteins as regulators of cell differentiation

The RNA-binding protein Quaking (QKI) has emerged as a potent regulator of cellular differentiation in developmental and pathological processes. The QKI gene is itself alternatively spliced to produce three major isoforms, QKI-5, QKI-6, and QKI-7, that possess very distinct functions. Here, we highlight roles of the different QKI isoforms in neuronal, vascular, muscle, and monocyte cell differentiation, and during epithelial-mesenchymal transition in cancer progression. QKI isoforms control cell differentiation through regulating alternative splicing, mRNA stability and translation, with activities in gene transcription now also becoming evident. These diverse functions of the QKI isoforms contribute to their broad influences on RNA metabolism and cellular differentiation. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Development.

The Clp1 R140H mutation alters tRNA metabolism and mRNA 3' processing in mouse models of pontocerebellar hypoplasia

Homozygous mutation of the RNA kinase CLP1 (cleavage factor polyribonucleotide kinase subunit 1) causes pontocerebellar hypoplasia type 10 (PCH10), a pediatric neurodegenerative disease. CLP1 is associated with the transfer RNA (tRNA) splicing endonuclease complex and the cleavage and polyadenylation machinery, but its function remains unclear. We generated two mouse models of PCH10: one homozygous for the disease-associated Clp1 mutation, R140H, and one heterozygous for this mutation and a null allele. Both models exhibit loss of lower motor neurons and neurons of the deep cerebellar nuclei. To explore whether Clp1 mutation impacts tRNA splicing, we profiled the products of intron-containing tRNA genes. While mature tRNAs were expressed at normal levels in mutant mice, numerous other products of intron-containing tRNA genes were dysregulated, with pre-tRNAs, introns, and certain tRNA fragments up-regulated, and other fragments down-regulated. However, the spatiotemporal patterns of dysregulation do not correlate with pathogenicity for most altered tRNA products. To elucidate the effect of Clp1 mutation on precursor messenger RNA (pre-mRNA) cleavage, we analyzed poly(A) site (PAS) usage and gene expression in Clp1 ^R140H/- spinal cord. PAS usage was shifted from proximal to distal sites in the mutant mouse, particularly in short and closely spaced genes. Many such genes were also expressed at lower levels in the Clp1 ^R140H/- mouse, possibly as a result of impaired transcript maturation. These findings are consistent with the hypothesis that select genes are particularly dependent upon CLP1 for proper pre-mRNA cleavage, suggesting that impaired mRNA 3' processing may contribute to pathogenesis in PCH10.