Answers.Body: <blockquote>
    <p>I was hoping to find some tool that already exists.</p>
  </blockquote>
  
  <p>The short answer is, <strong>none exist yet</strong>. I see a number of hurdles to overcome first.<sup>1</sup>  </p>
  
  <p>For example, <strong>no one has yet characterized saturation for transcript discovery</strong>, i.e., the reads per cell needed to maximize transcript information recovery.  This is a typical stepping-stone as a sequencing approach matures, and until it's done, it makes it unlikely that most scRNA-seq datasets are going to be appropriate vis-à-vis statistical power, even if there was a tool.  That is, if people are only aiming for cell-type discrimination via gene-level markers, datasets will likely be underpowered to measure transcript level differences.</p>
  
  <p>Another issue is the <strong>lack of any transcript-level scRNA-seq simulators</strong>. All the simulators out there that I know of produce gene-level matrices.  Having a transcript-level simulator is crucial for benchmarking and characterizing the performance of differential expression calling methodologies.  Without anything to test against, that leaves a huge blindspot in terms of estimating the uncertainty (false positive rates, sensitivity) that any tool might have.  However, since some scRNA-seq models (e.g., <a href="https://www.ncbi.nlm.nih.gov/pubmed/26040288" rel="nofollow noreferrer">Velten, et al. 2015</a>, or <a href="https://www.ncbi.nlm.nih.gov/pubmed/30504886" rel="nofollow noreferrer">Lopez, et al., 2018</a>) are generative, perhaps this is shouldn't be too difficult to adapt to generating FASTQs.</p>
  
  <p>Simulators are fine and all, but ground truth is probably the key hurdle of the enterprise right now. While there are definitely groups using smFISH to distinguish isoforms quantitatively on a single-cell level, <strong>there doesn't appear to be enough smFISH data out there to measure a transcriptome-wide scRNA-seq method against</strong>. Hopefully this will change in the coming year or two, as results from more highly-multiplexed single-molecule technologies become available.</p>
  
  <p><sup><strong>1</strong> This is just my opinion formed as I try to keep up with the literature, so I could be missing things.  Comments/corrections certainly welcome!</sup></p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-01-28T20:16:15.203Z
Answers.LastActivityDate: 2019-01-28T20:16:15.203Z
Answers.UserDisplayName: merv
Answers.Body: <p>There is a paper that takes advantage of transcript compatibility counts (TCCs). A TCC describes all of the potential transcripts that a cDNA alignment can be derived from.  The senior author, Lior Pachter also has a nice <a href="https://liorpachter.wordpress.com/2018/02/15/gde%C2%B2-dge%C2%B2-dtu%C2%B2-dte%E2%82%81%C2%B2-dte%E2%82%82%C2%B2/" rel="nofollow noreferrer">blog post</a> about differential transcript usage.</p>
  
  <p><a href="https://www.nature.com/articles/s41592-018-0303-9" rel="nofollow noreferrer">A discriminative learning approach to differential expression analysis for single-cell RNA-seq</a></p>
  
  <blockquote>
    <p>Single-cell RNA-seq makes it possible to characterize the
    transcriptomes of cell types across different conditions and to
    identify their transcriptional signatures via differential analysis.
    Our method detects changes in transcript dynamics and in overall gene
    abundance in large numbers of cells to determine differential
    expression. When applied to transcript compatibility counts obtained
    via pseudoalignment, our approach provides a quantification-free
    analysis of 3′ single-cell RNA-seq that can identify previously
    undetectable marker genes.</p>
  </blockquote>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-01-24T23:26:38.103Z
Answers.Text: That's a nice paper. Unfortunately, it's more of a proof-of-concept. I was hoping to find some tool that already exists.
Answers.UserDisplayName: burger
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2019-01-24T15:20:56.440Z
Answers.LastActivityDate: 2019-01-24T15:20:56.440Z
Answers.UserDisplayName: GWW
Body: <p>Many of the modern gene-quantification tools (Salmon/Kallisto) output transcript-level (as opposed to gene-level) data. All of the scRNA-seq analysis I have seen just uses the gene-level values. I suppose I can just use the transcript-level matrix, but then I lose any links between the different isoforms of the same gene. Maybe it does not really matter?</p>
  
  <p>There are a few options for bulk RNA-seq and there was a <a href="https://bioinformatics.stackexchange.com/questions/774/analysis-of-differential-transcript-usage-dtu">similar question here</a>, but are there any single-cell options?</p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-01-23T17:43:24.120Z
Comments.Text: Single cell data tends to be very spotty, I don't know that you'll have enough reads to do a transcript level analysis from a typical sc-RNA experiment.
Comments.UserDisplayName: swbarnes2
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-01-23T18:23:27.770Z
Comments.Text: @swbarnes2 True, but it's spotty even on gene level.
Comments.UserDisplayName: burger
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-01-24T14:39:10.190Z
Comments.Text: Could you please clarify the question? You want to know if using gene-level values in scRNA is fine with later methods or if it is fine to use transcript level values with later analysis? Or do you want to perform differential transcript usage and you don't know if it is fine (with what? other methods the conclusions the data?...) ?
Comments.UserDisplayName: llrs
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-01-24T23:25:13.653Z
Comments.Text: @llrs I was trying to leave the question a little open-ended. The real question is at the end: what are DTU options for scRNA-seq?
Comments.UserDisplayName: burger
Tags: <scrnaseq><single-cell><transcriptome>
Title: scRNA-seq differential transcript usage
ThreadText: # scRNA-seq differential transcript usage
  
  (Asked by: username_0 on 2019-01-23)
  
  Many of the modern gene-quantification tools (Salmon/Kallisto) output transcript-level (as opposed to gene-level) data. All of the scRNA-seq analysis I have seen just uses the gene-level values. I suppose I can just use the transcript-level matrix, but then I lose any links between the different isoforms of the same gene. Maybe it does not really matter?
  
  There are a few options for bulk RNA-seq and there was a [similar question here](https://bioinformatics.stackexchange.com/questions/774/analysis-of-differential-transcript-usage-dtu), but are there any single-cell options?
  
  **Tags:** scrnaseq, single-cell, transcriptome
  
  > **Comments:**
  > **username_1**: Single cell data tends to be very spotty, I don't know that you'll have enough reads to do a transcript level analysis from a typical sc-RNA experiment.  
  > **username_0**: @username_1 True, but it's spotty even on gene level.  
  > **username_2**: Could you please clarify the question? You want to know if using gene-level values in scRNA is fine with later methods or if it is fine to use transcript level values with later analysis? Or do you want to perform differential transcript usage and you don't know if it is fine (with what? other methods the conclusions the data?...) ?  
  > **username_0**: @username_2 I was trying to leave the question a little open-ended. The real question is at the end: what are DTU options for scRNA-seq?  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_3 on 2019-01-28. Score: 4
  
  > I was hoping to find some tool that already exists.
  
  The short answer is, **none exist yet**. I see a number of hurdles to overcome first.1
  
  For example, **no one has yet characterized saturation for transcript discovery**, i.e., the reads per cell needed to maximize transcript information recovery. This is a typical stepping-stone as a sequencing approach matures, and until it's done, it makes it unlikely that most scRNA-seq datasets are going to be appropriate vis-à-vis statistical power, even if there was a tool. That is, if people are only aiming for cell-type discrimination via gene-level markers, datasets will likely be underpowered to measure transcript level differences.
  
  Another issue is the **lack of any transcript-level scRNA-seq simulators**. All the simulators out there that I know of produce gene-level matrices. Having a transcript-level simulator is crucial for benchmarking and characterizing the performance of differential expression calling methodologies. Without anything to test against, that leaves a huge blindspot in terms of estimating the uncertainty (false positive rates, sensitivity) that any tool might have. However, since some scRNA-seq models (e.g., [Velten, et al. 2015](https://www.ncbi.nlm.nih.gov/pubmed/26040288), or [Lopez, et al., 2018](https://www.ncbi.nlm.nih.gov/pubmed/30504886)) are generative, perhaps this is shouldn't be too difficult to adapt to generating FASTQs.
  
  Simulators are fine and all, but ground truth is probably the key hurdle of the enterprise right now. While there are definitely groups using smFISH to distinguish isoforms quantitatively on a single-cell level, **there doesn't appear to be enough smFISH data out there to measure a transcriptome-wide scRNA-seq method against**. Hopefully this will change in the coming year or two, as results from more highly-multiplexed single-molecule technologies become available.
  
  **1** This is just my opinion formed as I try to keep up with the literature, so I could be missing things. Comments/corrections certainly welcome!
  
  
  ---
  
  ### Answer by username_4 on 2019-01-24. Score: 2
  
  There is a paper that takes advantage of transcript compatibility counts (TCCs). A TCC describes all of the potential transcripts that a cDNA alignment can be derived from. The senior author, Lior Pachter also has a nice [blog post](https://liorpachter.wordpress.com/2018/02/15/gde%C2%B2-dge%C2%B2-dtu%C2%B2-dte%E2%82%81%C2%B2-dte%E2%82%82%C2%B2/) about differential transcript usage.
  
  [A discriminative learning approach to differential expression analysis for single-cell RNA-seq](https://www.nature.com/articles/s41592-018-0303-9)
  
  > Single-cell RNA-seq makes it possible to characterize the
  > transcriptomes of cell types across different conditions and to
  > identify their transcriptional signatures via differential analysis.
  > Our method detects changes in transcript dynamics and in overall gene
  > abundance in large numbers of cells to determine differential
  > expression. When applied to transcript compatibility counts obtained
  > via pseudoalignment, our approach provides a quantification-free
  > analysis of 3′ single-cell RNA-seq that can identify previously
  > undetectable marker genes.
  
  > **Comments:**
  > **username_0**: That's a nice paper. Unfortunately, it's more of a proof-of-concept. I was hoping to find some tool that already exists.  
  
  ---
  
  