CCAR1 promotes DNA repair via alternative splicing.

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

The Alazami Syndrome-Associated Protein LARP7 Guides U6 Small Nuclear RNA Modification and Contributes to Splicing Robustness.

The La-related protein 7 (LARP7) forms a complex with the nuclear 7SK RNA to regulate RNA polymerase II transcription. It has been implicated in cancer and the Alazami syndrome, a severe developmental disorder. Here, we report a so far unknown role of this protein in RNA modification. We show that LARP7 physically connects the spliceosomal U6 small nuclear RNA (snRNA) with a distinct subset of box C/D small nucleolar RNAs (snoRNAs) guiding U6 2'-O-methylation. Consistently, these modifications are severely compromised in the absence of LARP7. Although general splicing remains largely unaffected, transcriptome-wide analysis revealed perturbations in alternative splicing in LARP7-depleted cells. Importantly, we identified defects in 2'-O-methylation of the U6 snRNA in Alazami syndrome siblings carrying a LARP7 mutation. Our data identify LARP7 as a bridging factor for snoRNA-guided modification of the U6 snRNA and suggest that alterations in splicing fidelity contribute to the etiology of the Alazami syndrome.

T helper cells exhibit a dynamic and reversible 3'-UTR landscape.

3' untranslated regions (3' UTRs) are critical elements of messenger RNAs, as they contain binding sites for RNA-binding proteins (RBPs) and microRNAs that affect various aspects of the RNA life cycle including transcript stability and cellular localization. In response to T cell receptor activation, T cells undergo massive expansion during the effector phase of the immune response and dynamically modify their 3' UTRs. Whether this serves to directly regulate the abundance of specific mRNAs or is a secondary effect of proliferation remains unclear. To study 3'-UTR dynamics in T helper cells, we investigated division-dependent alternative polyadenylation (APA). In addition, we generated 3' end UTR sequencing data from naive, activated, memory, and regulatory CD4+ T cells. 3'-UTR length changes were estimated using a nonnegative matrix factorization approach and were compared with those inferred from long-read PacBio sequencing. We found that APA events were transient and reverted after effector phase expansion. Using an orthogonal bulk RNA-seq data set, we did not find evidence of APA association with differential gene expression or transcript usage, indicating that APA has only a marginal effect on transcript abundance. 3'-UTR sequence analysis revealed conserved binding sites for T cell-relevant microRNAs and RBPs in the alternative 3' UTRs. These results indicate that poly(A) site usage could play an important role in the control of cell fate decisions and homeostasis.

Roquin Suppresses the PI3K-mTOR Signaling Pathway to Inhibit T Helper Cell Differentiation and Conversion of Treg to Tfr Cells.

Roquin proteins preclude spontaneous T cell activation and aberrant differentiation of T follicular helper (Tfh) or T helper 17 (Th17) cells. Here we showed that deletion of Roquin-encoding alleles specifically in regulatory T (Treg) cells also caused the activation of conventional T cells. Roquin-deficient Treg cells downregulated CD25, acquired a follicular Treg (Tfr) cell phenotype, and suppressed germinal center reactions but could not protect from colitis. Roquin inhibited the PI3K-mTOR signaling pathway by upregulation of Pten through interfering with miR-17∼92 binding to an overlapping cis-element in the Pten 3' UTR, and downregulated the Foxo1-specific E3 ubiquitin ligase Itch. Loss of Roquin enhanced Akt-mTOR signaling and protein synthesis, whereas inhibition of PI3K or mTOR in Roquin-deficient T cells corrected enhanced Tfh and Th17 or reduced iTreg cell differentiation. Thereby, Roquin-mediated control of PI3K-mTOR signaling prevents autoimmunity by restraining activation and differentiation of conventional T cells and specialization of Treg cells.

Impaired age-associated mitochondrial translation is mitigated by exercise and PGC-1α.

Sarcopenia, the age-related loss of skeletal muscle mass and function, can dramatically impinge on quality of life and mortality. While mitochondrial dysfunction and imbalanced proteostasis are recognized as hallmarks of sarcopenia, the regulatory and functional link between these processes is underappreciated and unresolved. We therefore investigated how mitochondrial proteostasis, a crucial process that coordinates the expression of nuclear- and mitochondrial-encoded mitochondrial proteins with supercomplex formation and respiratory activity, is affected in skeletal muscle aging. Intriguingly, a robust mitochondrial translation impairment was observed in sarcopenic muscle, which is regulated by the peroxisome proliferator-activated receptor γ coactivator 1 α (PGC-1α) with the estrogen-related receptor α (ERRα). Exercise, a potent inducer of PGC-1α activity, rectifies age-related reduction in mitochondrial translation, in conjunction with quality control pathways. These results highlight the importance of mitochondrial proteostasis in muscle aging, and elucidate regulatory interactions that underlie the powerful benefits of physical activity in this context.

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.

Alternative cleavage and polyadenylation in health and disease.

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP.

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

Prevention of dsRNA-induced interferon signaling by AGO1x is linked to breast cancer cell proliferation.

Translational readthrough, i.e., elongation of polypeptide chains beyond the stop codon, was initially reported for viral RNA, but later found also on eukaryotic transcripts, resulting in proteome diversification and protein-level modulation. Here, we report that AGO1x, an evolutionarily conserved translational readthrough isoform of Argonaute 1, is generated in highly proliferative breast cancer cells, where it curbs accumulation of double-stranded RNAs (dsRNAs) and consequent induction of interferon responses and apoptosis. In contrast to other mammalian Argonaute protein family members with primarily cytoplasmic functions, AGO1x exhibits nuclear localization in the vicinity of nucleoli. We identify AGO1x interaction with the polyribonucleotide nucleotidyltransferase 1 (PNPT1) and show that the depletion of this protein further augments dsRNA accumulation. Our study thus uncovers a novel function of an Argonaute protein in buffering the endogenous dsRNA-induced interferon responses, different than the canonical function of AGO proteins in the miRNA effector pathway. As AGO1x expression is tightly linked to breast cancer cell proliferation, our study thus suggests a new direction for limiting tumor growth.

CFIm-mediated alternative polyadenylation remodels cellular signaling and miRNA biogenesis.

The mammalian cleavage factor I (CFIm) has been implicated in alternative polyadenylation (APA) in a broad range of contexts, from cancers to learning deficits and parasite infections. To determine how the CFIm expression levels are translated into these diverse phenotypes, we carried out a multi-omics analysis of cell lines in which the CFIm25 (NUDT21) or CFIm68 (CPSF6) subunits were either repressed by siRNA-mediated knockdown or over-expressed from stably integrated constructs. We established that >800 genes undergo coherent APA in response to changes in CFIm levels, and they cluster in distinct functional classes related to protein metabolism. The activity of the ERK pathway traces the CFIm concentration, and explains some of the fluctuations in cell growth and metabolism that are observed upon CFIm perturbations. Furthermore, multiple transcripts encoding proteins from the miRNA pathway are targets of CFIm-dependent APA. This leads to an increased biogenesis and repressive activity of miRNAs at the same time as some 3' UTRs become shorter and presumably less sensitive to miRNA-mediated repression. Our study provides a first systematic assessment of a core set of APA targets that respond coherently to changes in CFIm protein subunit levels (CFIm25/CFIm68). We describe the elicited signaling pathways downstream of CFIm, which improve our understanding of the key role of CFIm in integrating RNA processing with other cellular activities.

Shortening of 3' UTRs in most cell types composing tumor tissues implicates alternative polyadenylation in protein metabolism.

During pre-mRNA maturation 3' end processing can occur at different polyadenylation sites in the 3' untranslated region (3' UTR) to give rise to transcript isoforms that differ in the length of their 3' UTRs. Longer 3' UTRs contain additional cis-regulatory elements that impact the fate of the transcript and/or of the resulting protein. Extensive alternative polyadenylation (APA) has been observed in cancers, but the mechanisms and roles remain elusive. In particular, it is unclear whether the APA occurs in the malignant cells or in other cell types that infiltrate the tumor. To resolve this, we developed a computational method, called SCUREL, that quantifies changes in 3' UTR length between groups of cells, including cells of the same type originating from tumor and control tissue. We used this method to study APA in human lung adenocarcinoma (LUAD). SCUREL relies solely on annotated 3' UTRs and on control systems such as T cell activation, and spermatogenesis gives qualitatively similar results at much greater sensitivity compared to the previously published scAPA method. In the LUAD samples, we find a general trend toward 3' UTR shortening not only in cancer cells compared to the cell type of origin, but also when comparing other cell types from the tumor vs. the control tissue environment. However, we also find high variability in the individual targets between patients. The findings help in understanding the extent and impact of APA in LUAD, which may support improvements in diagnosis and treatment.

Loss of TSC complex enhances gluconeogenesis via upregulation of Dlk1-Dio3 locus miRNAs.

Loss of the tumor suppressor tuberous sclerosis complex 1 (Tsc1) in the liver promotes gluconeogenesis and glucose intolerance. We asked whether this could be attributed to aberrant expression of small RNAs. We performed small-RNA sequencing on liver of Tsc1-knockout mice, and found that miRNAs of the delta-like homolog 1 (Dlk1)-deiodinase iodothyronine type III (Dio3) locus are up-regulated in an mTORC1-dependent manner. Sustained mTORC1 signaling during development prevented CpG methylation and silencing of the Dlk1-Dio3 locus, thereby increasing miRNA transcription. Deletion of miRNAs encoded by the Dlk1-Dio3 locus reduced gluconeogenesis, glucose intolerance, and fasting blood glucose levels. Thus, miRNAs contribute to the metabolic effects observed upon loss of TSC1 and hyperactivation of mTORC1 in the liver. Furthermore, we show that miRNA is a downstream effector of hyperactive mTORC1 signaling.

PolyASite 2.0: a consolidated atlas of polyadenylation sites from 3' end sequencing.

Generated by 3' end cleavage and polyadenylation at alternative polyadenylation (poly(A)) sites, alternative terminal exons account for much of the variation between human transcript isoforms. More than a dozen protocols have been developed so far for capturing and sequencing RNA 3' ends from a variety of cell types and species. In previous studies, we have used these data to uncover novel regulatory signals and cell type-specific isoforms. Here we present an update of the PolyASite (https://polyasite.unibas.ch) resource of poly(A) sites, constructed from publicly available human, mouse and worm 3' end sequencing datasets by enforcing uniform quality measures, including the flagging of putative internal priming sites. Through integrated processing of all data, we identified and clustered sites that are closely spaced and share polyadenylation signals, as these are likely the result of stochastic variations in processing. For each cluster, we identified the representative - most frequently processed - site and estimated the relative use in the transcriptome across all samples. We have established a modern web portal for efficient finding, exploration and export of data. Database generation is fully automated, greatly facilitating incorporation of new datasets and the updating of underlying genome resources.

Protein synthesis rates and ribosome occupancies reveal determinants of translation elongation rates.

Although protein synthesis dynamics has been studied both with theoretical models and by profiling ribosome footprints, the determinants of ribosome flux along open reading frames (ORFs) are not fully understood. Combining measurements of protein synthesis rate with ribosome footprinting data, we here inferred translation initiation and elongation rates for over a 1,000 ORFs in exponentially growing wild-type yeast cells. We found that the amino acid composition of synthesized proteins is as important a determinant of translation elongation rate as parameters related to codon and transfer RNA (tRNA) adaptation. We did not find evidence of ribosome collisions curbing the protein output of yeast transcripts, either in high translation conditions associated with exponential growth, or in strains in which deletion of individual ribosomal protein (RP) genes leads to globally increased or decreased translation. Slow translation elongation is characteristic of RP-encoding transcripts, which have markedly lower protein output compared with other transcripts with equally high ribosome densities.

Reconstitution of CPSF active in polyadenylation: recognition of the polyadenylation signal by WDR33.

Cleavage and polyadenylation specificity factor (CPSF) is the central component of the 3' processing machinery for polyadenylated mRNAs in metazoans: CPSF recognizes the polyadenylation signal AAUAAA, providing sequence specificity in both pre-mRNA cleavage and polyadenylation, and catalyzes pre-mRNA cleavage. Here we show that of the seven polypeptides that have been proposed to constitute CPSF, only four (CPSF160, CPSF30, hFip1, and WDR33) are necessary and sufficient to reconstitute a CPSF subcomplex active in AAUAAA-dependent polyadenylation, whereas CPSF100, CPSF73, and symplekin are dispensable. WDR33 is required for binding of reconstituted CPSF to AAUAAA-containing RNA and can be specifically UV cross-linked to such RNAs, as can CPSF30. Transcriptome-wide identification of WDR33 targets by photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) showed that WDR33 binds in and very close to the AAUAAA signal in vivo with high specificity. Thus, our data indicate that the large CPSF subunit participating in recognition of the polyadenylation signal is WDR33 and not CPSF160, as suggested by previous studies.

The transcriptional landscape of a hepatoma cell line grown on scaffolds of extracellular matrix proteins.

BACKGROUND: The behavior of cells in vivo is complex and highly dynamic, as it results from an interplay between intercellular matrix proteins with surface receptors and other microenvironmental cues. Although the effects of the cellular niche have been investigated for a number of cell types using different molecular approaches, comprehensive assessments of how the global transcriptome responds to 3D scaffolds composed of various extracellular matrix (ECM) constituents at different concentrations are still lacking. RESULTS: In this study, we explored the effects of two diverse extracellular matrix (ECM) components, Collagen I and Matrigel, on the transcriptional profile of cells in a cell culture system. Culturing Huh-7 cells on traditional cell culture plates (Control) or on the ECM components at different concentrations to modulate microenvironment properties, we have generated transcriptomics data that may be further explored to understand the differentiation and growth potential of this cell type for the development of 3D cultures. Our analysis infers transcription factors that are most responsible for the transcriptome response to the extracellular cues. CONCLUSION: Our data indicates that the Collagen I substrate induces a robust transcriptional response in the Huh-7 cells, distinct from that induced by Matrigel. Enhanced hepatocyte markers (ALB and miR-122) reveal a potentially robust remodelling towards primary hepatocytes. Our results aid in defining the appropriate culture and transcription pathways while using hepatoma cell lines. As systems mimicking the in vivo structure and function of liver cells are still being developed, our study could potentially circumvent bottlenecks of limited availability of primary hepatocytes for preclinical studies of drug targets.

Terminal exon characterization with TECtool reveals an abundance of cell-specific isoforms.

Sequencing of RNA 3' ends has uncovered numerous sites that do not correspond to the termination sites of known transcripts. Through their 3' untranslated regions, protein-coding RNAs interact with RNA-binding proteins and microRNAs, which regulate many properties, including RNA stability and subcellular localization. We developed the terminal exon characterization (TEC) tool ( http://tectool.unibas.ch ), which can be used with RNA-sequencing data from any species for which a genome annotation that includes sites of RNA cleavage and polyadenylation is available. We discovered hundreds of previously unknown isoforms and cell-type-specific terminal exons in human cells. Ribosome profiling data revealed that many of these isoforms were translated. By applying TECtool to single-cell sequencing data, we found that the newly identified isoforms were expressed in subpopulations of cells. Thus, TECtool enables the identification of previously unknown isoforms in well-studied cell systems and in rare cell types.

TUT-DIS3L2 is a mammalian surveillance pathway for aberrant structured non-coding RNAs.

Uridylation of various cellular RNA species at the 3' end has been generally linked to RNA degradation. In mammals, uridylated pre-let-7 miRNAs and mRNAs are targeted by the 3' to 5' exoribonuclease DIS3L2. Mutations in DIS3L2 have been associated with Perlman syndrome and with Wilms tumor susceptibility. Using in vivo cross-linking and immunoprecipitation (CLIP) method, we discovered the DIS3L2-dependent cytoplasmic uridylome of human cells. We found a broad spectrum of uridylated RNAs including rRNAs, snRNAs, snoRNAs, tRNAs, vault, 7SL, Y RNAs, mRNAs, lncRNAs, and transcripts from pseudogenes. The unifying features of most of these identified RNAs are aberrant processing and the presence of stable secondary structures. Most importantly, we demonstrate that uridylation mediates DIS3L2 degradation of short RNA polymerase II-derived RNAs. Our findings establish the role of DIS3L2 and oligouridylation as the cytoplasmic quality control for highly structured ncRNAs.