Multiplexed single-cell RNA-seq analysis of multiple samples using pooling is a promising experimental design, offering increased throughput while allowing to overcome batch variation. To reconstruct the sample identify of each cell, genetic variants that segregate between the samples in the pool have been proposed as natural barcode for cell demultiplexing. Existing demultiplexing strategies rely on availability of complete genotype data from the pooled samples, which limits the applicability of such methods, in particular when genetic variation is not the primary object of study. To address this, we here present Vireo, a computationally efficient Bayesian model to demultiplex single-cell data from pooled experimental designs. Uniquely, our model can be applied in settings when only partial or no genotype information is available. Using pools based on synthetic mixtures and results on real data, we demonstrate the robustness of Vireo and illustrate the utility of multiplexed experimental designs for common expression analyses.

Site-specific gene addition can allow stable transgene expression for gene therapy. When possible, this is preferred over the use of promiscuously integrating vectors, which are sometimes associated with clonal expansion1 and oncogenesis2. Site-specific endonucleases that can induce high rates of targeted genome editing are finding increasing applications in biological discovery and gene therapy3. However, two safety concerns persist: (1) endonuclease-associated adverse effects, both on4 and off-target5,6; and (2) oncogene activation caused by promoter integration, even without nucleases7. Here, we perform recombinant adeno-associated virus (rAAV) mediated promoterless gene targeting without nucleases and demonstrate amelioration of the bleeding diathesis in haemophilia B mice. In particular, we target a promoterless human coagulation factor IX (hF9) gene to the liver-expressed albumin (Alb) locus. hF9 is targeted, along with a preceding 2A-peptide coding sequence, to be integrated just upstream to the Alb stop codon. While hF9 is fused to Alb at the DNA and RNA levels, two separate proteins are synthesized by way of ribosomal skipping. Thus, hF9 expression is linked to robust hepatic albumin expression without disrupting it. We injected an AAV8-hF9 vector into neonatal and adult mice and achieved on-target integration into ~0.5% of the albumin alleles in hepatocytes. We established that hF9 was produced only from on-target integration, and ribosomal skipping was highly efficient. Stable hF9 plasma levels at 7–20% of normal were obtained, and treated factor IX deficient mice had normal coagulation times. In conclusion, transgene integration as a 2A-fusion to a highly expressed endogenous gene may obviate the requirement for nucleases and/or vector-borne promoters. This method may allow for safe and efficacious gene targeting in both infants and adults by greatly diminishing off-target effects while still providing therapeutic levels of expression from integration.

The X‐linked genetic bleeding disorder caused by deficiency of coagulator factor IX, hemophilia B, is a disease ideally suited for gene therapy with genome editing technology. Here, we identify a family with hemophilia B carrying a novel mutation, Y371D, in the human F9 gene. The CRISPR/Cas9 system was used to generate distinct genetically modified mouse models and confirmed that the novel Y371D mutation resulted in a more severe hemophilia B phenotype than the previously identified Y371S mutation. To develop therapeutic strategies targeting this mutation, we subsequently compared naked DNA constructs versus adenoviral vectors to deliver Cas9 components targeting the F9 Y371D mutation in adult mice. After treatment, hemophilia B mice receiving naked DNA constructs exhibited correction of over 0.56% of F9 alleles in hepatocytes, which was sufficient to restore hemostasis. In contrast, the adenoviral delivery system resulted in a higher corrective efficiency but no therapeutic effects due to severe hepatic toxicity. Our studies suggest that CRISPR/Cas‐mediated in situ genome editing could be a feasible therapeutic strategy for human hereditary diseases, although an efficient and clinically relevant delivery system is required for further clinical studies.

Single-cell RNA-seq (scRNA-seq) provides a comprehensive measurement of stochasticity in transcription, but the limitations of the technology have prevented its application to dissect variability in RNA processing events such as splicing. Here, we present BRIE (Bayesian regression for isoform estimation), a Bayesian hierarchical model that resolves these problems by learning an informative prior distribution from sequence features. We show that BRIE yields reproducible estimates of exon inclusion ratios in single cells and provides an effective tool for differential isoform quantification between scRNA-seq data sets. BRIE, therefore, expands the scope of scRNA-seq experiments to probe the stochasticity of RNA processing.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1248-5) contains supplementary material, which is available to authorized users.