Coverage-dependent bias creates the appearance of binary splicing in single cells
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.
-
Defining the early steps of cardiovascularlineage segregation by single cell RNA-seqNCBI Gene Expression Omnibus, GSE100471.
-
Pseudo-temporal ordering of individual cells reveals regulators of differentiationNCBI Gene Expression Omnibus, GSE52529.
-
Olfactory stem cell differentiation: horizontal basal cell (HBC) lineageNCBI Gene Expression Omnibus, GSE95601.
Article and author information
Author details
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
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)
Further reading
-
- Cell Biology
- Chromosomes and Gene Expression
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.
-
- Cancer Biology
- Chromosomes and Gene Expression
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.