Differentially expressed genes analysis in Seurat

Bioinformatics Asked on September 30, 2021

For the differentially expressed genes analysis, is it possible to check for DEGs based on the levels already identified in the object? For example, my dataset contains cells from 11 subjects, which are marked as so (orig.ident). I want to analyze the differentially expressed markers between them. Will it be accurate to do so? Or should I follow the clustering analysis in Seurat? I want to see the difference in genes between the cells according to the subjects but I’m afraid that may not be the right way to do so.

  head([email protected])
         orig.ident nCount_RNA nFeature_RNA percent.mito
 BC01_02       BC01   998944.0         9189 0.2387108038
 BC01_03       BC01   999696.4         9548 0.2925345423
 BC01_04       BC01   999057.5         7440 0.0009893724

Is clustering analysis absolutely necessary to see the DEGs between them? or can I directly move to FindMarkers based on the levels (orig.ident)?

One Answer

I'm not quite sure what you are asking, but you can use anything in the metadata as a grouping with FindMarkers.

DefaultAssay(combined.all) <- "RNA"
Idents(combined.all) <- "orig.ident"

for (sample in unique(Idents(combined.all))) {
  clustermarkers <- FindMarkers(combined.all, ident.1 = sample)
  print(paste("Cluster", sample))
  print(head(clustermarkers, n = 10))

Answered by swbarnes2 on September 30, 2021

Add your own answers!

Related Questions

Viral Metagenomics

1  Asked on November 1, 2020 by l-r-joshi


Issues with AutoDock Vina

0  Asked on October 18, 2020 by ibio_rep1


Swapping to effect increasing allele in case/control studies

0  Asked on October 9, 2020 by dale-handley


Error When Using biocLite as an installer in rpy2 python library

1  Asked on September 27, 2020 by abiologist


Seurat DE t.test

1  Asked on August 11, 2020 by vdu12345


Number of reactions per metabolic pathway

0  Asked on August 11, 2020 by mmphysics


Find all the bases for given reference position

0  Asked on August 8, 2020 by diesel__100


Convert VCF to genotype table

1  Asked on July 30, 2020 by snowflake


Ask a Question

Get help from others!

© 2023 All rights reserved. Sites we Love: PCI Database, MenuIva, UKBizDB, Menu Kuliner, Sharing RPP, SolveDir