<p>The answer is behind these papers:<br>
<a href="http://science.sciencemag.org/content/early/2018/04/25/science.aar5780" rel="nofollow noreferrer">The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution</a><br>
 <a href="http://science.sciencemag.org/content/early/2018/04/25/science.aar4362.full" rel="nofollow noreferrer">Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo</a></p>

<p>Also <a href="http://science.sciencemag.org/content/sci/suppl/2018/04/25/science.aar3131.DC1/aar3131-Farrell-SM.pdf" rel="nofollow noreferrer">this</a>.</p>

CC BY-SA 4.0
2018-09-06T05:51:13.963Z
Please avoid link only answers, provide summary text from the links.
zx8754
CC BY-SA 4.0
2018-09-07T10:38:08.953Z
The papers are too difficult to follow but I just resealized that finding parent child clusters likely would work.
Zizogolu
CC BY-SA 4.0
2018-09-03T08:49:18.620Z
2018-09-30T14:22:31.487Z
2018-09-30T14:22:31.487Z
Zizogolu
<p>This is what packages like <a href="https://www.ncbi.nlm.nih.gov/pubmed/29630061" rel="nofollow noreferrer">FateID</a> and <a href="http://cole-trapnell-lab.github.io/monocle-release/docs/" rel="nofollow noreferrer">Monocle</a> are for, namely taking single-cell RNAseq data and inferring differentiation trajectories from it. Don't try to reinvent the wheel on this, there are a number of packages out there to do this sort of thing and they're going to have pretty complicated methods that you're really not going to want to reinvent unless you absolutely have to. </p>

CC BY-SA 4.0
2018-08-23T07:02:26.960Z
2018-08-23T07:02:26.960Z
Devon Ryan
<p>I have single cell RNA-seq data on about 2000 cells in 9 time point. I have clustered my cells in each time point by Seurat. I am seeing in some time points I have 3 clusters while in another time points I have 2 clusters of cells (please look at the plot).</p>
<p>One hypothesis would be that, a new cell type is being emerged at last time point. But maybe if we had another further time point this third cluster being disappeared because that has been arose just because of undifferentiated cells from the later time point (second hypothesis).<br />
Another hypothesis would be this new third cluster is a new cells type that only temporarily expresses a set of genes but if we had cells for another further time point, we would not see the expression of these genes anymore and we had again 2 clusters. I am really struggling with this how I could test my hypothesises.<br />
The final hypothesis would be that the third cluster is really a new cell type born from old cell types during the differentiation.</p>
<p>Like UAD algorithm addresses my question. I just tried FateID algorithm on my data but got stopped by</p>
<pre><code>&gt; pr &lt;- prcurve(y,fb,dr,k=2,m=&quot;tsne&quot;,trthr=0.4,start=2)
Error in svd(xstar) : a dimension is zero
</code></pre>
<p>Any suggestion would be highly appreciated.</p>

<r><scrnaseq><seurat>
How I can test my hypothesises computationally
# How I can test my hypothesises computationally

(Asked by: username_0 on 2018-08-22)

I have single cell RNA-seq data on about 2000 cells in 9 time point. I have clustered my cells in each time point by Seurat. I am seeing in some time points I have 3 clusters while in another time points I have 2 clusters of cells (please look at the plot).
One hypothesis would be that, a new cell type is being emerged at last time point. But maybe if we had another further time point this third cluster being disappeared because that has been arose just because of undifferentiated cells from the later time point (second hypothesis).
Another hypothesis would be this new third cluster is a new cells type that only temporarily expresses a set of genes but if we had cells for another further time point, we would not see the expression of these genes anymore and we had again 2 clusters. I am really struggling with this how I could test my hypothesises.
The final hypothesis would be that the third cluster is really a new cell type born from old cell types during the differentiation.
Like UAD algorithm addresses my question. I just tried FateID algorithm on my data but got stopped by
> pr <- prcurve(y,fb,dr,k=2,m="tsne",trthr=0.4,start=2)
Error in svd(xstar) : a dimension is zero

Any suggestion would be highly appreciated.

**Tags:** r, scrnaseq, seurat

---

## 2 answers

### Answer by username_0 on 2018-09-03. Score: 1 | Accepted answer

The answer is behind these papers:  
[The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution](http://science.sciencemag.org/content/early/2018/04/25/science.aar5780)  
[Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo](http://science.sciencemag.org/content/early/2018/04/25/science.aar4362.full)

Also [this](http://science.sciencemag.org/content/sci/suppl/2018/04/25/science.aar3131.DC1/aar3131-Farrell-SM.pdf).

> **Comments:**
> **username_1**: Please avoid link only answers, provide summary text from the links.  
> **username_0**: The papers are too difficult to follow but I just resealized that finding parent child clusters likely would work.  

---

### Answer by username_2 on 2018-08-23. Score: 3

This is what packages like [FateID](https://www.ncbi.nlm.nih.gov/pubmed/29630061) and [Monocle](http://cole-trapnell-lab.github.io/monocle-release/docs/) are for, namely taking single-cell RNAseq data and inferring differentiation trajectories from it. Don't try to reinvent the wheel on this, there are a number of packages out there to do this sort of thing and they're going to have pretty complicated methods that you're really not going to want to reinvent unless you absolutely have to.


---

