RESUMEN
Alternative polyadenylation (APA) is emerging as a widespread regulatory layer since the majority of human protein-coding genes contain several polyadenylation (p(A)) sites in their 3'UTRs. By generating isoforms with different 3'UTR length, APA potentially affects mRNA stability, translation efficiency, nuclear export, and cellular localization. Polyadenylation sites are regulated by adjacent RNA cis-regulatory elements, the principals among them are the polyadenylation signal (PAS) AAUAAA and its main variant AUUAAA, typically located ~20-nt upstream of the p(A) site. Mutations in PAS and other auxiliary poly(A) cis-elements in the 3'UTR of several genes have been shown to cause human Mendelian diseases, and to date, only a few common SNPs that regulate APA were associated with complex diseases. Here, we systematically searched for SNPs that affect gene expression and human traits by modulation of 3'UTR APA. First, focusing on the variants most likely to exert the strongest effect, we identified 2,305 SNPs that interrupt the canonical PAS or its main variant. Implementing pA-QTL tests using GTEx RNA-seq data, we identified 330 PAS SNPs (called PAS pA-QTLs) that were significantly associated with the usage of their p(A) site. As expected, PAS-interrupting alleles were mostly linked with decreased cleavage at their p(A) site and the consequential 3'UTR lengthening. However, interestingly, in ~10% of the cases, the PAS-interrupting allele was associated with increased usage of an upstream p(A) site and 3'UTR shortening. As an indication of the functional effects of these PAS pA-QTLs on gene expression and complex human traits, we observed for few dozens of them marked colocalization with eQTL and/or GWAS signals. The PAS-interrupting alleles linked with 3'UTR lengthening were also strongly associated with decreased gene expression, indicating that shorter isoforms generated by APA are generally more stable than longer ones. Last, we carried out an extended, genome-wide analysis of 3'UTR variants and detected thousands of additional pA-QTLs having weaker effects compared to the PAS pA-QTLs.
Asunto(s)
Polimorfismo de Nucleótido Simple/genética , Señales de Poliadenilación de ARN 3'/genética , Estabilidad del ARN/genética , Regulación de la Expresión Génica/genética , Humanos , Poli A , Poliadenilación/genética , ARN Mensajero/genéticaRESUMEN
Alternative polyadenylation (APA) is emerging as an important layer of gene regulation because the majority of mammalian protein-coding genes contain multiple polyadenylation (pA) sites in their 3' UTR. By alteration of 3' UTR length, APA can considerably affect post-transcriptional gene regulation. Yet, our understanding of APA remains rudimentary. Novel single-cell RNA sequencing (scRNA-seq) techniques allow molecular characterization of different cell types to an unprecedented degree. Notably, the most popular scRNA-seq protocols specifically sequence the 3' end of transcripts. Building on this property, we implemented a method for analysing patterns of APA regulation from such data. Analyzing multiple datasets from diverse tissues, we identified widespread modulation of APA in different cell types resulting in global 3' UTR shortening/lengthening and enhanced cleavage at intronic pA sites. Our results provide a proof-of-concept demonstration that the huge volume of scRNA-seq data that accumulates in the public domain offers a unique resource for the exploration of APA based on a very broad collection of cell types and biological conditions.
Asunto(s)
Poliadenilación/genética , ARN Mensajero/genética , Análisis de Secuencia de ARN/métodos , Análisis de la Célula Individual/métodos , Transcriptoma/genética , Regiones no Traducidas 3'/genética , Animales , Linaje de la Célula/genética , Regulación de la Expresión Génica/genética , HumanosRESUMEN
Genome-wide association studies (GWAS) have identified numerous common genetic variants associated with complex human traits and diseases. However, the translation of GWAS discoveries into biological and clinical insights is highly challenging. In this study, we present a novel bioinformatics approach for enhancing the functional interpretation of GWAS signals, based on their integration with single-cell (sc)RNA-seq datasets that examine developmental processes. Our approach performs three tasks: (1) Identification of links between cell differentiation trajectories and traits; (2) Elucidation of biological processes and molecular pathways that underlie such trajectory-trait links; and (3) Prioritization of target genes that carry the links between trajectories, pathways and traits. We applied our method to a set of 11 traits of various pathologies, and 12 scRNA-seq datasets of diverse developmental processes, and it readily detected well-established biological connections, including those between the maturation of cortical inhibitory interneurons and schizophrenia, hepatocytes and cholesterol levels, and pancreatic beta-islet cells and type-2 diabetes. For each of these associations, our method pinpointed top candidate genes that are strongly associated with both the kinetics of the differentiation trajectory and the disease's genetic risk. By the identification of trajectory-disease links, molecular pathways that underlie them and prioritizing candidate risk genes, our method improves the understanding of the etiology of complex diseases, and thus holds promise for enhancing rational drug development that is aimed at targeting specific biological processes that mediate the genetic predisposition to diseases.