Question: Is there a way to compare single-cell RNA seq and Ribotag-seq?
gravatar for Jon
21 days ago by
United States
Jon90 wrote:

Hello there,

I have two different dataset(different library), one is single cell RNA sequencing of 60 cells. And other one is 10 Ribotag IP sequencing (traslatome profiling). Is there any way to compare and correlate these two datasets?

When I try to plot PCA they cluster at different places but both datasets are belong to same cell types.

Thanks Jon

ADD COMMENTlink modified 12 days ago • written 21 days ago by Jon90

You'll have to think about what you really want to get out of the data. Most likely you'll want to look at how DE genes in your Ribotag data change in different cells in your scRNA-seq data (with only 60 cells, I expect you explicitly took different cell types).

ADD REPLYlink written 12 days ago by Devon Ryan89k

Hi Devon,

I plotted the cell trajectory using monocle. When I tried to plot, I used different different datasets of single cell seq and my ribotag seq data. I used normalizebetweenarrays function from limma, to normalize the datasets, to have same column sums. All the single cells are looking good, that they are staying in acceptable positions in the trajectory, but my ribotag seq data stays at the end of trajectory, but they must be in the middle of trajectory(based on their cell maturation level).

And how can I do this

how DE genes in your Ribotag data change in different cells in your scRNA-seq data

ADD REPLYlink modified 11 days ago • written 11 days ago by Jon90

Combining those two datasets in a trajectory won't produce sensible results.

ADD REPLYlink written 11 days ago by Devon Ryan89k

Okay, as the single cell datasets have the information of maturation. Can i show particular ribotag-seq belongs to particular maturation level? Is it possible?

ADD REPLYlink written 10 days ago by Jon90

That's possible but quite difficult. Your setup doesn't lend itself to directly comparing the two types of datasets. What you might be able to do is look at genes changing around the time point you want in the trajectory and then look at how they're expressed in your Ribotag data. Reviewers aren't going to have plenty of objections to that, but if the results fit in with a broader story that sort of procedure should work.

ADD REPLYlink written 10 days ago by Devon Ryan89k

okay, Can you tell me some examples/R packages/algorithms to do.


ADD REPLYlink written 9 days ago by Jon90

You'll have to invent this.

ADD REPLYlink written 8 days ago by Devon Ryan89k

That's exciting, thanks for your discussions.

ADD REPLYlink written 8 days ago by Jon90

Hi Devon, I got Dynamically regulated genes across the maturation trajectory(for single cell data), around 1000 genes, with 24 pseudo-time points.

I have calculated correlation (spearman's/Pearson) between each ribotag data and all the pseudotime points. (to find which pseudo time expression pattern is similar to ribotag data)

All the ribotag data highly correlated with more mature state. but in reality, they are not.

Is it a right way of doing? do you have any suggestions?

ADD REPLYlink written 3 days ago by Jon90

That's probably as good as can be done. You might have to tweak filtering to get things to align a bit better with your background knowledge of how things should look.

ADD REPLYlink written 2 hours ago by Devon Ryan89k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1505 users visited in the last hour