Coverage-dependent bias creates the appearance of binary splicing in single cells

  1. Carlos F Buen Abad Najar
  2. Nir Yosef  Is a corresponding author
  3. Liana F Lareau  Is a corresponding author
  1. University of California, Berkeley, United States

Abstract

Single cell RNA sequencing provides powerful insight into the factors that determine each cell's unique identity. Previous studies led to the surprising observation that alternative splicing among single cells is highly variable and follows a bimodal pattern: a given cell consistently produces either one or the other isoform for a particular splicing choice, with few cells producing both isoforms. Here we show that this pattern arises almost entirely from technical limitations. We analyze alternative splicing in human and mouse single cell RNA-seq datasets, and model them with a probabilistic simulator. Our simulations show that low gene expression and low capture efficiency distort the observed distribution of isoforms. This gives the appearance of binary splicing outcomes, even when the underlying reality is consistent with more than one isoform per cell. We show that accounting for the true amount of information recovered can produce biologically meaningful measurements of splicing in single cells.

Data availability

All sequencing data reanalyzed in this study were acquired from GEO.

The following previously published data sets were used

Article and author information

Author details

  1. Carlos F Buen Abad Najar

    Center for Computational Biology, University of California, Berkeley, Berkeley, United States
    Competing interests
    The authors declare that no competing interests exist.
  2. Nir Yosef

    Department of Electrical Engineering and Computer Science and the Center for Computational Biology, University of California, Berkeley, Berkeley, United States
    For correspondence
    niryosef@berkeley.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9004-1225
  3. Liana F Lareau

    Department of Bioengineering, University of California, Berkeley, Berkeley, United States
    For correspondence
    lareau@berkeley.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3223-3426

Funding

UC MEXUS-Conacyt (Doctoral Fellowship)

  • Carlos F Buen Abad Najar

Chan Zuckerberg Biohub

  • Nir Yosef

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2020, Buen Abad Najar et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,719
    views
  • 583
    downloads
  • 38
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Carlos F Buen Abad Najar
  2. Nir Yosef
  3. Liana F Lareau
(2020)
Coverage-dependent bias creates the appearance of binary splicing in single cells
eLife 9:e54603.
https://doi.org/10.7554/eLife.54603

Share this article

https://doi.org/10.7554/eLife.54603

Further reading

    1. Cell Biology
    2. Chromosomes and Gene Expression
    Bhumil Patel, Maryke Grobler ... Needhi Bhalla
    Research Article

    Meiotic crossover recombination is essential for both accurate chromosome segregation and the generation of new haplotypes for natural selection to act upon. This requirement is known as crossover assurance and is one example of crossover control. While the conserved role of the ATPase, PCH-2, during meiotic prophase has been enigmatic, a universal phenotype when pch-2 or its orthologs are mutated is a change in the number and distribution of meiotic crossovers. Here, we show that PCH-2 controls the number and distribution of crossovers by antagonizing their formation. This antagonism produces different effects at different stages of meiotic prophase: early in meiotic prophase, PCH-2 prevents double-strand breaks from becoming crossover-eligible intermediates, limiting crossover formation at sites of initial double-strand break formation and homolog interactions. Later in meiotic prophase, PCH-2 winnows the number of crossover-eligible intermediates, contributing to the designation of crossovers and ultimately, crossover assurance. We also demonstrate that PCH-2 accomplishes this regulation through the meiotic HORMAD, HIM-3. Our data strongly support a model in which PCH-2’s conserved role is to remodel meiotic HORMADs throughout meiotic prophase to destabilize crossover-eligible precursors and coordinate meiotic recombination with synapsis, ensuring the progressive implementation of meiotic recombination and explaining its function in the pachytene checkpoint and crossover control.

    1. Cancer Biology
    2. Chromosomes and Gene Expression
    Ashley L Cook, Surojit Sur ... Nicolas Wyhs
    Research Article

    Despite exciting developments in cancer immunotherapy, its broad application is limited by the paucity of targetable antigens on the tumor cell surface. As an intrinsic cellular pathway, nonsense-mediated decay (NMD) conceals neoantigens through the destruction of the RNA products from genes harboring truncating mutations. We developed and conducted a high-throughput screen, based on the ratiometric analysis of transcripts, to identify critical mediators of NMD in human cells. This screen implicated disruption of kinase SMG1’s phosphorylation of UPF1 as a potential disruptor of NMD. This led us to design a novel SMG1 inhibitor, KVS0001, that elevates the expression of transcripts and proteins resulting from human and murine truncating mutations in vitro and murine cells in vivo. Most importantly, KVS0001 concomitantly increased the presentation of immune-targetable human leukocyte antigens (HLA) class I-associated peptides from NMD-downregulated proteins on the surface of human cancer cells. KVS0001 provides new opportunities for studying NMD and the diseases in which NMD plays a role, including cancer and inherited diseases.