Published June 8, 2021
| Published + Supplemental Material + Submitted
Journal Article
Open
BUTTERFLY: addressing the pooled amplification paradox with unique molecular identifiers in single-cell RNA-seq
Chicago
An error occurred while generating the citation.
Abstract
The incorporation of unique molecular identifiers (UMIs) in single-cell RNA-seq assays makes possible the identification of duplicated molecules, thereby facilitating the counting of distinct molecules from sequenced reads. However, we show that the naïve removal of duplicates can lead to a bias due to a "pooled amplification paradox," and we propose an improved quantification method based on unseen species modeling. Our correction called BUTTERFLY uses a zero truncated negative binomial estimator implemented in the kallisto bustools workflow. We demonstrate its efficacy across cell types and genes and show that in some cases it can invert the relative abundance of genes.
Additional Information
© The Author(s). 2021. This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data. Received 27 August 2020; Accepted 21 May 2021; Published 08 June 2021. We thank Pall Melsted, Sina Booeshaghi, and Joseph Min for helpful suggestions on the project and on the integration of BUTTERFLY in bustools. Availability of data and materials: Means to access the datasets analyzed during the current study are listed in Additional file 1: Table S2. The source code as well as Jupyter notebooks for generating the figures is available in GitHub [36], as well as the source code for the branch of bustools used in this project [37]. Snapshots of the repositories are available in Zenodo [38]. Jupiter notebooks have not been produced for Fig. 5 and the data generation for Additional file 1: Fig. S6, S24, and S25, due to difficulties in setting up the right versions of R packages, but the code is directly available in the Github repository. All code is released under the BSD 2-clause license, except for the R file "modZTNB.R," which is released under a GPL3 license to comply with the GPL3 license of PreseqR. Review history: The review history is available as Additional file 3. Peer review information: Barbara Cheifet was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team. This work was supported by funding from the Knut and Alice Wallenberg foundation (J.N.), the National Cancer Institute of the National Institutes of Health under award number F32CA220848 (J.R.), and NIH U19MH114830 (L.P.). Author Contributions: Conceptualization, J.G., L.P.; Methodology, J.G, L.P.; software, J.G. ; writing—original draft, J.G., L.P.; writing—review and editing, J.G., L.P., J.R., J.N.; supervision, L.P., J.R., J.N.; funding acquisition, L.P., J.R., J.N. The authors read and approved the final manuscript. Ethics approval and consent to participate: Not applicable. Consent for publication: Not applicable. The authors declare that they have no competing interests.Attached Files
Published - s13059-021-02386-z.pdf
Submitted - 2020.07.06.188003v2.full.pdf
Supplemental Material - 13059_2021_2386_MOESM1_ESM.pdf
Supplemental Material - 13059_2021_2386_MOESM2_ESM.xlsx
Supplemental Material - 13059_2021_2386_MOESM3_ESM.docx
Files
13059_2021_2386_MOESM1_ESM.pdf
Files
(13.1 MB)
| Name | Size | Download all |
|---|---|---|
|
md5:5cbf238064c1e0e537ed242323aed98d
|
7.9 MB | Preview Download |
|
md5:b4a9f4d0ccf45b163ce584d8931ecc7a
|
1.7 MB | Download |
|
md5:495f28f9e2f8e8f32d3dc337fdf1ccdf
|
20.2 kB | Download |
|
md5:90ac1e6162ee3a93c27dba5c9837ff4f
|
971.0 kB | Preview Download |
|
md5:4e72842c055945bb9c37c59244809c15
|
2.6 MB | Preview Download |
Additional details
- Alternative title
- Addressing the pooled amplification paradox with unique molecular identifiers in single-cell RNA-seq
- PMCID
- PMC8188791
- Eprint ID
- 104254
- Resolver ID
- CaltechAUTHORS:20200707-114817234
- Knut and Alice Wallenberg Foundation
- NIH Postdoctoral Fellowship
- F32CA220848
- NIH
- U19MH114830
- Created
-
2020-07-07Created from EPrint's datestamp field
- Updated
-
2021-07-12Created from EPrint's last_modified field
- Caltech groups
- Division of Biology and Biological Engineering (BBE)