Single-cell Trajectory analysis using Monocle3 and Seurat | Step-by-step tutorial

[@legttv97]

This channel will be a respectable resource in our bioinformatics community

[@kimayatekade-cu1ds]

If it hasnt been said enough already, I will repeat : what you are doing is AWESOME!

[@beyzaozdemir394]

You are making my learning journey much more easier for an intermediate use and making follow vignettes a breeze! Glad to have found your SEURAT videos!

[@zlj8435]

Thanks a lot! The class is easy to follow with the github code and great to know how to do that and why we need to do that. Really appreciate it!!

[@christinaz3573]

Honestly, this was so precisely what I needed in a previously very frustrating situation. The way you manage to explain the background of what you are doing connected to how it is done is absolutely awesome. Thank you so so much!

[@waynejanuzaj9679]

Your tutorials have been very timely, informative and helpful. We will appreciate it if you could also make a tutorial on SCENIC. Thank you so much!

[@Bioinformagician]

Will surely consider making a video on it! Thanks for the suggestion :)

[@jintaoxu6393]

This playlist is so helpful. Would it be possible to have a Step-by-step tutorial for cell-cell interaction?

[@junorzhu4313]

very detailed tutorial, I found that myself more susceptible to your pattern compaired to others

[@tautarius]

thanks for making these awesome videos. I really appreciate it !!!

[@yairprawer4405]

so helpful!! saved me hours of of reading error messages

[@jianwu4593]

Thanks

[@5602KK]

EXCELLENT video thank you very much. Excellent clear explanations accompanied with great demonstrations & slides I really mean that

[@Aquaria-]

Thank you so much! Really helpful for me to convert Seurat object for Monocle3

[@melinabordone8243]

First, thank you a lot. Your tutorial videos on scRNAseq analysis are a treasure! I would like to know if you already have a video on differential cell-type composition analysis or if you are planning to record one. Thank you so much!

[@babil1908]

I am learning so much from you, thank you so much!

[@SwedishRagers]

Hi! there is not tutorial on YT on combining microarray and rna sequencing data under one analysis. this would be a very helpful tutorial! thanks again for all your amazing content!

[@Bioinformagician]

Since microarray and RNA-Seq are different technologies, I am not sure of any use cases where there would be a requirement to combine data from both these technologies.

[@SwedishRagers]

@@Bioinformagician I was thinking of when one does a meta-analysis on data from GEO. Some studies have done this for example: doi: 10.3389/fgene.2021.663787 however, it is difficult to learn how to replicate this. Thank you though.

[@navyanandhanaofficial]

Very nice Presentation! I liked it... Appreciate you!

[@Phsjd.2736]

Thanks so much. Very informative videos!

[@cocnfkdjsnnnbb]

Excellent tutorial 👏

[@LazizAsamov-u9h]

you are the best. Thank you

[@abdullahaltulea142]

Excellent! Thank you very much

[@youvikasingh7955]

Hi @Bioinfomagician..Loved your tutorial, and much-needed channel for budding bioinformaticians. I wanted to know if we don't have information on root cluster then how do we order cells according to Pseudotime?For eg in your case ProB cells are from cluster 5 but In my case, I just have CD4T cells in 4 clusters...how do I order them? Thanks in advance🙏

[@Bioinformagician]

I am afraid pseudotime can be determined accurately for other cells if we don't know the cells it originated from i.e. the root cluster.

[@youvikasingh7955]

@@Bioinformagicianthanks for your reply😍 is there any strategy to find out root cluster? I was doing this like..looking for lymphoid differentiation markers and setting those clusters or cells as root I am not sure if this way is right I appreciate your help! Thanks a lot... you are doing a great job..whenever i am stuck somewhere i look into your channel for help😊

[@Bioinformagician]

@@youvikasingh7955 Unfortunately I am not aware of any reliable methods to find out root cluster. Your best bet is to dig up literature and research papers to find out whether any lineage studies have been done previously for the cells you are studying which can give you clues for the origin cell type and associated markers. I am glad to hear my channel is serving as a useful resource :)

[@rps2930]

Hi Khushbu, first of all I really appreciate the effort you put in sharing your knowledge with us. I want to know from where I can get dataset for trajectory analysis which can be further published?

[@Jenkins-f7s]

Thank you! I am not seeing any '^MT- genes showing up in the data seuratobj I made from the downloaded data. Just get error: Error in validObject(object = x) : invalid class “Seurat” object: 1: all cells in assays must be in the same order as the Seurat object invalid class “Seurat” object: 2: 'active.idents' must be named with cell names Thanks if you have seen this error before

[@louisgervais5059]

I have the same error message

[@xinyuqu4407]

I have the same issue. I use the AddMetaData function from Seurat to solve the problem. rownames(metadata)

[@abhimanyu--3741]

Thank you for the tutorial. Can you extend this to "plot_genes_in_pseudotime" after the "graph_test"? that is one confusing part of the code on their tutorial

[@ryanpark7323]

Thank you! So so helpful :)

[@tushardhyani3931]

Thank you for this video !!

[@Pato-py7wt]

Amazing clear helpful video

[@2007dolphinlover]

Why do you also want to compare the cluster before trajectory?

[@miguelcuevas2976]

What is the benefit of monocle over other trajectory analysis tools such as slingshot?

[@hassanbalawar5379]

You are doing a great job. What if the data is not annotated?

[@kiranpatil2924]

Thank you,, Please Make same video for Spatial Transcriptomics Analysis

[@mayaziv8088]

Thank you great videos!!

[@alexbrown2272]

Thank you so much for this great video! Curious if you have found a way to plot individual genes in pseudo time? There is a way in Monocle, but when using this integration with Seurat there is not an expression family that is needed to plot specific genes in pseudo time or to find which genes are driving the trajectory. Any advice would be greatly appreciated! Thanks again!

[@brooke2100]

I am wondering this as well. Thanks in advance!

[@Bioinformagician]

You can plot specific genes in pseudotime by running cds

[@meetukaur0909]

Ma'am, can you please let me know where can I get gene and cell meta data for a particular single cell analysis dataset?

[@qiongyujadehe5778]

Thanks alot !very useful!

[@potterae]

Your videos are awesome, very helpful! For setting the clusters in the cell data set, do you think it could be a viable option to use the defined subclusters or is it better to use the clusters which seurat found with FindClusters? Thank you!

[@prasadchaskar8542]

Thank you so much. Best!!!!

[@adelyntsai1116]

This is great!! Thank you so much. Will it be possible if you can make tutorials on intercellular interaction tools in R? Thank you again!

[@Bioinformagician]

Thanks for the suggestion, will surely considering making videos covering this topic.

[@msumode4493]

Hello Khushbu. 2 points I am confused about. 1. I was wondering when converting from Seurat to Monocle3 cds object, are the counts used for downstream analysis in your tutorial or the data from Seurat's data slot? Counts are not normalized and wouldn't that lead to erroneous calculations if not normalized in Monocle3 later on? I am not sure if this was done in this tutorial. 2. Why take Seurat's UMAP coordinates for analysis in Monocle3? Wouldn't it be more informative to know how Monocle3 projects the Seurat generated clusters in 2D space? Otherwise, isn't Monocle just drawing or projecting pseudotime over Seurat based proximity? Thank you.

[@lukesimpson1507]

Hi Bioinfomagician, great video as always. You spoke about the importance of having the optimal cluster resolution, is there any objective way to determine this? Is there any tools you would recommend? Thank you!

[@Bioinformagician]

What I usually do is run clustering with multiple clustering resolutions like this: FindClusters(object = seurat_integrated, resolution = c(0.2, 0.5, 0.8, 1.0)) and then plot each of them to visually see which resolutions provide optimal separation. If you have known cell types in your data, it might be of help to visualize the markers for those cell types too for different resolutions to make sure different cell types, group into separate clusters.

[@kitdordkhar4964]

Question? From where did you get the redefined_cluster? I always get confused with the term "Subset" in scRNA analysis.

[@Bioinformagician]

'redefined_cluster' was already present in the data I am using, as the data was already annotated. The term 'subset' is to just get certain group of cells from the entire dataset.

[@kitdordkhar4964]

@@Bioinformagician Gotcha! It was already in the ref database I guess. Thanks!

[@anaburgos4969]

Thank you SO MUCH :)

[@stacygenovese1761]

As always these are wonderful. Quick question - if we're doing a trajectory analysis following integration do we just run as.cell_data_set() on our existing Seurat object? Or do we have to re-create the object transforming the expression matrix?

[@Bioinformagician]

I think you could read the integrated object into as.cell_data_set() and use the clustering information from Seurat's UMAP to cluster cells and learn trajectory on those embeddings.

[@judep6077]

Make sure to change the active assay to 'RNA' from 'integrated' before you run the monocle or you will get an error.

[@yongdongwu8472]

ZLJ 3周前 Thanks a lot! you are so excellent !!and i have a question,the s5" filter(status=="OK")",what's difference between “OK" and "failure"，i cannot figure out the meaning and difference.

[@Bioinformagician]

graph_test() results outputs status "FAIL" for those genes whose p values, Moran I and Moran test statistic values are NA

[@jasminlakia]

Would you consider a tutorial for TI with slingshot? that would be so helpful for me!

[@sreejas1302]

Hi, thank you so much for this video. How to do trajectory analysis in diffusion maps instead of principal components?

[@tizianatocci1176]

Hey, your videos are truly amazing! Thank you so much! I have a question concerning point 3 "# 3 Learn trajectory graph -----------". I obtain a trajectory different from the one that you obtain in the video and the same happens if I copy your code from github. Do you have any idea of what I am doing wrong? Thank you again!! This is the error I get when I run the plot of the third point: Warning message: ggrepel: 7 unlabeled data points (too many overlaps). Consider increasing max.overlaps

[@adampassman]

I have hundred of laser-microdissected samples, is there a pipeline for trajectory analysis for this kind of data? Do you think that I could just format my libraries as if they were single cells instead of single microdissections and trick it into working?

[@erin8360]

Thank you for sharing this! Is there a way to compare the trajectories between two conditions/treatments?

[@Bioinformagician]

You could build separate trajectories for each condition/treatment, get top 50 or 100 genes that change expression over pseudotime from both groups and compare them.

[@asadpasha8084]

Hello Miss, i am unable to download the data, the website is not opening

[@aimanbarki9968]

link for data is not working :(.. Can you please cite another source

[@Bioinformagician]

It seems to be working fine for me. Here's the link to the associated publication: academic.oup.com/nsr/article/8/3/nwaa180/5896476#267860186 They have their data deposited in NCBI GEO - GSE137864 and GSE149938

[@zeeman5007]

@@Bioinformagician Hi. Love these videos! The link to the data is not working for me as well. I am able to get the expression matrix file from NCBI GEO. But having difficulty finding the gene and cell metadata files. Any assistance would be greatly appreciated! :)

[@Bioinformagician]

@@zeeman5007 I have uploaded these files here: drive.google.com/file/d/1CJ9VSrUCoqPsUI1jrdm2nrLRawI04xZ1/view?usp=sharing

[@eramartinez8790]

Thanks for the video you are saving my PhD, but I have the same problem, I can´t install SeuratWrappers, this is the Warning in install.packages :package ‘SeuratWrappers’ is not available for this version of R , I tried to install from but devtools::install_github("satijalab/seurat-wrappers") and remotes::install_github('satijalab/seurat-wrappers')nothing, any recommendation? I have R version 4.2.2

[@TheVillka]

Hi, I had the same problem and installed r.utils. Then it worked with SeuratWrappers

[@hsinchiehwu4397]

Thank you for all the Rstudio vedios. They're so clear and helpful !!!! I've failed to instal SeuratWrappers (Mac Monterey 12.5, Rstudio IDE 4.2) using remotes or devtools; troubleshooting follow the github site was also nothing worked. Any suggestion would be greatly appreciated !!! Many thanks,

[@Bioinformagician]

What is the error you see?

[@mmhmmh8545]

I had the same problem, I was able to install SeuratWrappers after installing R.utils package install.packages("R.utils"). try doing this and see.

[@张凯-z4w]

Sorry,my monocle3 cannot find funtion as.cell_data_set(). Can you help me?

[@merrinleong1388]

Error in as.cell_data_set(b.seu) : could not find function "as.cell_data_set"

[@Bioinformagician]

Have you loaded SeuratWrappers library before running that command?

[@moinicholas3828]

Hi, Following yout script I find that mitopercent is always 0, any ideas why this might be?

[@moinicholas3828]

Is it that the sample size is small enough that there are actually 0 mitochondrial genes?

[@moinicholas3828]

I also find that seu.obj.filtered@metadata has the same number of rows as the orignal, is this right?

[@MsZhang666]

The voice and the video do not match on 25:00-33:00 min.

[@marinafernandez6778]

Hi, Great video! thanks! I got this error any idea how to solve it? seu.obj$mitopercent

[@xinyuqu4407]

I have the same issue. I use the AddMetaData function from Seurat to solve the problem. rownames(metadata)

[@soburali5430]

@@xinyuqu4407 It worked, Thank you so much!

[@felixoluwasegun8567]

@@xinyuqu4407 Hello, I have tried this solution multiple times; it still shows the same error. is there anything else I can do?

[@burakkahveci4123]

Hi, thank you for this great video, but I have a problem. When I install monocle3, I get a error. ** byte-compile and prepare package for lazy loading Hata: package or namespace load failed for 'SummarizedExperiment' in library.dynam(lib, package, package.lib): DLL 'DelayedArray' not found: maybe not installed for this architecture? Ek olarak: Warning messages: 1: package 'matrixStats' was built under R version 4.1.3 2: package 'GenomicRanges' was built under R version 4.1.2 3: package 'S4Vectors' was built under R version 4.1.3 4: package 'GenomeInfoDb' was built under R version 4.1.2 Çalıştırma durduruldu ERROR: lazy loading failed for package 'monocle3' * removing 'C:/Users/Burak/OneDrive - Dokuz Eylül Üniversitesi/Belgeler/R/win-library/4.1/monocle3' Warning message: In i.p(...) : installation of package ‘C:/Users/Burak/AppData/Local/Temp/RtmpC69TpC/file28646e1844d/monocle3_1.2.9.tar.gz’ had non-zero exit status