This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Cell-type-specific markers across datasets

Hi,

I'm working with two snRNA-seq datasets and I need some advice on the best strategy to compare cell-type-specific markers across datasets.

A previous analysis, generated with 10X, focusing on a large part of the brainstem, allowed us to identify specific gene markers of a population of neurons amoug other populations using FindAllMarkers in Seurat. We now have a new sequenced dataset generated with FLASH-seq, restricted to one part of the previous analysis. Since both datasets originate from overlapping regions, one of our questions is whether the genetic signature identified in the first dataset can be recovered in this new one. Therefore, validating the new method we are testing in addition to the quality control analysis. Moreover, since we just wanted to test FLASH-seq method, we drove the analysis on a small batch (32 cells). Because of that, we cannot use the following Seurat fonctions : UMAP, Clustering and FindAllMarkers.

Using FindAllMarkers on the 10X dataset gives us a cluster marker list that reflect internal differences within the dataset. Therefore, comparing this list to the top AverageExpression genes of the second dataset doesn't seem biologically meaningful to study shared biological identity across datasets. Due to differences in dataset size and maybe method as well, I doubt that we can merge both datasets.

Has anyone dealt with a similar cross-dataset cell-type validation question?

Thank you !

gene snrnaseq seurat findallmarkers

2 answers

From my experience working across datasets generated with different technologies, the safest approach is usually to analyze each dataset independently first, then compare the results to identify the consensus signal rather than trying to force a direct merge from the beginning.

For your case, I would separate the problem into two steps.

  1. Assign cell-type identity consistently in both datasets

First, try to label the cells in each dataset using known markers for the neuronal populations you care about. You can find some helpful markers here: https://omnibusx.com/celltypes

If you have already identified a marker signature from the 10X dataset, that signature can be used as a reference when examining the FLASH-seq cells. Of course, some genes may be missing or weakly detected in the FLASH-seq dataset because of the technology, coverage, or the small number of cells, so perfect agreement is not necessarily expected. At that point, the interpretation becomes biological versus technical, and that needs to be judged carefully.

  1. Identify markers for the target population independently in each dataset, then intersect them

The shared genes are usually the most convincing evidence that the same biological identity is being recovered across datasets.

With only 32 cells in the FLASH-seq dataset, I would be cautious about relying on clustering-based workflows. In this setting, a targeted validation of a known signature is often more meaningful than trying to rediscover the full structure de novo. So rather than asking whether the FLASH-seq dataset reproduces the entire clustering landscape, I would ask whether it still captures the expected marker program of the neuronal population you identified in the larger 10X dataset.

I hope this helps.

I would start with just mapping your FLASH-seq dataset onto your original 10X dataset, which is pretty straightforward with Seurat. Should let you annotate cell types and visualize them side-by-side in a clear way.

Log in to answer this question.