Mapping and modeling the genomic basis of differential RNA isoform expression at single-cell resolution with LR-Split-seq
Open accessPublished version of a preprint
Abstract
The rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
Emerging long-read sequencing approaches enable profiling cellular transcriptomes at isoform resolution.
“Emerging long-read sequencing approaches enable us to profile cellular transcriptomes at isoform resolution,8,31,32,33 thus providing new opportunities to model the relationships between isoforms using co-expression network analysis.”
From hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data (2023) · cites this paper as [31]
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
The self-check needs the paper’s full text.
Cites
UMAP: Uniform Manifold Approximation and Projection
McInnes, L., Healy, J., Saul, N., & Großberger, L. (2018) · Journal of Open Source Software · doi:10.21105/joss.00861 · Cited by 9,826 · cited 1× · see citing works
“After loading the loom file back into Seurat with the ReadVelocity function from the SeuratWrappers package, SCTransform (v. 0.3.1) was used to regress percent mitochondrial reads, number of genes, and sublibrary, followed by UMAP dimensionality reduction [56, 57].”
MethodInfluentialEnrichr: interactive and collaborative HTML5 gene list enrichment analysis tool
Chen, E. Y., Tan, C. M., Kou, Y., Duan, Q., Wang, Z., Meirelles, G. V., Clark, N. R., & Ma’ayan, A. (2013) · BMC Bioinformatics · doi:10.1186/1471-2105-14-128 · Cited by 7,814 · see citing works
Hafemeister, C., & Satija, R. (2019) · Genome Biology · doi:10.1186/s13059-019-1874-1 · Cited by 4,479 · cited 1× · see citing works
“After loading the loom file back into Seurat with the ReadVelocity function from the SeuratWrappers package, SCTransform (v. 0.3.1) was used to regress percent mitochondrial reads, number of genes, and sublibrary, followed by UMAP dimensionality reduction [56, 57].”
MethodInfluentialOpportunities and challenges in long-read sequencing data analysis
Amarasinghe, S. L., Su, S., Dong, X., Zappia, L., Ritchie, M. E., & Gouil, Q. (2020) · Genome Biology · doi:10.1186/s13059-020-1935-5 · Cited by 2,530 · see citing works
Failure of MBNL1-dependent post-natal splicing transitions in myotonic dystrophy
Lin, X., Miller, J. W., Mankodi, A., Kanadia, R. N., Yuan, Y., Moxley, R. T., Swanson, M. S., & Thornton, C. A. (2006) · Human Molecular Genetics · doi:10.1093/hmg/ddl132 · Cited by 510 · see citing works
Fibronectin Regulates Wnt7a Signaling and Satellite Cell Expansion
Bentzinger, C. F., Wang, Y. X., von Maltzahn, J., Soleimani, V. D., Yin, H., & Rudnicki, M. A. (2013) · Cell Stem Cell · doi:10.1016/j.stem.2012.09.015 · Cited by 386 · see citing works
Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells
Gupta, I., Collier, P. G., Haase, B., Mahfouz, A., Joglekar, A., Floyd, T., Koopmans, F., Barres, B., Smit, A. B., Sloan, S. A., et al. (2018) · Nature Biotechnology · doi:10.1038/nbt.4259 · Cited by 343 · see citing works
A spatially resolved brain region- and cell type-specific isoform atlas of the postnatal mouse brain
Joglekar, A., Prjibelski, A., Mahfouz, A., Collier, P., Lin, S., Schlusche, A. K., Marrocco, J., Williams, S. R., Haase, B., Hayes, A., et al. (2021) · Nature Communications · doi:10.1038/s41467-020-20343-5 · Cited by 229 · see citing works
Flagged first, then by citations
Research works printed but not linked: 58
- [1]Volden R, Palmer T, Byrne A, Cole C, Schmitz RJ, Green RE, et al. Improving nanopore read accuracy with the R2C2 method enables the sequencing of highly multiplexed full-length single-cell cDNA. Proc Natl Acad Sci. 2018;115(39). 10.1073/pnas.1806447115. PMC6166824 30201725 · DOI
- [3]Fan X, Tang D, Liao Y, Li P, Zhang Y, Wang M, et al. Single-cell RNA-seq analysis of mouse preimplantation embryos by third-generation sequencing. PLoS Biol. 2020;18(12). 10.1371/journal.pbio.3001017. PMC7773192 33378329 · DOI
- [5]Wang Z Grabowski PJ Cell-and stage-specific splicing events resolved in specialized neurons of the rat cerebellum Rna 1996 2 12 1241 1253 8972773 PMC1369451
- [6]Smith MA, O’Dowd DK. Cell-specific regulation of agrin RNA splicing in the chick ciliary ganglion. Neuron. 1994;12(4). 10.1016/0896-6273(94)90332-8. 8161451 · DOI
- [7]Zipursky SL Wojtowicz WM Hattori D Got diversity? Wiring the fly brain with Dscam Trends Biochem Sci 2006 31 10 581 588 10.1016/j.tibs.2006.08.003 16919957 · DOI
- [8]Frank M Kemler R Protocadherins Curr Opin Cell Biol 2002 14 5 557 562 10.1016/S0955-0674(02)00365-4 12231349 · DOI
- [9]Kwan T Benovoy D Dias C Gurd S Provencher C Beaulieu P Genome-wide analysis of transcript isoform variation in humans Nat Genet 2008 40 2 225 231 10.1038/ng.2007.57 18193047 · DOI
- [10]Osborne RJ Thornton CA RNA-dominant diseases Hum Mol Genet 2006 15 2 R162 R169 10.1093/hmg/ddl181 16987879 · DOI
- [11]Philips AV, Timchenko LT, Cooper TA. Disruption of splicing regulated by a CUG-binding protein in myotonic dystrophy. Science. 1998;280(5364). 10.1126/science.280.5364.737. 9563950 · DOI
- [13]Dos Santos M, Backer S, Saintpierre B, Izac B, Andrieu M, Letourneur F, et al. Single-nucleus RNA-seq and FISH identify coordinated transcriptional activity in mammalian myofibers. Nat Commun. 2020;11(1). 10.1038/s41467-020-18789-8. PMC7547110 33037211 · DOI · Paper page
- [14]Petrany MJ, Swoboda CO, Sun C, Chetal K, Chen X, Weirauch MT, et al. Single-nucleus RNA-seq identifies transcriptional heterogeneity in multinucleated skeletal myofibers. Nat Commun. 2020;11(1). 10.1038/s41467-020-20063-w. PMC7733460 33311464 · DOI
- [15]Cornelison DDW. Context matters: In vivo and in vitro influences on muscle satellite cell activity. J Cell Biochem. 2008;105(3). 10.1002/jcb.21892. PMC3577053 18759329 · DOI
- [16]Trapnell C, Cacchiarelli D, Grimsby J, Pokharel P, Li S, Morse M, et al. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol. 2014;32(4). 10.1038/nbt.2859. PMC4122333 24658644 · DOI · Paper page
- [17]Zeng W, Jiang S, Kong X, El-Ali N, Ball AR, Ma CI-H, et al. Single-nucleus RNA-seq of differentiating human myoblasts reveals the extent of fate heterogeneity. Nucleic Acids Res. 2016. 10.1093/nar/gkw739. PMC5137429 27566152 · DOI
- [18]Olguin HC, Yang Z, Tapscott SJ, Olwin BB. Reciprocal inhibition between Pax7 and muscle regulatory factors modulates myogenic cell fate determination. J Cell Biol. 2007;177(5). 10.1083/jcb.200608122. PMC2064278 17548510 · DOI
- [19]Motohashi N, Asakura A. Muscle satellite cell heterogeneity and self-renewal. Front Cell Dev Biol. 2014;2. 10.3389/fcell.2014.00001. PMC4206996 25364710 · DOI
- [20]Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5). 10.1038/nbt.1621. PMC3146043 20436464 · DOI · Paper page
- [21]Harada Y, Nakamura M, Asano A. Temporally distinctive changes of alternative splicing patterns during myogenic differentiation of C2C12 cells. J Biochem. 1995;118(4). 10.1093/oxfordjournals.jbchem.a124980. 8576093 · DOI
- [22]Rosenberg AB Roco CM Muscat RA Kuchina A Sample P Yao Z Graybuck LT Peeler DJ Mukherjee S Chen W Pun SH Sellers DL Tasic B Seelig G Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding Science 2018 360 6385 176 182 10.1126/science.aam8999 29545511 PMC7643870 · DOI · Paper page
- [23]Wyman D Balderrama-Gutierrez G Reese F Jiang S Rahmanian S Zeng W Williams B Trout D England W Chu S Spitale RC Tenner A Wold B Mortazavi A A technology-agnostic long-read analysis pipeline for transcriptome discovery and quantification Biorxiv 2019 1 672931
- [24]Tardaguila M, de la Fuente L, Marti C, Pereira C, Pardo-Palacios FJ, del Risco H, et al. SQANTI: extensive characterization of long-read transcript sequences for quality control in full-length transcriptome identification and quantification. Genome Res. 2018;28(3). 10.1101/gr.222976.117. PMC5848618 29440222 · DOI
- [25]La Manno G Soldatov R Zeisel A Braun E Hochgerner H Petukhov V Lidschreiber K Kastriti ME Lönnerberg P Furlan A Fan J Borm LE Liu Z van Bruggen D Guo J He X Barker R Sundström E Castelo-Branco G RNA velocity of single cells Nature 2018 560 7719 494 498 10.1038/s41586-018-0414-6 30089906 PMC6130801 · DOI · Paper page
- [26]Reese, F., Mortazavi, A. (2020). Swan: a library for the analysis and visualization of long-read transcriptomes. Bioinformatics, btaa836. 10.1093/bioinformatics/btaa836 PMC8189675 32991665 · DOI · Paper page
- [27]Neganova I Vilella F Atkinson SP Lloret M Passos JF Von Zglinicki T An important role for CDK2 in G1 to S checkpoint activation and DNA damage response in human embryonic stem cells Stem Cells 2011 29 4 651 659 10.1002/stem.620 21319273 · DOI
- [28]Strzalka W Ziemienowicz A Proliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation Ann Bot 2011 107 7 1127 1140 10.1093/aob/mcq243 21169293 PMC3091797 · DOI
- [29]Chen T Sun Y Ji P Kopetz S Zhang W Topoisomerase IIα in chromosome instability and personalized cancer therapy Oncogene 2015 34 31 4019 4031 10.1038/onc.2014.332 25328138 PMC4404185 · DOI
- [30]Ren H Yin P Duan C IGFBP-5 regulates muscle cell differentiation by binding to IGF-II and switching on the IGF-II auto-regulation loop J Cell Biol 2008 182 5 979 991 10.1083/jcb.200712110 18762576 PMC2528583 · DOI
- [31]Rossi G, Antonini S, Bonfanti C, Monteverde S, Vezzali C, Tajbakhsh S, et al. Nfix regulates temporal progression of muscle regeneration through modulation of myostatin expression. Cell Rep. 2016;14(9). 10.1016/j.celrep.2016.02.014. PMC4793149 26923583 · DOI
- [33]Thomas K Engler AJ Meyer GA Extracellular matrix regulation in the muscle satellite cell niche Connect Tissue Res 2015 56 1 1 8 10.3109/03008207.2014.947369 25047058 PMC4464813 · DOI
- [34]Choi M, Ryu S, Hao R, Wang B, Kapur M, Fan C, et al. HDAC4 promotes Pax7-dependent satellite cell activation and muscle regeneration. EMBO Rep. 2014;15(11). 10.15252/embr.201439195. PMC4253491 25205686 · DOI
- [35]Oprescu SN, Yue F, Qiu J, Brito LF, Kuang S. Temporal dynamics and heterogeneity of cell populations during skeletal muscle regeneration. Science. 2020;23(4). 10.1016/j.isci.2020.100993. PMC7125354 32248062 · DOI
- [36]Yoshida T MCAT elements and the TEF-1 family of transcription factors in muscle development and disease Arterioscler Thromb Vasc Biol 2008 28 1 8 17 10.1161/ATVBAHA.107.155788 17962623 · DOI
- [37]Stuart T, Srivastava A, Lareau C, Satija R. Multimodal single-cell chromatin analysis with Signac. Biorxiv. 2020. 10.1101/2020.11.09.373613. · DOI
- [38]Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive integration of single-cell data. Cell. 2019;177(7). 10.1016/j.cell.2019.05.031. PMC6687398 31178118 · DOI · Paper page
- [39]Wei B, Jin J-P. TNNT1, TNNT2, and TNNT3: isoform genes, regulation, and structure–function relationships. Gene. 2016;582(1). 10.1016/j.gene.2016.01.006. PMC5325693 26774798 · DOI
- [40]Hestand MS, Klingenhoff A, Scherf M, Ariyurek Y, Ramos Y, van Workum W, et al. Tissue-specific transcript annotation and expression profiling with complementary next-generation sequencing technologies. Nucleic Acids Res. 2010;38(16). 10.1093/nar/gkq602. PMC2938216 20615900 · DOI
- [41]The ENCODE Project Consortium, Moore JE, Purcaro MJ, Pratt HE, Epstein CB, Shoresh N, et al. Expanded encyclopaedias of DNA elements in the human and mouse genomes. Nature. 2020;583(7818). 10.1038/s41586-020-2493-4. PMC7410828 32728249 · DOI · Paper page
- [42]Kivelä R, Salmela I, Nguyen YH, Petrova TV, Koistinen HA, Wiener Z, et al. The transcription factor Prox1 is essential for satellite cell differentiation and muscle fibre-type regulation. Nat Commun. 2016;7(1). 10.1038/ncomms13124. PMC5064023 27731315 · DOI
- [43]Feng X, Wang Z, Wang F, Lu T, Xu J, Ma X, et al. Dual function of VGLL4 in muscle regeneration. EMBO J. 2019;38(17). 10.15252/embj.2018101051. PMC6717915 31328806 · DOI
- [49]Li H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018;34(18). 10.1093/bioinformatics/bty191. PMC6137996 29750242 · DOI · Paper page
- [50]Wyman D, Mortazavi A. TranscriptClean: variant-aware correction of indels, mismatches and splice junctions in long-read transcripts. Bioinformatics. 2019;35(2). 10.1093/bioinformatics/bty483. PMC6329999 29912287 · DOI
- [53]Wolf FA, Angerer P, Theis FJ. SCANPY: large-scale single-cell gene expression data analysis. Genome Biol. 2018;19(1). 10.1186/s13059-017-1382-0. PMC5802054 29409532 · DOI · Paper page
- [54]Dobin A Davis CA Schlesinger F Drenkow J Zaleski C Jha S Batut P Chaisson M Gingeras TR STAR: ultrafast universal RNA-seq aligner Bioinformatics 2013 29 1 15 21 10.1093/bioinformatics/bts635 23104886 PMC3530905 · DOI · Paper page
- [55]Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33(5). 10.1038/nbt.3192. PMC4430369 25867923 · DOI · Paper page
- [58]Traag VA, Waltman L, van Eck NJ. From Louvain to Leiden: guaranteeing well-connected communities. Sci Rep. 2019;9(1). 10.1038/s41598-019-41695-z. PMC6435756 30914743 · DOI · Paper page
- [60]Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25(14). 10.1093/bioinformatics/btp324. PMC2705234 19451168 · DOI · Paper page
- [61]Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based Analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9). 10.1186/gb-2008-9-9-r137. PMC2592715 18798982 · DOI · Paper page
- [62]Yu W, Uzun Y, Zhu Q, Chen C, Tan K. scATAC-pro: a comprehensive workbench for single-cell chromatin accessibility sequencing data. Genome Biol. 2020;21(1). 10.1186/s13059-020-02008-0. PMC7169039 32312293 · DOI
- [63]Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42(W1). 10.1093/nar/gku365. PMC4086134 24799436 · DOI · Paper page
- [64]McLean CY Bristor D Hiller M Clarke SL Schaar BT Lowe CB GREAT improves functional interpretation of cis-regulatory regions Nat Biotechnol 2010 28 5 495 501 10.1038/nbt.1630 20436461 PMC4840234 · DOI · Paper page
- [66]Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6). 10.1093/bioinformatics/btq033. PMC2832824 20110278 · DOI · Paper page
- [67]Martin M Cutadapt removes adapter sequences from high-throughput sequencing reads EMBnet J 2011 17 1 10 12 10.14806/ej.17.1.200 · DOI · Paper page
- [68]Hestand MS, ‘t Hoen PA. DeepCAGE and DeepSAGE with proliferating and differentiated C2C12 mouse myoblasts. GEO. 2010. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE21580. Accessed Sept 2021.
- [69]Neph S Kuehn MS Reynolds AP Haugen E Thurman RE Johnson AK Rynes E Maurano MT Vierstra J Thomas S Sandstrom R Humbert R Stamatoyannopoulos JA BEDOPS: high-performance genomic feature operations Bioinformatics 2012 28 14 1919 1920 10.1093/bioinformatics/bts277 22576172 PMC3389768 · DOI
- [70]Li W, You B, Hoque M, Zheng D, Luo W, Ji Z, et al. Systematic profiling of poly(A) + transcripts modulated by core 3’ end processing and splicing factors reveals regulatory rules of alternative cleavage and polyadenylation. PLoS Genet. 2015;11(4). 10.1371/journal.pgen.1005166. PMC4407891 25906188 · DOI
- [71]Li W, You B, Zheng D, Hoque M, Luo W, Ji Z, et al. Regulation of alternative cleavage and polyadenylation by 3’ end processing and splicing factors. GEO. 2015. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE62001. Accessed Sept 2021.
- [72]Hinrichs AS, Karolchik D, Baertsch R, Barber GP, Bejerano G, Clawson H, et al. The UCSC Genome Browser Database: update 2006. Nucleic Acids Res. 2006;34(90001). 10.1093/nar/gkj144. PMC1347506 16381938 · DOI · Paper page
- [73]Frankish A, Diekhans M, Ferreira A-M, Johnson R, Jungreis I, Loveland J, et al. GENCODE reference annotation for the human and mouse genomes. Nucleic Acids Res. 2019;47(D1). 10.1093/nar/gky955. PMC6323946 30357393 · DOI
Cited by
Sorted by citations
hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data
Morabito, S., Reese, F., Rahimzadeh, N., Miyoshi, E., & Swarup, V. (2023) · Cell Reports Methods · doi:10.1016/j.crmeth.2023.100498 · Cited by 773 · see citing works
Self-citationCardiac splicing as a diagnostic and therapeutic target
Gotthardt, M., Badillo-Lisakowski, V., Parikh, V. N., Ashley, E., Furtado, M., Carmo-Fonseca, M., Schudy, S., Meder, B., Grosch, M., Steinmetz, L., et al. (2023) · Nature Reviews Cardiology · doi:10.1038/s41569-022-00828-0 · Cited by 46 · see citing works
- Authors
- Rebboah, E., Reese, F., Williams, K., Balderrama-Gutierrez, G., McGill, C., Trout, D., Rodriguez, I., Liang, H., Wold, B. J., & Mortazavi, A.
- Affiliations
- University of California, Irvine (US)
- California Institute of Technology (US)
- Corresponding author
- Mortazavi, Ali
- ORCID
- Venue
- Genome Biology
- Published
- 2021-10-07
- Volume
- 22
- Issue
- 1
- Pages
- 286-286
- ISSN
- 1465-6906 (ISSN-L)
- 1474-760X
- PMC
- PMC8495978
- OpenAlex
- W3202223582
- Cited by
- 58
- References
- 8
- Citations by year
- 2026: 12 · 2025: 10 · 2024: 17 · 2023: 11 · 2022: 7 · 2021: 1
- FWCI
- 3.14
- Type
- journal article
- Language
- en
- Open access
- gold
- doi.org
- Locations
- Genome biology (publishedVersion, cc-by)
- bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
- Europe PMC (PubMed Central) (publishedVersion, other-oa)
- PubMed Central (publishedVersion, cc-by)
- PubMed (publishedVersion)
- eScholarship (California Digital Library) (submittedVersion, cc-by)
- CaltechAUTHORS (California Institute of Technology) (submittedVersion, cc-by)
- DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
- Funders
- National Human Genome Research Institute (UM1HG009443)
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Research and Splicing
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
- Keywords
- combinatorial barcoding, alternative transcription start sites, alternative exon usage, C2C12 myoblasts, cell differentiation, single-cell genomics
- MeSH terms
- RNA-Seq, Animals, Cell Differentiation, Cell Line, Cell Nucleus, Chromatin, Models, Genetic, Transcription, Genetic, Myogenin, Genomics, Transcription Initiation Site, Mice, PAX7 Transcription Factor, Single-Cell Analysis, RNA Isoforms
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review