Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome

Abstract

Alternative splicing (AS) is a fundamental step in eukaryotic mRNA biogenesis. Here, we develop an efficient and reproducible pipeline for the discovery of genetic variants that affect AS (splicing QTLs, sQTLs). We use it to analyze the GTEx dataset, generating a comprehensive catalog of sQTLs in the human genome. Downstream analysis of this catalog provides insight into the mechanisms underlying splicing regulation. We report that a core set of sQTLs is shared across multiple tissues. sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many also affect the expression of the same or other genes, uncovering regulatory loci that act through different mechanisms. sQTLs tend to be located in post-transcriptionally spliced introns, which would function as hotspots for splicing regulation. While many variants affect splicing patterns by altering the sequence of splice sites, many more modify the binding sites of RNA-binding proteins. Genetic variants affecting splicing can have a stronger phenotypic impact than those affecting gene expression.

References 87

1. 1.

   Nilsen TW, Graveley BR. Expansion of the eukaryotic proteome by alternative splicing. Nature. 2010;463:457-463.

   Not in our library10.1038/nature08909

2. 2.

   Keren H, Lev-Maor G, Ast G. Alternative splicing and evolution: diversification, exon definition and function. Nat. Rev. Genet. 2010;11:345-55.

   Not in our library10.1038/nrg2776

3. 3.

   Scotti MM, Swanson MS. RNA mis-splicing in disease. Nat. Rev. Genet. 2016;17:19-32.

   Not in our library10.1038/nrg.2015.3

4. 4.

   Chen M, Manley JL. Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat. Rev. Mol. Cell Biol. 2009;10:741-754.

   Not in our library10.1038/nrm2777

5. 5.

   Fu X-D, Ares M. Context-dependent control of alternative splicing by RNA-binding proteins. Nat. Rev. Genet. 2014;15:689-701.

   Not in our library10.1038/nrg3778

6. 6.

   Kornblihtt AR, et al. Alternative splicing: a pivotal step between eukaryotic transcription and translation. Nat. Rev. Mol. Cell Biol. 2013;14:153-165.

   Not in our library10.1038/nrm3525

7. 7.

   Lappalainen T, et al. Transcriptome and genome sequencing uncovers functional variation in humans. Nature. 2013;501:506-11.

   Not in our library10.1038/nature12531

8. 8.

   Battle A, et al. Characterizing the genetic basis of transcriptome diversity through RNA-sequencing of 922 individuals. Genome Res. 2014;24:14-24.

   Not in our library10.1101/gr.155192.113

9. 9.
   Genome-wide identification of splicing QTLs in the human brain and their enrichment among schizophrenia-associated loci · 2017

   No notice
10. 10.

    Raj T, et al. Integrative transcriptome analyses of the aging brain implicate altered splicing in Alzheimer’s disease susceptibility. Nat. Genet. 2018;50:1584-1592.

    Not in our library10.1038/s41588-018-0238-1

11. 11.

    Tian J, et al. CancerSplicingQTL: a database for genome-wide identification of splicing QTLs in human cancer. Nucleic Acids Res. 2019;47:D909-D916.

    Not in our library10.1093/nar/gky954

12. 12.

    The GTEx Consortium. The GTEx Consortium atlas of genetic regulatory effects across human tissues. Science369, 1318–1330 (2020).

    Not in our library10.1126/science.aaz1776

13. 13.

    Ma L, Jia P, Zhao Z. Splicing QTL of human adipose-related traits. Sci. Rep. 2018;8:318.

    Not in our library10.1038/s41598-017-18767-z

14. 14.

    Caswell JL, et al. Multiple breast cancer risk variants are associated with differential transcript isoform expression in tumors. Human Mol. Genet. 2015;24:7421-31.

    Not in our library10.1093/hmg/ddv432

15. 15.

    Li YI, et al. RNA splicing is a primary link between genetic variation and disease. Science. 2016;352:600-604.

    Not in our library10.1126/science.aad9417

16. 16.

    Rotival M, Quach H, Quintana-Murci L. Defining the genetic and evolutionary architecture of alternative splicing in response to infection. Nat. Commun. 2019;10:1671.

    Not in our library10.1038/s41467-019-09689-7

17. 17.

    Ongen H, Dermitzakis ET. Alternative splicing QTLs in European and African populations. Am. J. Human Genet. 2015;97:567-75.

    Not in our library10.1016/j.ajhg.2015.09.004

18. 18.

    Li YI, et al. Annotation-free quantification of RNA splicing using LeafCutter. Nat. Genet. 2018;50:151-158.

    Not in our library10.1038/s41588-017-0004-9

19. 19.

    López-Bigas N, Audit B, Ouzounis C, Parra G, Guigó R. Are splicing mutations the most frequent cause of hereditary disease? FEBS Lett. 2005;579:1900-1903.

    Not in our library10.1016/j.febslet.2005.02.047

20. 20.

    Monlong J, Calvo M, Ferreira PG, Guigó R. Identification of genetic variants associated with alternative splicing using sQTLseekeR. Nat. Commun. 2014;5:4698.

    Not in our library10.1038/ncomms5698

21. 21.

    Anderson M. A new method for non-parametric multivariate analysis of variance. Australian Ecol. 2001;26:32-46.

    Not in our library

22. 22.

    Anderson MJ, Robinson J. Generalized discriminant analysis based on distances. Australian N. Zealand J. Stat. 2003;45:301-318.

    Not in our library10.1111/1467-842x.00285

23. 23.

    Di Tommaso P, et al. Nextflow enables reproducible computational workflows. Nat. Biotechnol. 2017;35:316-319.

    Not in our library10.1038/nbt.3820

24. 24.
    Genetic effects on gene expression across human tissues · 2017

    Corrected
25. 25.
    xCell: digitally portraying the tissue cellular heterogeneity landscape · 2017

    No notice
26. 26.

    Jangi M, Sharp P. Building robust transcriptomes with master splicing factors. Cell. 2014;159:487-498.

    Not in our library10.1016/j.cell.2014.09.054

27. 27.
    Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells · 2016

    No notice
28. 28.
    Alternative start and termination sites of transcription drive most transcript isoform differences across human tissues · 2017

    No notice
29. 29.

    The GTEx Consortium. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science348, 648–60 (2015).

    Not in our library10.1126/science.1262110

30. 30.

    Urbut SM, Wang G, Carbonetto P, Stephens M. Flexible statistical methods for estimating and testing effects in genomic studies with multiple conditions. Nat. Genet. 2019;51:187-195.

    Not in our library10.1038/s41588-018-0268-8

31. 31.

    Yanai I, et al. Genome-wide midrange transcription profiles reveal expression level relationships in human tissue specification. Bioinformatics. 2005;21:650-659.

    Not in our library10.1093/bioinformatics/bti042

32. 32.

    Adzhubei IA, et al. A method and server for predicting damaging missense mutations. Nat. Methods. 2010;7:248-249.

    Not in our library10.1038/nmeth0410-248

33. 33.

    Sim N-L, et al. SIFT web server: predicting effects of amino acid substitutions on proteins. Nucleic Acids Res. 2012;40:W452-W457.

    Not in our library10.1093/nar/gks539

34. 34.

    Van Nostrand EL, et al. Robust transcriptome-wide discovery of RNA-binding protein binding sites with enhanced CLIP (eCLIP). Nat. Methods. 2016;13:508-514.

    Not in our library10.1038/nmeth.3810

35. 35.
    Enhanced Regulatory Sequence Prediction Using Gapped k-mer Features · 2014

    Corrected
36. 36.

    Lee D, et al. A method to predict the impact of regulatory variants from DNA sequence. Nat. Genet. 2015;47:955-961.

    Not in our library10.1038/ng.3331

37. 37.

    Ghandi M, et al. gkmSVM: an R package for gapped-kmer SVM. Bioinformatics. 2016;32:2205-7.

    Not in our library10.1093/bioinformatics/btw203

38. 38.

    Yang E-W, et al. Allele-specific binding of RNA-binding proteins reveals functional genetic variants in the RNA. Nat. Commun. 2019;10:1338.

    Not in our library10.1038/s41467-019-09292-w

39. 39.

    Wickramasinghe VO, et al. Regulation of constitutive and alternative mRNA splicing across the human transcriptome by PRPF8 is determined by 5’ splice site strength. Genome Biol. 2015;16:201.

    Not in our library10.1186/s13059-015-0749-3

40. 40.

    Tilgner H, et al. Deep sequencing of subcellular RNA fractions shows splicing to be predominantly co-transcriptional in the human genome but inefficient for lncRNAs. Genome Res. 2012;22:1616-1625.

    Not in our library10.1101/gr.134445.111

Cited by papers we hold 3

• Genetic control of RNA splicing and its distinct role in complex trait variation · 2022
  No notice
• Extensive Variation in Gene Expression is Revealed in 13 Fertility-Related Genes Using RNA-Seq, ISO-Seq, and CAGE-Seq From Brahman Cattle · 2022
  Not checked
• Clinical implementation of RNA sequencing for Mendelian disease diagnostics · 2022
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.