RNA isoform-resolved multiplexed sequencing with bioorthogonal barcoding

RNA isoform dysregulation drives disease pathogenesis and is the target of FDA-approved splice-switching therapeutics. However, multiplexed sequencing methods discard splice junction information because only 3' termini are barcoded and counted. Here, we repurpose acylation and click chemistries to conjugate bioorthogonal barcodes (bobcodes) directly onto multiple internal positions along cellular RNAs. Bobcoded RNAs from multiple samples are pooled for multiplexed cDNA synthesis, during which reverse transcriptase switches from each RNA template onto its tethered bobcode with greater than 99% accuracy in species mixing experiments. Bobcode attachment intervals set cDNA insert sizes without a library fragmentation step, and priming with poly(dT) or random hexamers selects between 3'-end counting and full-length isoform capture. A bioorthogonal barcode-sequencing (BOB-seq v0.1) drug screen identifies transcriptome-wide on- and off-target RNA splicing effects and outperforms existing multiplexing RNA sequencing methods in workflow simplicity, sample-to-sample variability, and barcoding accuracy. Bobcodes add isoform resolution to scalable multiplexed RNA sequencing.